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PREFACE 


This  technical  report  consists  of  publications,  manuscripts,  and 

supplemental  information  related  to  the  scattering  of  light  by  bubbles 

in  liquids  and  in  glass.  The  emphasis  of  the  report  is  on  the  scattering 

by  air  bubbles  in  water  into  the  angular  regions  where  diffraction  makes 

essential  corrections  to  elementary  ray  optics.  For  bubbles  these  regions 

are  backward  scattering,  forward  scattering,  and  critical  angle  scattering. 

Most  of  the  previous  literature  on  light  scattering  is  concerned  with 

drop-like  objects  where  the  refractive  index  of  the  scatterer  exceeds 

that  of  the  surroundings.  The  present  results  should  be  useful  for  the 

1  2 

optical  characterization  of  microbubbles  at  sea  and  in  water  tunnels  , 
or  for  the  prediction  of  the  optical  properties  of  bubbly  media.  They 
should  also  be  useful  for  1 aser-Doppler  anemometry  with  microbubbles  and 
the  detection  of  microbubbles  in  glass. 

Certain  phenomena  described  here  should  be  present  for  other  spherical 
or  nearly  spherical  scatterers  where  the  refractive  index  is  less  than  that 
of  the  surroundings,  for  example:  air  bubbles  in  ice,  water  inclusions  in 
fused  quartz  (such  as  those  present  in  varieties  of  opal  glass),  ice  spheres 
in  water,  and  superheated  drops  in  a  host  liquid.  The  coarse  structure  in 
the  critical  scattering  region  can  be  present  even  if  the  scatterer  is  not 
spherical . 

The  first  paper  summarizes  experiments,  models,  and  results  of  Mie 
computations  for  both  backward  and  critical -angle  scattering.  Theoretical 
results  presented  indicate  that  axial  focussing  enhances  the  scattering  by 
an  air  bubble  in  water  into  the  near  forward  and  near  backward  regions.  This 
enhancement  for  the  forward  region,  (forward  optical  glory),  is  more  signifi- 
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cant  for  depolarized  than  for  polarized  scattering  since  the  latter  is 
dominated  by  ordinary  forward  diffraction.  At  present  the  experimental 
observations  of  axial  focussing  from  single  bubbles  are  limited  to  the 
backward  glory  of  air  bubbles  1‘n  a  viscous  polymer  liquid.  (These  are 
described  in  detail  in  the  seventh  paper).  The  first  paper  also  discusses 
the  scattering  efficiency  of  bubbles  and  the  relevance  of  far-field  compu¬ 
tations  to  near-field  observations  of  bubbles. 

The  second  paper  Describes  the  first  observations  of,  and  model  for, 
the  scattering  in  the  critical  region  which  is  near  83  degrees  f or  an  air 
bubble  in  water.  This  paper  is  included  here  for  completeness  though  it 
was  not  supported  by  the  present  ONR  contract.  It  was  supported  in  part 
by  ONR  Contract  N00014-76-C-0527  (R.  E.  Apfel  principal  investigator). 

The  third  paper  describes  an  improved  model  (which  includes  both 
diffraction  and  interference)  for  the  coarse  structure  near  the  critical 
scattering  region. 

The  fourth  paper  is  the  first  to  be  published  which  correctly  gives 
the  results  of  Mie  theory  for  angular  structure  in  the  scattering  by  bubbles. 

This  paper  substantiates  the  model  of  the  coarse  structure  developed  in  the 
third  paper.  The  paper  was  excerpted  from  the  M;  ster's  Degree  Thesis  of 
D.  L.  Kingsbury.  A  supplement  attached  here  gives  examples  which  were  not 
published  for  reasons  of  brevity. 

The  fifth  paper  gives  Mie  theory  and  model  results  for  air  bubbles  in 
fused  silica  glass.  The  paper  is  the  first  to  discuss  the  Brewster  scattering 
angle  for  bubbles  and  the  first  to  give  valid  Mie  results  for  bubbles  over 
the  entire  (0  to  180  degrees)  range  of  scattering  angles. 

The  sixth  paper  describes  computer  codes  used  for  some  of  the  Mie  and 
model  computations.  The  codes  shown  here  have  been  adapted  fov’  use  on  an  HP  1000 
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minicomputer  system.  This  paper  is  excerpted  from  the  Master's  Degree 

Thesis  of  D.  L.  Kingsbury.  For  a  more  complete  description  of  the  Mie 

3 

scattering  algorithm,  consult  the  paper  by  Wiscombe  which  is  the  basis 
of  these  codes. 

The  seventh  paper  describes  the  first  observations  of  backscatteri ng 
from  air  bubbles  in  liquids.  It  also  describes  a  physical -optics  approx¬ 
imation  for  the  backward  axial  focussing  of  scattered  light.  This  model 

shows  that  the  intensity  of  individual  axially  focussed  rays  are  propor- 
3 

tional  to  a  ,  while  the  simply  reflected  scattering  is  proportional  to 
2 

a  where  a  is  the  bubble  radius.  (The  total  backscattering  is  not 

3 

simply  proportional  to  a  due  to  the  interference  of  various  paths; 

see  paper  1).  Polarization  and  quasi -periodic  properties  of  the  glory  are 

discussed. 

The  approximations  and  physical  models  described  here  for  bubbles 
should  also  be  useful  for  certain  cases  in  the  scattering  of  sound  from 

4 

fluid  or  elastic  spheres.  This  has  been  verified  by  Marston  and  Kingsbury 

for  scattering  in  the  critical  region  from  fluid  spheres.  Backward  axial 

focussing  or  "acoustic  glory"  has  also  been  proposed  for  elastic  spheres 
5 

in  water. 

Care  should  be  taken  when  applying  the  results  of  this  report  to 
light  scattering  by  mircobubbles  in  seawater.  Our  computations  assume 
that  the  bubbles  are  spherical  and  that  the  refractive  index  of  the 
scatterer  is  homogeneous.  There  is  a  natural  tendancy  for  the  gas  in 
microbubbles  to  disolve  into  the  surrounding  liquid. ^  There  is  some 
recent  evidence  that  microbubbles  in  seawater  can  be  stabilized  by 
substances  sorbed  onto  their  surface.^  The  presences  of  a  sorbed  film 
or  small  deviations  from  sphericity  may  alter  the  details  of  the 
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scattering  patterns  described  here,  but  it  should  not  obliterate  the 

existence  of  axial  focusing  and  the  coarse  and  fine  structures. 

Philip  L.  Marston 
Principal  Investigator 
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Paper  No.  1 

Light  scattering  by  bubbles  in  liquids:  Mie  theory,  physical -optics 
approximations,  and  experiments  (P.  L.  Marston,  D.  S.  Langley, 
and  D.L.  Kingsbury)  to  be  published  in  the  journal  Applied 
Scientific  Research  in  the  proceedinos  of  the  IUTAM  Sumposium 
on  the  Mechanics  and  Physics  of  Bubbles  in  Liquids.  This 
paper  was  presented  at  the  Symposium  (June,  1981). 


Abstract 

Angular  structures  in  the  far-field  scattering  from  bubbles  are 
observed  and  modeled.  Mie  theory  supports  a  model  of  diffraction  arid 
interference  near  the  critical  scattering  angle.  A  new  expression  for 
the  angular  spacing  of  fine  structure  is  derived.  Photographs  of  scat¬ 
tering  show  some  of  the  predicted  features.  Application  of  these  struc¬ 
tures  to  bubble  sizing  and  detection  are  summarized  and  the  theoretical 
extinction  coefficient  in  water  is  plotted. 

Mie  computations  for  bubbles  in  water  also  reveal  backward  and  for¬ 
ward  glory  effects.  These  are  partially  manifested  as  cross-polarized 
scattering.  Observed  scattering  from  bubbles  in  ^he  near  backward  di¬ 
rection  is  found  to  have  a  strong  cross-polarized  component. 
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1.  Introduction 

Equations  for  the  scattering  of  plane  electromagnetic  waves  by  a  di¬ 
electric  sphere  were  given  by  Mie  [1]  in  1908  and  the  resulting  features 
of  the  angular  scattering  pattern  of  drops  are  well  known  [2-4].  The  Mie 
solution,  though  exact,  does  not  give  insight  into  either  the  scattering 
process,  or  changes  in  the  pattern  resulting  from  changes  in  shape,  re¬ 
fractive  index,  c.  profile  of  the  incident  wave  fronts.  Reviews  of  the 
literature  on  light  scattering  [2-4]  reveal  a  paucity  of  information  about 
the  scattering  pattern  of  bubbles  in  liquids  where  the  refractive  index 

of  the  scatterer  n.  is  less  than  that  of  the  surroundings  n  .  Conse- 

l  o 

quently  we  have  begun  [5-10]  a  systematic  study  of  the  scattering  of  light 
by  gas  bubbles  in  liquids.  Aspects  of  the  study  are:  (1)  the  computation 
of  Mie  scattering;  (2)  the  development  of  simple  physical  models  which  give 
insight;  and  (3)  observations  of  features  in  the  scattering  which  differ 
significantly  from  both  the  scattering  by  drops  and  the  scattering  pre¬ 
dicted  by  geometrical  optics  [11].  In  this  paper  we  summarize  the  main 
features  of  the  scattering  with  an  emphasis  on  the  critical  [5-7]  and  back- 
scatter  [10]  regions.  New  experiments  and  applications  will  be  described. 

These  are  the  first  detailed  observations  of  scattering  by  bubbles. 

Mie's  solution  [1-3]  to  the  problem  of  the  scattering  efficiency  and 
pattern  of  a  dielectric  sphere  is  a  function  of  the  ratio  m  -  n^/nQ  .  It 
is  usually  expressed  as  a  function  of  x  ::  ka  -  2 ir.-i/A  where  a  is  the  sphere 
radius,  k  and  A^  are  the  wavenumber  and  wavelength  in  the  outer  dielectric. 
Optical  sources  are  typically  characterized  by  the  wavelength  in  a  vacuum 
we  note  that  k  =  2nno/Av.  For  an  air  filled  bubble,  n .  =  1.00029 
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and  it  is  usually  a  good  approximation  to  take  in  =  1/n  .  Most  of  this 
paper  will  deal  with  scattering  by  spherical  bubbles  with  plane  incident 
waves.  We  will  gain  some  insight  however,  into  which  features  of  the 
scattering  pattern  should  be  sensitive  to  deviations  from  sphericity. 

2.  Scattering  Efficiency 

The  efficiency  factor  Q  is  the  ratio  of  the  total  scattering 

SCd 

2 

cross  section  to  the  geometric  cross  section  nu  .  When  n.  is  real 

valued,  Q  also  gives  tfie  extinction  efficiency  [2-4,12].  Figure  1 
sea 

compares  Q  for  drops  of  water  in  air  with  m  -  4/3  with  that  for  a 

bubble  in  water  with  in  =  3/4.  The  computations  were  performed  using  a 

slightly  modified  version  [9]  of  Wiscombe's  MIEVO  Mie  scattering  algorithm 

[13].  A  table  of  Q  for  bubbles  exists  [14]  which  is  consistent  with 

Fig.  1.  For  a  fixed  Ay,  and  a  given  value  of  x  ,  the  drop's  radius  is 

larger  than  the  bubble's  radius  by  a  factor  of  4/3.  For  the  case  of 

light  from  a  He-Ne  laser,  \  =  0.6328  pm;  for  bubbles  in  water,  x  of  10, 

100,  1000,  and  10  000  give  a  of  0.75,  7.55,  75.5  ,  and  755  pin,  respectively. 

The  salient  feature  of  Fig.  1  is  that  for  drops  Q  exhibits  a 

sed 

fine  "ripple"  structure  [3,12]  but  that  our  calculations  of  Q  for 

sea 

bubbles  do  not  reveal  a  ripple  structure.  For  both  drops  and  bubbles,  Q 

SCd 

has  a  broad  undulation  with  a  quasi-period  Ax  ~  ir/|m  -  1|.  This  approxi¬ 
mation,  which  has  been  derived  from  the  theory  of  "anomalous  diffraction" 
[2,3],  appears  to  be  useful  for  both  bubbles  and  drops.  The  ripple  struc¬ 
ture  present  for  drops  is  due  to  optical  resonances  [12]  which  are  attributed, 
ir,  part,  to  internal  surface  waves.  The  absence  of  such  structure  for 
bubbles  is  probably  because  m<l  does  not  favor  the  entrapment  of  internal 
surface  waves. 
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3.  Critical  and  Brewster  Angle  Scattering 


For  scattering  by  drops,  diffraction  is  important  for  the  description 
of  the  forward,  backward ,  and  rainbow  regions  [2,3].  For  bubbles,  there  is 
no  longer  a  rainbcw;  however,  a  new  region  appears  known  as  the  critical 
scattering  region  [5],  Diffraction  is  important  in  this  region  because  of 
an  abrupt  change  in  the  amplitude  of  the  reflected  wave  as  the  local  angle 
of  incidence  0  changes  from  6  <  0c  for  small  impact  parameters  to  0  >  0c 
for  large  ones.  Here  0c  =  arcsin(in'^)  which  is  the  critical  angle  for  a 
plane  surface.  Figure  2  illustrates  several  ray  paths  which  lead  to  a 
scattering  angle  <J>  (the  deviation  from  the  forward  direction)  of  50°. 


The  number  on  the  left  specifies  the  number  of  internal  chords 
denote  the  angle  of  incidence  and  refraction  of  the  pth  ray. 


p;  0  and  n 
P  ''P 

The  reflected 


ray  (which  has  p  =  0  )  has  a  scattering  angle  <J>  =  tt  -  26_.  For  an  air 

u 

bubble  in  water  with  m  =  3/4,  the  critical  scattering  angle  is  <t  =  it  -  20  = 

c  c 

82.82°.  Geometric  optics  [5,  6,  11]  predicts  that  |dl./d<|>|  ■*  «  as  ch  ap- 

J 

proachcs  r  from  above  4  .  Here  I.  is  the  normalized  scattered  intensity 
c  c  j 

defined  as  follows:  the  actual  j-polarized  intensity  at  a  distance  R  >>  a 


from  the  bubble's  center  is  the  incident  j-polarized  intensity  multiplied  by 
2 

I,  (a/R)  /4.  For  the  electric  vector  perpendicular  to  the  scattering  plane, 

J 

j  =  1,  for  the  parallel  case,  j  =  2.  This  normalization  is  appropriate  for 


bubbles  since  geometric  optics  predicts  tnat  if  the  intensity  of  the  p  =  0 

reflection  could  be  considered  by  itself  I .  ($  <  <j>  )  =  1  due  to  total  reflection. 

J  c 

Our  observations  [5],  model  [5,  6],  and  Mie  computations  [7-9]  demonstrate 

that  instead  of  a  divergence  of  |dl./d<j>|  ,  the  rise  in  I.  is  spread  out  over 

J  J 

the  region  |<|>-4>c|  -  SI  where  Cl,  =  arcsin  [0.8  (l-m2)_Js(X  /a)’5]  *  144x_iideg 
for  water.  The  model  makes  use  of  a  physical  optics  approximation  which 
is  to  (a)  use  ray  optics 
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along  with  the  reflection  coefficients  of  a  plane  surface  to  compute  the 
complex  amplitudes  of  virtual  waves  which  simulate  the  scattering,  and  (b) 
use  Fraunhofer's  approximation  to  compute  the  diffraction  of  the  virtual 
wave  to  the  far  field  where  R  >>  xa.  Step  (b)  yields  the  improvements 
over  the  geometrical  predictions  [11].  This  procedure  is  analogous  to 
Airy's  model  of  diffraction  neat'  the  rainbow  [2]  .  Additional  approximations 
used  in  the  model  are:  only  the  two  most  intense  virtual  waves  for  <f.  -  cp 
are  included  (p  =  0  and  1);  and  the  divergence  of  the  derivative  of  the 
reflectivity  is  simulated  by  truncating  the  p  =  0  reflection  when  °o  <  °c  ' 

To  compare  the  model  results  with  Mie  theory,*  we  consider  examples  of 
bubbles  in  water  which  complement  those  previously  published  [7].  Figure  3 
shows  both  theories  for  x  =  100  and  j  =  1.  It  shows  that  there  is  a  coarse 
structure  in  the  Mie  scattering  for  <p  <  ^  which  has  a  quasi-period  <.  SI  ~  14.4". 
This  structure  is  described  by  the  model  excent  in  the  near  forward  direction 
(4>  £  20°). 

The  Mie  result  also  has  a  superposed  fine  structure  where  the  magnitude 

of  the  quasi-period  <  \Ja  rad.  =  360°x~^  .  Near  ,  this  fine  structure 

0  c 

arises  primarily  from  the  interference  of  the  p  =  0  wave  with  the  wave  due 
to  the  p  ~  2'  ray  in  Fig.  2.  Its  quasi-period  .nay  be  estimated  from  the 
lateral  separation  bQ  +  b?l  of  virtual  point  sources  which  would  simulate  the 
scattering  at  a  given  ft  .  As  shown  in  Fig.  2,  this  separation  is  the  sum  of 
the  impact  parameters  for  these  rays.  Standard  relations  for  the  far- 
field  interference  applied  to  these  sources  gives: 

nf(<|>)  =  arcsinUo/(bQ  +  , ) ]  (1) 


*The  conversion  from  the  Mie  amplitudes  3-  [2,3  13]  to  the 

UlSjl/x)2. 


I  .  is 
J 
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where  b^/a  =  sinO^  »  °q  =  U  ~  t)/2*  <t>  =  Tr  +  2  ( 0  ^ ,  -  2p^ , )  an1!  from  Snell's 

law,  msinpp  =  sinO^  .  Lq.  (1)  is  only  an  approximation  because  it  fails 
to  include  other  (such  as  p  =  1)  virtual  waves  and  it  does  not  completely 
account  for  the  longitudinal  spacing  of  virtual  sources  and  for  how  those 
sources  vary  with  >|>.  With  cf  =  <{>  and  m  =  3/4,  b^  +  b^,  ~  a/0.825  and 
fir  -  arcsin  (5.18  x  ' ) •  This  gives  2.97°  at  x  •-  100  which  is  in  reasonable 
agreement  wi th  fig.  3 . 

Figure  4  shows  the  Mie  and  phys ica I -optics  results  for  x  -  10  000  and 
j  =  1.  The  fine  structure  quasi-period  is  greatly  reduced;  Fq.  (1)  gives 
ftfUc)  ~  0.0297°  which  agrees  with  the  Mie  result  of  0.03°.  The  amplitudes 
of  the  coarse  undulations  (n  ~  1.4°)  decrease  slightly  with  decreasing  <|i  until 
4)  ~  75°.  A  graph  of  1^  for  |>  •:  75°  shows  coarse  undulations  increasing 
in  amplitude  with  decreasing  4.  The  explanation  is  that  for  x  =  10  000, 
with  <J>  %  75°  the  coarse  structure  is  primarily  due  to  diffraction  of  the 
p  =  0  wave  [5]  but  for  <j>  ^  75°  it  is  primarily  from  the  interference  of  the 
p  =  0  wave  with  p  >  0  waves  [6].  Plots  [7]  of  for  this  x  also  show 
thi s.  transi t ion  and  that  fine  structure  is  significantly  weaker  when  j  -  2. 

Figure  3  does  not  show  the  transition  because  the  diffraction  region  (•}>  -  4  ,<  Cl) 

overlaps  the  region  where  interference  with  the  p  -  1  wave  is  significant. 

As  0g  approaches  tfie  Brewster  angle,  0^  -  arctan(m),  the  reflectivity 

of  the  p  =  0,  j  =  2  polarized  ray  vanishes  [2].  Consequently,  according  to 

geometric  optics,  there  is  no  contribution  to  from  this  ray  at  the  Brewster 

scattering  angle  4. ..  =  it  -  2up  =  2  arctan(m  ^ ) .  Because  the  reflectivity  varies 

t)  b 

slowly  near  <p ^ ,  diffraction  is  less  important  than  near  4  .  A  coarse  minimum 
is  evident  in  the  Mie  near  4  for  bubbles  in  glass  [8]  and  for  bubbles 
in  water  where  4^  ~  105.3°.  When  x  drops  below  5,  its  location  shifts  toward 
90°;  the  scattering  pattern  approaches  of  o  dipole  radiator  predicted  by 

Rayleigh  scattering  theory  [2,  3,  8]. 


13 


6 


4.  Observations  of  Critical  Angle  Scattering 

Previous  observations  [5]  were  limited  to  the  <j>  region  where  the 
coarse  undulations  decreased  with  decreasing  .  We  have  a  new  apparatus 
similar  to  that  in  [5]  except  that  it  permits  observations  with  4)  down  to 
70°.  As  in  [5],  bubbles  were  attached  to  a  vertical  needle  in  distilled 
water.  They  were  illuminated  with  a  plane  wave  from  a  He  -  Tie  laser  with 
-  632.8  mil.  They  were  photographed  via  a  window  with  a  camera  focused 
on  infinity  to  place  the  film  in  the  far  field.  Because  of  a  technical  problem 
we  have  not  obtained  precise  direct  measurements  of  ■:  in  the  (horizontal) 
scattering  plane  along  with  the  photographs,  however,  we  have  demonstrated 
that  )  from  Eq.  (1)  gives  an  a  which  is  consistent  with  the  observed 

and  modeled  coarse  structure.  Thus  the  n  values  quoted  in  F'g,  5  were 
determined  from  the  coarse  and  fine  structures  and  not.  from  direct  observations 
Figure  5  serves  more  to  illustrate  phenomena  described  in  Sec.  3  than  to  rig¬ 
orously  test  the  models  with  real  bubbles.  The  angle  scale  is  for  ?.  t}> 

change  in  water.  The  horizontal  displacement  is  not  exactly  linear  in  <j>  due 
to  refractive  corrections  at  the  water-window-aii  interfaces.  In  Fig.  5(a), 
the  broad  bright  region  on  the  left  is  the  first  coarse  maximum  with  <},  <  ,j)  . 

On  the  right,  the  coarse  undulations  increase  in  amplitude  as  <J>  decreases, 
which  is  the  modeled  behavior.  Figure  5(b)  is  for  a  smaller  bubble  (verified 
Vfith  direct  observations).  The  coarse  undulations  are  spread  out  and  the  fine 
structure  is  clearly  visible  with  -•  0.07°. 

5.  Backward  and  Forward  Glory 

Backscattering  from  drops  is  known  to  be  enhanced  by  the  axial  focusing 
of  certain  rays  [2,  3],  A  complete  description  of  this  "glory"  for  water  drops 
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is  complicated  due  to  the  ner  ssity  of  including  surface  waves  [15].  We 

have  observed  and  modeled  backscattering  from  air  bubbles  [10]  and  find  that 

some  aspects  can  be  explained  with  a  physical -optics  approximation  that  does 

not  include  surface  waves,  for  the  following  discussion,  it  is  convenient 

to  define  y  =  it  -  <{>;  y  -  0  corresponds  to  exact  backscattering . 

Computed  Mie  intensities  show  that  backscattering  from  bubbles  may  be 

significantly  larger  than  1 ha t  expected  from  reflection  of  the  p  =  0  ray. 

For  y  -  0,  the  plane  of  scattering  is  no  longer  defined  and  1^  =  =  I. 

! igu'c  6  shows  several  broad  peaks  in  1  which  ar>  significantly  larger  than 

the  normalized  intensity  of  the  p  =  0  reflection  which,  by  itself,  is  (m  -l)'/ 

2 

(m  +  1 )  =  0.02.  Other  axial  (e.g.  p  2)  rays  are  too  weak  to  explain  the 

magni tude  of  I . 

For  polarized  incident  light,  scattering  near  y  -  0  may  be  described  in 

part  with  (non-cross )  polarized  1^  and  cross-polarized  1  ^  normalized  in- 
(2) 

tensities.  For  r  ,  the  plane  of  polarization  is  rotated  by  90°  from  incident 

polarization;  symmetry  gives  =  0)  -  0  while  l^(y  -  0)  =  I.(y  =  0). 

in 

For  y  ^  0,  the  I  depend  on  both  y  and  the  angle  f,  which  the  scattering 
plane  make,  with  the  incident  electric  vector.*  Figure  7(a)  illustrates  the 
computed  scattering.  Geometric  models  of  l^(y  ~  0)  predict  that  the  con¬ 
tributions  of  off-axis  (p  -  3  and  4)  rays  diverge  as  y  -►  0  because  of  a  factor 
which  accounts  for  focusing.  This  divergence  is  present  in  previous  ray  models 
[11],  but  was  not  discussed;  it  is  also  present  for  drops  [2]  for  certain  m  • 
Because  of  this  divergence,  diffraction  provides  an  essential  correction  to 
ray  optics. 


*The  conversions  from  the  complex  Mie  amplitudes  S.  to  the  I^(y«10°,F, 

.2,2  ■  .(2)  .  .2,  2 


are  I 


(1 


5  -  S  00^2',  i  / x  and  I 


S  sin2f, !  /x  where  S' 


S-j  1  5^.  Evidently  the  dependence  of  Eg.  (2)  is  of  the  correct  type. 
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The  scattering  due  to  the  p  -  3  ray  in  Fig.  /(a)  (insert)  is  a  back- 
facing  toroidal  wave  front  which  appears  to  originate  at  a  virtual  ring¬ 
like  source  known  as  the  focal  circle  in  the  analogous  [2]  p  -  ?.  scattering 
from  drops.  This  source  at  point  F  is  ringlike  because  the  figure  may  be 
rotated  around  the  optic  axis  through  the  center  C.  For  largo  x,  the 
stationary-phase  approximation  of  far-field  diffraction  integrals  gives  the 
following  proportionalities  for  otf-axis  (p>2)  contributions  to  the  1^ 
taken  separately: 

(  1  ^  9 

I(p)  cx  xt(c1  +  c2)  00(u)  +  (c]  -  c2)  d,,(u)  C0S2/JL  I 


l{pj  x[(Cl  -  c2)  J2(j)  sinZ;.]? 


where  the  c.  are  combined  Fresnel  reflectivities  and  transmissivities  of 


the  j-polarized  fields,  u  =  kb^siny  ,  and  bp 


is  the  impact  parameter 


of  the  exactly  backscettered  pth  ray.  Because  l|  j  <•>  x  ,  the  unnormal i zed 

intensity  is  proportional  to  k./  while  that  for  the  p  =  0  ray  is  proportional 
2 

to  ,/  .  Consequently,  off-axis  waves  are  dominant  when  <i  is  sufficiently 

large.  It  is  evident  from  Fig.  b  and  7(a)  that  the  modeling  of  rl  ^  may  re¬ 
quire  summation  of  electric  fields  from  several  virtual  sources.  (In  physical 
optics,  intensities  do  not  add  but  fields  do.)  Interference  of  various  ring¬ 
like  and  axial  virtual  sources  depends  on  x  and  could  lead  to  the  broad  struc¬ 
tures  in  Fig.  6.  Detailed  computation  of  ij^j  (y  =  0)  give  magnitudes  suf¬ 
ficient  to  explain  the  Hie  scattering.  Due  to  the  Bessel  functions  in  Fq.  (2), 
the  actual  scattering  should  be  peaked  at  or  near  y  0. 

We  have  observed  backscattering  from  single  air  bubbles  in  a  dimethyl- 
si  loxane-polymcr  liquid  which  were  nearly  immobilized  by  viscosity.  In  these 
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experiments ,  m'1  =  1.403,  .\  =  632.8  nm,  a  -  0.3  -  0.8  mm  and-  x  -  4000  - 

11  000.  The  far-field  cross-polarised  intensity  had  a  dependence  on  y  and 
£  similar  to  that  predicted  by  Eq.  (2b)  with  p  -  3.  Equation  (2b)  predicts 
that  when  both  siny  y  and  u  :•>  1,  the  minima  in  1^  should  be  spaced 
by  Ay  ~  it  (kb^T1  radians  where  b.'{  -  0.44/a.  The  prediction  correctly 
described  the  observed  f  -  2  scattering  which  is  apparently  dominate^  here  by 
the  p  =  3  virtual  source.  The  focal  circle  was  viewed  by  focusing  the  camera  on 
the  bubble.  The£=  1  scatteriny  was  not  dominated  by  a  single  class  of  rays. 
Cross-polarized  scattering  from  a  polydi spersicn  of  bubbles  in  water  ha s  also 
been  seen. 

Figure  7(b)  shows  an  enhancement  of  near-forward  cross-polarized  scatteriny 
due  to  axial  focusing.  Focal  circles  due  to  the  p  =  ?.  and  3  rays  largely 
contribute  to  the  C  ~  2  scattering.  The  forward  y  =  1  scattering  is  domi¬ 
nated  by  ordinary  diffraction  [2'J  when  x  is  large.  Hie  theory  gives  I^~  10^ 
when  ^  =  0  for  the  x  and  m  of  Fig.  7.  Forward  optical  nlory  has  also  been 
displayed  in  Mi"  scattering  by  water  drops  [16].  Its  description  is  complicated 
by  surface  waves. 

5.  Discussion 

Figures  3-7  demonstrate  that  for  several  angular  regions,  she  intensity 
exceeds  that  of  geometric  scattering  from  a  perfectly  reflecting  sphere  of  the 
same  a  (for  which  K  =  1).  This  information  should  be  useful  in  optical  de¬ 
vices  which  size  or  detect  buobles  [7,  17,  18].  Due  to  diffraction  near  , 
it  is  preferable  to  detect  scattering  with  .(>  <  vj>  -  n  than  to  observe  it 
with  <J)  =  90°  which  is  the  usual  practice  [1/]  for  bubbles  in  water.  For  an 
unpolarized  source,  the  normalized  total  scat  taring  is  i!^  1  l ,,)/?. 
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Though  this  paper  has  emphasized  the  far-field  |  dependence,  the 
results  are  also  applicable  to  the  imaging  of  bubbles.  For  example,  if  it 
desired  to  resolve  the  virtual  sources  of  the  p  =  0  and  1  rays  in  Fig.  2, 
the  aperture  of  the  imaging  system  should  have  an  angular  width  >  Q  .  It 

*\j 

is  easier  to  resolve  the  p  -  0  and  2'  virtual  sources  since  the  angular 
width  requirements  are  reduced  to  Resolution  of  identifiable  virtual 

sources  reveals  the  size  of  a  bubble. 

Fine  structure  in  the  far-field  scattering,  such  as  that  shown  in  Fig.  5(b), 
arises  from  the  interference  of  widely  spaced  rays.  C)nsequenuy  the  positions 
of  the  maxima  will  be  quite  sensitive  to  changes  in  the  bubble's  shape.  That 
part  of  the  coarse  structure  due  to  the  interference  o\  p  -  0  and  1  rays 
will  also  depend  on  the  shape  but  more  weakly.  As  is  the  case  for  drops  [19], 
these  shape  dependences  may  be  useful  for  detecting  mechanical  resonances. 

Scattered  intensities  for  bubbles  in  water  purported  to  be  from  Mie  theory 
have  beer,  used  in  a  study  of  cavitation  nuclei  [1/].  Comparison  of  our  Hie 
results  and  model  of  the  coarse  structure  (which  are  consistent  for  the  a  in 
question)  with  those  in  [17],  show  that  the  latter  err  significantly.  Plots 
in  [17]  predict  coarse  maxima  at  $  =  100°  for  a  =  5  pm  and  <p  ~  120"'  for 
a  ■  7.5  pm  which  are  not  present  in  our  Mie  computations.  We  cannot  find  any 
physical  justification  for  coaise  maxima  at  these  \  for  x  near  100. 

We  are  grateful  to  W.  J.  Wiscombe  for  providing  the  initial  comouter 
program  from  which  the  program  used  here  was  derived.  This  research  was 
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Committee  and  by  the  Office  of  Naval  Research,  P.  L.  Mars  ton  is  an  Alfred 
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Figure  Captions 


Fig.  1  Scattering  efficiency  from  Mie  theory  as  a  function  of  the  size 
parameter  x  .  The.  solid  curve  is  for  drops  and  the  dashed  curve  is  for 
bubbles . 

Fig.  2  Ray  paths  for  a  bubole  with  in  =  3/4  .  The  number  of  internal 
chords  p  is  given  a  prime  if  the  ray  enters  the  bubble  below  the  center- 
1  ine. 

Fig.  3  Logarithm  (base  10)  of  the  normalized  scattered  intensity  pre¬ 
dicted  by  Mie  theory  (solid  curve)  for  x  =  100.  The  dashed  curve  is  the 
physical -optics  approximation  [6]  of  the  coarse  structure. 

Fig.  4  Like  Figure  3  but  with  a  linear  scale  and  x  -  10  000. 

Fig.  5  Photographs  of  far-field  scattering  for  (a)  the  bubble  radius 

'{  ~  400  i,n  (x  '  6360,  •  1.3°)  and  j  =  2 ;  (b)  ,i  330  pm  (;,  ■  43/0,  2,2 

and  j  -  1 .  The  critical  scattering  angle  ■!■  is  near  the  left  edge  of  the 

photograph;  <\>  decreases  from  left  to  right.  Coarse  structure  is  manifest 
as  broad  vertical  bands.  The  rings  are  artifacts. 

Fig.  6  Normalized  Mie  backscatteri ng  from  a  spherical  air  bubble  in  water 
(m  -  3/4)  plotted  as  a  function  of  the  size  parameter  x;  for  green  light, 
a  ••  200  pin. 

Fig.  1  Normalized  Mie  neac-backscottering  (a)  and  near-forward  scattering 
(b)  for  m  =  3/4  and  x  =  3040.  For  clarity,  the  computed  I'  '  was  multiplied 
by  two  in  (a)  before  plotting.  The  inserts  illustrate  some  of  the  significant 
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The  intensity  of  light  scattered  liv  an  air  bubble  n  water  m  predicted  by  the  geometric-optics 
calculation  >f  I)avi»  <l!lr>5>  to  have  a  divergent  .mn.br  derivative  as  the  i  ritual  suittenng  unglu 
<J>,  is  approached  FH' a  Is  of  dilVrartion  i  t  the  angular  region  m  ar  •!>  are  (tenant  ed  here.  The 
Fraunhofer  diffraction  for  scatt*  ring  angl-s  .f. ■  ■'>  estimateil  using  a  id’ii, <lifiei|  physienlstpticii 
approximation  A  ringing  and  Oec.iv  of  mi"  i  i r  field  vitrnsitx  i  i  nrediete.t  itiat  ;a  foiomUv  :  inula,-  to 
the  neurfkdd  diffraction  at  a  straight  ‘  lire  Ole-  v  Hum  rf  I'allimeter  radius  h.ihhlc.  ■  it  water  with 
collimated  mor.ocht amin’ic  liluniin.ition  eotifin.t  tte-  existence  of  this  ringing  wl  .  h  has  a  qua  li  froth.*) 

25  mrad.  Ttie  (hflVaction  rulrul.ituni  ;'ar,  an  a|i|iroxi:nat<‘  description  of  tin  relative  of  the 
olmerved  rnuximn  and  .nie.nna  Fringes  with  a  lower  com  mat  and  spacing  •  0  I  mrnd  were  also 
observed;  they  appear  to  he  caused  by  :!i  n.u  rferetn  e  of  ravs  witli  distinct  paths.  Implications  for 
the  critical  angle  scattering  of  white  light  are  discussed. 


INTRODUCTION 

The  purpose  of  this  paper  is  to  present  mi  ,i  ■  i i «  n».  ids 
and  an  approximate  theorelical  iIcm  ripi  an  •  .1  dillr  tnm 
phenomena  observed  when  lit;  hi  i  •  si  ,o  lei  i  d  m  ar  t  lie  cnl  ii  ,il 
angle  Irom  a  Inililtlewilh  an  in'erior  mlrai  five  index  h  than 
that  of  t  lie  outer  media.  The  opt  r  a  I  Wax  el.  1'i'th  is  assumed 
to  he  very  niueh  less  than  t  hi  i  ad  it  is  nl  the  m  mi  Iv  pin  m  .d 
bubble.  Approximation:,  given  here  i  oin  erjiutg  tin  oehavior 


ol  the  reflect  imi  cm  llii  lent  and  the  extent  ■  d'  >.  r.iT.'iit::  wave 
front  \  ii  Id  an  i  xpri  .sum  for  the  seal  I  via  m;  in  I  ae  tar  livid  that 
m  sin i ilar  in  lorut  to  the  Fresnel  invar  fn-ldl  dill'aelton  by  a 
t  might  edge  Oli  .arx  al  e  his  ri  ported  here  ol  l  lie  •  in  i  ‘  i ti ir 
In.  a  i  in  I  (Id  i  m  Water  siip|  ■<  n  I  SI  line  aspects  ol  l  'ns  Ml  null  I  led 
llieotv.  The  I  real  men'  mav  be  readilv  extended  to  include 
scattering  trout  other  curved  Mirl.ii  c  such  as  t Ii  , i  id  a  rylin 
di  r.  prnv.deil  that  the  oi  on  p!e  radii  o|  T  an  vaMire  are  much 
;;  i  eat  i  I  '  hail  t  lie  waxi  del  lilt  h  ol  I  Ig  lit . 

i  1 '.1711  Optical  Sm  iet,  nl  A . run  I '.'05 
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MG.  1  To!  »l  nflorhon  in  H>«* 
plant*  of  s«  .attoiUkj  nl  .1  r . iy  with  .m 
impact  p.u.mictiM  h  l.iujor  ih.ut 
th«  cntu.ai  iinp.M  t  p.uanmtf'r  t\ 
Only  onr  ijun-li.int  of  Hu*  i.foss 
suction  nl  tho  tnihhlo  i:-  shown 


Ha  vis1  calculated  the  angular  disl  i  i  In  it  inn  nl  the  intensity 
of  I  if! lit  sea  tie  red  hv  an  air  Imiilile  in  \\  itei  I  runi  a  i  nil  minted 
incident  lieatn.  rite  ruleulat  inn  w  as  ha  e.  I  nl!  ;;i  mini  r n  <  ■  ( *<  tes 
and  the  Fresnel  relleet  inn  me  I  lien  nt nl  In  In  at  (-lane  miei 
filet’s;  ealeulat  lulls  nl  I  Ills  ( \  [ *«•  eon  nnili  an  .ipprn .  min'  mu 
that  is  uselnl  uiilv  when  the  radius  nl  the  bubble  a  i  very 
tIUH'h  greater  than  the  \\ aveli  n;;th  III  1  lie  outer  ineil.a  Let 
t  h(‘  seat  t eri I e.;  ally le  i/i  nl  a  I  a x  denote  I E ie  dex  in t  mi ,  I II;: le  nl 
the  seat  t  ered  I,  i '  I  ruin  the  dire.  I  mo  nl  I  is  ,  :i.  t  >  i  •  nl  pat  allei 
beam  nl  1  ays  1  si  1  I  i>:  1  ).  A  si,:ml  n  ant  h  a  I  m  t  "I  I  Mix  1  ■  <  a  I 
eulalinli  is  tii  i|  (here  should  l.i  an  ilnnpl  de.  mm  in  the 
sent  I  ered  intensitv  when  the  >.i  at  ter in",  a  n;t  !*■  .  ■  imml, 
where 


Ifi,  IT  -  1 III  .  Ill 

sink,  -■  n  ('.!) 

and  m  =  (n„/n,  I  >  I,  n„  and  n,  heini;,  n  . i mt  1 1 \  >  I ■.  the  re 
frnctive  indues  nl  the  niiler  and  the  min  i  media. 

Figure  I  illustrates  the  physical  origin  nl  this  predicted 
decrease.  Consider  how  f  he  .unde  nl  im  idem  e  at  the  htihb'i  \ 
surface  0  varies  a-  the  impai  I  parameter  nl  the  rax  >  on  lease  ,. 
Hays  with  an  impact  paraniefl'i  in  excess  nl  a  entn  al  value  h, 
have  (I  >  ft,  where  It, ,  t;iven  liv  I'  U f.1).  is  t  lie  critical  aiipje  Ini 
a  plcim'  interlace.  Aeeardini:  to  the  nenmetrm  optics  ap 
proximal  inn,  raxs  with  /■  •  l>,  ii  c.  M  will  be  totally 

relectetl,  Hays  with  h  s  b,  will  In  partialis  tellcitid  ae 
en, din;'  to  the  ina^nilude  ol  the  I  resm  l  relleet mn  eoetlinents 
and  I  his  partial  relleet  ion  leads  to  tin  abropi  di-erea  -c  in  i  Im 
scattered  intensity  lor  I  >,r.  i ealeulati  d  al ti  ruif 

intensity  also  me  I  tides  I  lie  multiple  I  el  Ii  .  I  inn  ami  1 1  iiimiii  • 
stun  ol  the  rays  with  impact  p.ii.iiiu  b  i-, /•  ■  / 1,  Tin  ravs  that 
penel rate  t lie  sphere  an  I  ate  part  i.ills  reflected  enni  ■  ilmte  to 
seal  ti  ii|;  tor  tin  entire  ranee  o|  however,  I  lies  do  not  oh 
score  the  ilist mi  t  decrease  in  the  sent  let  "d  mteie.it  vlnr * 
•A,  This  is  evident  'mm  i  nit  in;;  pi  rdn  i  mn  ( J-T  •  ISnlKel 
II  for  the  in  tens  1 1  s  nl  an  an  bubble  ins  e  .  i  I  b  ii  be.  ~ 
H'TH° ,  it  i,  also  es  idi  nl  in  Lie  mie  m  n  ■  1  lie.;  pi  e.  in  I  nne  I  nr 
bubbles  in  other  Inpiids  I )•  i ■  • .  .  n|  Kei  'll  The  imcnspv 
coat  t  ibul  ion  ol  I  he  mull  ipl\  i  e Ur.  led  and  rel  i  .a  1 1  d  rax  smII 
be  omitted  in  the  treatment  presented  lu  re  sun  e  «"  are  <  mi 
sideritlfj  only  the  elleets  ol  dillr.u  limi  on  the  seallered  in 
tensity  for  small  |.,i  -  | 

The  calculations  ot  I  'avis  p  red  it  l  t  hat  .//.  !<i .  is  de  emit  in 
nous  at  i.'i,  where /,  1 1, ’>1  is  I  lie  tot  ,i  .ilu-ie.  I  rati  nsilv.  I  M 
tract  ion  ot  bj'.lit  will  um  peru.il  ill ,  i  n  be  d 1  ,i  mil  min  ms. 
The  method  presented  here  IS  a  aiaph  nil  )  le  d  Im  e  .l  i.n.il  iiij; 
the  elici  ts  of  ddlraeliuii  ii  the  nation  w  h,  n  u  .  |  i .  .ui.ill 

I'  irst,  till1  approx i mate  shape  ol  l  lie  \  iiloal  w  ,ue  I  n>nt  is  <  al 
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ciliated  peniiu  l  rn  ally  Im  those  ravs  with  b  ">  b,  Hays  wit  It 

b  <  h,  will  la . nl  led  i  ee  Si  <  III.  The  virt  ual  wave  I  nail 

is  I  ben  tnndil  mil  by  a  pli.ee  i  m  reel  inn  due  to  t  he  phase  slid  i 
ot  the  uimmibil.il  !  n  m  l  .m.plil  ude  rellection  i  ocllii  lenls 
Im  h  •  It,  Tim  inn. blmd  wa\i-  Irmit  describes  the  apparenl 
mirei  viewisl  In  hi  •  ■  1 1  cut  lai  Irni'.i  (lie  hnhhle.  (lie 
I'- r . n 1 1 1 1 1 n 1 1  i  (iilli.ni  ,ii  nl  tin.  wave  fixes  the  desired  ap 
prnxiiii.il  inn  Im  iln  ..  ittciim;  atiiptiliide  and  intensity 

The  xv, e  e  I mni  i  .  ■  .  lined  In  extend  over  a  hall  plane  and 
is  not  limited  b\  t  lie  linite  Si/e  ot  t  lie  I  nibble.  Tim  !  real  unlit 
resembles  Aux  -  Ullr.ii  limi  correct  inn  to  I  iesearti  - '  ,;eo 
inetrii  optus  t  re.it  inont  ot  the  rainbow  caustic.1  '  Keen 
thnupji  Airs '  •  .ippinxim.il  mn  t  niisiders  the  virtual  wave  trout 
to  hi  n!  ml  .mti  e\ti  nl  i  .mil  a  di  up  iimsl  be  finite),  it  is  iisel  ill 
Im  the  dniinii  mt  pnlari.Ml  inn  w  In  n  I  lie  radius  nl  t be  sphere 
e  \  •  led  |  he  w  .m  h  1 1  ;t  U  iiv  i  la.  I  nr  nl  approx  i  tua  I  el  x  'Mill.'1  ' 
III  nd  a  i  pa  l  ei  I  i '  u  .i  mu  n  si  1 1.  .pent  criterion  max'  delineate 
t  lie  ranee  nt  n  .a., In  .  .In.,  approx  i  mat  intis  presented  here 
lur  dill i aei mu  U\  i  !,■  nl.le  in  ilm  repimi  nt  critical  aiipjc 

M  e.i  ni  l  ine i  '  .;  .  .  , '  I  .a  i  ii;;  i  it  n/iite  huh  I  bv  an  air 

In il iltle  in  'X al ."  u  i  •  c e. I  1 1\  1  i I,,  -i  i  (  Fi;;  I  nl  Ucl  7 1  I  lu  x 

-ub  i  a  ut  mt  .•  i 1  pi ...  (n  !  :■  .  •  .1  ,  .■■  * ,  . i  ri e  opt  ies  nl  a  rapid  de 

.  re, I  ...  in  /.  :  .|  He  ise  .  .1  v.  I , Oe  lipht  ami  I  lie  lllilll '  .I 

a 1 1 1 ; ) i L 1 1  1 1  ■  .i' n  mn  .  I  i  hi  ei.  .  -in i  meiit  ■  make  I  b.em  ill! -nila hie 
I n r  1  n 1 1 1 1  .ei  .1,  x  ,  1  i.  •!  "  m.'I  mn  |  In  nrv.  The  nl  iserval  mu 

xx  it  Ii  la -a  a  i  II .  i  ie  . . mi  h  u  i  ed  here  in  Sc.  .  I V  md  ii  ale  1 1  ml 

the  aiuplt  bid. i  I .  dii  nix  fixes  a  lair  description  oi  the 

telai  ive  pi.  it  i  .I  :  1  ■  m  r.  nna  and  minima  ol  t  lie  inteie.it  v 
in  I  he  ie  .inn  .  •!  ,  1 1: ,,  ,,l  an; be  si  atteritl|; 


I.  CURVATURE.  OF  1  HE  VIRTUAL  WAVE  FRONT 

Fit" in  Tim  '  i,ii.  .  ,i  ( bn  1 1  s|,m  iTinrilmale  system  iisi-d  ill 
I  lie  de-.  i  iptmu  nl  i  he  x  u  in.il  w  ix  e  I  rmil  I  he  oruy.n  is  ehir-en 
to  in  I  lie  p,  uni  xs  lii  i  e  I  lie  i  ax  xx  it  Ii  /■  b,  ei  ill  tacts  I  lie  sphere 
and  t  be  x  axis  I  •  p.n  dli  i  <  o  I  lie  inch  lent  rax  I ’units  (  v  I  nil 

tlm  xirtn.il  .'..ix.  I.  .nl  an  obtained  by  extendi  ip  t In  pH 

met  r  n  ..'.lx  del  ,a  ni  "d  .  Otei  i  d  r.ax  baekxv.ard  a  dial  am  ■  b 

trnlll  I  In  pm.il  .  '  I ; ,  I  I,  I  .  .  a  I.  1  I  e  t  h|s  cm  inlitl.ll  e  SV  t ,  111 . 

tin  sot  I.  ii  e  ..I  1 1..  I  ml  mil  "i  i  lie  plain  al'  si  at  term)!  is  de-.i  i  died 
by  t be  | mm i  T  i  .x  in  i e  It  m  ■  />,' ) 1  -  (ii  b  l !  ,  x 

It  b,  .  and  i '  i i  athus  .  .M  lie  Imhlile.  <  ’.ah  iil.it  mti  nt  t  In 
slope  m  of  the  si  altered  m  ,pves  the  I ollnxviui'  p.ii.niuTru 
eipi.it Ion  lor  lie  virtual  ".axe  front: 

t  -  ti 1 1  lit  ni  h  1  "].  CO 

x  •  -  "ib ( I  t  in  ')  !  Ml 


I  W  ♦  4 

?  *<•  \  ; 

i; 

i' 

L  :r 


t'r 


I  / 

•  i  / 

!'  ** 


f-K'i  ?  IIUj  ‘Ttimn  ui  liio  M  fo*  lo<:«Tlin<)  ttw  ci*wdin.iWr; 

\*.V)  of  u  potol  Onr  Mi.*  «n v«  I  v ■  r l •  t  t'  v,  »v*’  Ituftl  Tfio  uppi'f  hflr  of 
/•  is  par.iH*'!  ttt  Mi. ■  »  , j v . /i ,  ,  .1 to.*  rl.-.i.fnt  »*  it,->  ptTint  of  contai.t 

^tlfi  Iho  fmi'f’M*  ■  i  I  f  *•  | I.M  f  t  I  [ft  •,»[  ititl  Mv*  point  L)  The  rotat.-  1 
f.iK>f  i  » *  *•*'.  t  x  t’f  Hit:  !■.  »‘vt  hi  f.  ii  ilit.iU*  .1  laylui  »*»p.»nsi(  *o  ollhn 

W.iVU  Iff/''* 


1 


do 


Philip  1,  Marston 


m  =  titm/>  -  tan['.!  ms  1  (/»/«  )|.  (,r>) 

The  calculal  inn  of  I  he  dill  on  I  nm  i .  Lit  dilated  hv  tic.,  nhmg 
virtual  wave  trouts  using  coordinates  lu.r)  which  have  the 
same  origin  hut  are  rulatril  •  lockvv.a-  hv  an  ancle  t  ?2) 
so  that  the  posit  is  e  e  axis  is  l  lie  pal  h  nl  the  seal  I  ereil  ent  leal 
ray  predieleil  hv  geometric  opli.  .  Equation  (I)  and  the 
rotation  traiislormatioti  give  tin  Inflowing  expression  tor  the 
new  coordinates  of  the  virtual  wave  front 

u  -  x  sin '.’A,  t  y  m  'H>, ,  (*‘> ) 

v  -  —x  eos ‘ill,  4  x’  sinVtt,  (7) 

i’ll  obtain  a  I'avlor  expansion  of  etn)  ahoiit  u  (t,  the  Id  I 
lowing  derivatives  were  computed  for  0  -  H,  dr/tln  -  / *  "i,- 
linti  (I  -V/i/u  -  -  (/  i/p'i/H  t  /’ll  '</ 1;  t.'i/tt  whi  le  pill)  ill /till 
and  (,(9)  =  dii/dH  The  caletilat ion  of  ;>  and  >/  used  the  re!  i 
tioli  i,  “  »i  sind  winch  is  evident  from  Fig  I  I  lie  let  in  ot  t  he 
Tliylor  series  which  Is  linear  in  u  vanishes  liciaii.-.e  ol  the 
coordinate  rotation,  this  is  evident  since  p1'1,  i  n  anil  i/h',  i 
--  a  eos  II,  ( )m  is-  loll  ol  collie  and  III;:  I  nr  older  In  Ills  in  the 
Taylor  series  gives 

v  ■■  mi2,  (Hi 

nt  ■  —  l/o  corf',..  (VI) 

For  the  case  of  an  air  luihhlc  in  water,  nnmerical  calculations 
show  that  K(|.  (H)  is  within  10'  i  ol  tile  value  of  t  (u  >  predict cd 
by  F,qs.  (I )  (7)  provided  that  0  •  ii  '  titi.'ii  In  the  <  ah  ula 
tionof  the  diffraction  that  lollows,  Kq  ihi  will  be  used  tor  the 
region  u  •:  0,  (Airy's  treatment  ol  the  rainbow  also  retains 
only  the  leading  term  in  approximating  the  virtual  wave  trout 
but  in  that  ease  it  is  ii  cubic  term.)  I'he  sign  "in  indicates  that 
the  reflected  wave  is  divergent 

II.  AMPLITUDE  OF  THE  VIRTUAL  WAVE  FRONT 

The  physical  optics  approximation1  involves:  m  tin 
specification  of the  polari/nl  mil  ol  the  mm  lent  wave,  and,  i  ii) 
the  use  of  the  reflection  coefficient**  h|r  a  phuir  surface  to 
approximate  the  refleeled  lields  at  each  point  on  a  curved 
surface  In  the  application  ol  the  teehnnpic  that  follow.,  the 
incident  wave  has  an  elect  rn  field  polar  i/a  I  mu  cithci  cut  irclv 
perpendicular  to  the  plane  of  .ittcriiii' or  entirely  parallel 
to  it;  a  subscript  )  will  be  as  .igned  t  he  value  I  iii  t  lie  fn  a  e. isl¬ 
and  2  in  the  second  t  A  linear  com  Inna  t  ion  ol  reflected  In- ids 
may  be  used  lor  other  incident  polari/.il  mils  )  I’o  a  di  a  ant 
observer,  the  ret  lei' ted  wave  appears  l  o  c  ■nm  Ii ,  ,m  virtual 
wove  describeil  hv  the  follovving  i  oiaplex  ampin  ode  at  points 
on  the  it  axis: 

I’di.f)  --  r,,“'  >■■■"  “>',  till) 

where:  t  =  \  —  I  ,  in  is  tin*  angular  Irequeiu  v,  /  is  the  tone.  /■ 
=  2rr/x,  \  is  the  wav elengl  Ii  in  I  lie  m in  i  medium,  and  r .  is  i  lie 
amplitude  reflection  eoeltu  leal  ol  a  plane  siii'aec  tor  tin  angle 
of  incidence  I1'': 

cos  il  --in  -  sin  ft  t 1 

r,  -  ,  1)1) 

cos  II  -i  In  -  on  ID ' 

-II  '  -  Cos  II  4  I II  -  -III 

ri  =  -  . .  ■■  ■  -  !  1 ,  1 

n  '  cost)  4  ( u  '  —  sin  'll  1 1 

Reflection  lor  h  '  b.  where  b.  -  /i  •  >  n  i  an  In  lies,  ribed  hv  t  lie 
following  approx  i  mat  ion  -  that  ::iv  e  1 1  •  ■  leading  dcpcu-ti  n.  i 

1207  .1  Opt  Sue.  Am  .  Vol  t'.'i.  *i  Si  pi*  nilnr  |;i/'.i 


of  r  i  and  r  ,  on  i  ~  II,  ~  II  >  0: 

Il  s  i  -  (8</n  cost/,  ) ' (Rfl 

r->  ~  —  1  +  n2(8</rr  eos 0,  ) 1  /2,  (14) 

and  in  each  case  \ilrhlll |  ■■  <  1/2  When  0  >  I),  ,  the  reflection 
eoeltu  units  are  unimodular  and  complex;  they  are  r,  - 
expt  -i(i|l  and  r  <  -  —  exp(  —  iiV.)  wher"'-': 

tan(A|/2)  =  (sin-ff  -  n  ’ ) 1  /2/t* osf),  (IS) 

tunt  A.,/2)  =  n  2  taut  A  |  /2).  (18) 

Reflect r  hi  for  b  V*  h,  where  h  -  b,  «.  n  can  he  described  by  the 
Ini'  . wing  approximations  that  describe  the  leading  depen 
deuces  of^i  and  «/•  on  >'  -  H  —  0,  >  tt: 

<\|  v.  (Hi'/ri  (osf),  O'2,  1 17) 

h>»  ~  ’j2(8< 7/1  cost),  )' (|8i 

Approximations  for  r  in  terms  of  a  may  be  obtained  using 
ii  n  u/ipii,  )  For  u  -f  (i,  Kqs.  (17)  and  (18)  give  ht  ~~ 
i'jIi1  '  where 

Ii,  (8/  rm ) 1  '2/(  osf), ,  (19) 

(i-i  ~  n2li  .  1 20) 

The  spec  it  lent  uni  ot  '.'{it  I  >  is  I  urt  her  simplified  hv  t  a  kite;  r  \ 
r  ■  =  o  I  or  H  0, ,  and  equivalent  I  v  for  u  <0;  the  use  of  this 

appi  ix i in.i I  nm  is  nml  ivaled  by  t  he  divergences  |  dr /dll  \  and 
•  .’/r as  «  -  I  tl  and  <;■  <t\  •  H),  respectively 


III.  FRAUNHOFER  DIFFRACTION  PATTERN 

T  he  scattering  amplitude  at  large  distances  from  the  bubble 
is  proportional  to  the  Fourier  transform  of  r(  exp(-/7.*<*u  T: 

/,(())-•  f  rl<u)vd»-"“‘-"‘idu,  (21) 

i'  •  *■ 

w  here  1/  -  <;■  (an  angle  of  the  observation  point  relative 

to  )  is  a '.sumed  to  be  small  so  tha,  .sin r;  ~  t /;  and  constant 
rmiginttidi  mult iplnative  factors  tii.it  precede  the  integral 
have  been  omitted.  The  approximation  described  in  flic 
preceding  ..c<  turns  yield  the  following  exp  essioti  for  /, 

/_.  -  i  (  c  •(*<•■■* ♦  90 1 ♦  i*,u '  bdit,  (22) 

»/u 

vvticic,  here  and  m  the  equations  that  follow,  *.he  upper  sign 
is  lor  poluri/ui  mu  i  I  and  I  he  lower  one  lor  )  —  2. 

Kvalualioii  of  Kq  is  faeilitalei*  hv  a  change  of  varialile- 
*  ha  t  (nmpleti's  the  partial  square  formed  by  (<iu2  )  i/u  ) 

I  )i Tine  .  i  lo  M 1  ii  -  ii  where 

tc  -  i) |(o /A  i  ( osf1,  | 1  i '.>:() 

We  will  approxitnale  the  integral  otilv  for  0  <  •/.  <  whit  It 
corresponds  to  the  region  predicted  hv  geometric  optics  to  he 
strongly  illuminated ;  this  rest  riel  ion  gives  te  >0.  Fa  pint  ion 
(22)  may  I'len  lie  written  as  lollows:  f,  -  4  ( ,\o  cost),  1 1 '  'g/2, 
where 

j-  ,al  4,1"  a  1 1  (t‘/  |,.|  ^  ,.-iv'io.iil(,il</'.,li1(/;  (24) 

Vf  l.'.K  •  *  d/.’  ’"|tie  *  .-I1  '  -  ii  1  '")  ( X  (i  eosd,  )'  1  1251 

Philip  1.  Marslon  1  Jo') 


In  the  exponent nil  tli.it  precedes  t lit*  integral,  ^(ie,(>)  is 
evaluated  hv  the  a|i|iro|irmle  variable  substitution  In  F.ip  t'JM, 
(0  replaces  ir  and  ir  rcplat  es  I;  phvsicallx ,  v1 ,  be.O)  is  the 
approximate  phase  shill  associated  with  r,  til  the  "eometric 
ray  scattered  to  the  obsei  vat  ion  ancle  o  The  second  temi  in 
this  exponent  nil  arises  I  mm  I  be  coinplet  ion  ol  the  sipiare. 

We  now  examine  the  condition  tor  oinitlinn  the 
exp|tV  .a  )|  lai  lor  w  it  Inn  t  lie  inle;.i  al;  when  t  lie,  sinipldi 
cation  is  permissible,  the  inlcpral  mav  be  estimated  usmt; 
Fresnel  integrals.  II  ^  ,  I  .\n  I  is  a  soil  n  lent  I  v  slow  I  v  xarviii;; 
function  of the  principle  i  onlrihui  mu  of  the  inteprnnd  oc¬ 
curs  when  | .’  |  7  beeaiisi*  the  phase  ol  e\p|i  ( ,t  /'_!  I  c"  |  is  sta¬ 
tionary  at  =••(!.  t  For  a  hn  id  discussion  ol  stilt  inlinrv  phase 
aiiproxini.it  ions,  see  Seeel'l  Korn  •  0,  y  ic ,  and  j  >  |c|. 
it  is  evident  that  | y'  /  I  :  n  l|  •:  ^  ,!:•  .01.  >  nnsetpientlv  the  factor 
in  (pies t inn  mav  lie  omitted  it  the  radius  ol  the  bobble  is  sot 
fieiently  larne  t1  it  i/ .r.!,il)  «  I,  For  e-,,imple-  water  has  o 
I  ititll  and  visible  li|;li|  from  a  lie  W  laser  has  a  wavelength 
in  water  of  n  1 1> dtitT-t  pm.  For  pol.u  i/at  ion  1  the.  eomln  ion 
nit  i^,(-.bl  tpxes  pi/.a',l  pmd  1  >s  I  and  tor  pnlari'/al  mn  d  it  stives 
(n/filio  pm)1  1  ■>  I  where  I'.||S.  t’Ji,  (  I'll,  C.’tl),  and  (dal  have  also 

liven  Used 

Alter  some  rear r .t tip ■  ru; ,  the  result  ol  this  approximation 
max  he  written  a-  tallows 


K  ■■  «"l  )  -  /-'I  ■  •  >|,  c.;t;i 

where  /•'(«)  is  Fresnel's  inieeral1"- 


and  f’l--")  =  St  ',)  s-  >/,.  For  any  state  of 'incident  pn 
lari/al  inn.  the  billow  mi;  fact  nr  pixes  the  dependence  of  the 
intensity  on  i/  where  n  *  is  the  complex  'onpipate  ol  r: 

RR*  =  ((’(11  I  -I  I •'  I  |S(l.  I  t  Va| CIO) 

The  horimili/at  ion  used  in  F.i|s.  t I  )  and  l. lib  is  such  that  rh* 

*  ‘2  as  b/i,  -  i/d  heeotnes  lar-’e  and  positive  As  expected,  hr* 

*  0  as  </>  ••  i/i,  becomes  lari'.e  ami  positive,  however.  I  lie  rej'icm 
te  <  II  is  out  o|  l he  ranee  of  lie  validity  I esl  1 1 >r  'he  approx i 
mathm  used  in  denv  i n p  !•'< |  i.lio 

Davis'  a  nil  I  v  sis'  ol  I  he  jicoinet  rie  a  at  lei  nip  mav  be  used  to 
obtain  the  m , i  mall/, il  ion  lecpiip  d  to  expre  the  absolute 
intensity  /  associated  with  the  totlcctcd  rav  at  a  distance  /» 
»  a  from  the  center  ol  the  hobble: 


/  /„l(l/AT  /"’VS,  (111  ) 

where  A>  is  lie  incident  intensitv  and  rr*/H  *  '/,  for  lurpe  n 
so  tl'.at  ha  |  i.l  I )  c  m  in  ides  vvu  li  l-.c  | .  I  a  I  ol  I!  el  1  in  the  appro 
priate  limit  Furl  liermnre.  il  the  Intel  t  c  lencc  helween  tin* 
totally  reflivled  lavs  and  the  multiply  relict  led  ravs  are 
olnil'eil,  the  total  intensitv  /,  bceomes 

/,  /,,<(»  * /c  )  'l.m-'/M  I  (,'  i:n>) 

where  li'ICit  is  the  total  rot  it  n  I  ailion  hv  the  limit  iplv  reflected 
rays.  Fur  an  ar  Imlililc  m  water,  lilunnn  ited  bv  impel. irt.ed 
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t-ICl  I  View  ol  1hi*  appar atus  in  the  pi. me  ut  '.c.at’.eriiH)  which  was  hoo- 
/ontal  m  llm  lah  -nib  ry  Ihn  tmMom  ol  the  bubble  was  in  ,  oillai.l  with  tl 
n.ieille  lh.il  pi  everivn  llm  h'ihh!i  loini  llu.c.ini)  out  ot  the  '.culieono  plane 
The  top  ot  tlie  bubble  (wlm  i.  i  .  the  snte  He, played  here)  was  viewed  by  a 
microscope!  The  wavclenqth  m  an  ol  tlie  illumination  was  0J2  8  nrn 


lij’ht .  Table  1 1 1  ol  Ifei  1  shows  that  (>'  is  a  smoot  b  In  net  ion 
ol  </i  lor  <:>  ne.-ii  ■  a i id  i  h. 1 1  t 1 1,'  ~  H'.’.it")  =  O.UO'.’T  and  fi 
(7.ri°)  ~  ft  i »i t  ('  (  onseipiently  (/  lor  an  air  bubble  in  water 
is  approximately  p,  in  "I  if  t  lie  asymptotic  value  til  /,'!,' V  M 

Kcpiation  1  ab  i  .iihil.il  m  |.*i  hi  lo  I  he  well-known1"  Fresnel 
( lien  r  fir  hi)  dil  1 1  ,e  i  no  bv  a  si  Me;  hi  eilpe  where  the  para  in 
elei  u  is,  in  tli.ii  .  a-.e,  proportional  to  the  lateral  displaeemenl 
of  the  ol.-  tax  al  ion  pmiu  In  in  the  cdpc  of  the  pculilctrii 
shadow  For  hot  h  i  an  m  I  lie  ref,  ion  u  >0tlv,  intensity  is 
predict ed  to  o  a  dl.de  wil  l\  an  osi  illalion  ampl.itude  that  tie 
creases  rapid! ,  with  nn  re.e  mp  re.  This  behavior  is  predicted 
for  the  / m  iirhi  ;  1 1  <  iisp  v  of  a  Imbble  because  of  the  abrupt 
variation  m  the  m  plitude  of  the  virtual  wave  front  near  the 
critical  a  ti|'  le  and  Ins  a  use  of  t  he  '  iirriil  HIT  of  this  waive  I  rout 
inti  need  by  tie  retie,  turn  I  rotn  the  curved  surface. 

From  the  alorrtoenl  to.  ed  bound  on  |*|  for  the  principle 
eont  it  hot  ion-,  oi  the  mtei  i  al  in  Fa|  (’.fit,  we  can  est  in  ml  e  bn 
th"  ic  rep  ton  of  i  til,  rest .  the  errors  introduced  in  t  lie  ealciila 
tii. n  nl  the  'i(  Iron,  Fn  .  ( lb)  and  (.’HI.  In  the  interval  (I  ■  n 
<  it.  K.(|  (ittb  pi  ed  ic  I  ■  I  xvi  i  intensity  maxima.  In  this  reiimn 
the  pi  rieiple  i  out  i  ibip  b  ms  occur  with  (I  •'  it  ‘.i  ■  i-  t 

ic)t-  A '  lot1  'with  in  -  it.  With  the  atoiemeiitioni  :l  A 
and  with  n  I  n  il.  the,  estimate  pives  ~  Ml  pin.  A 

coin  pa  i  ison  ol  F.,|  i  1  i  xx  il  b  A  i  ealeillated  llsilip  K<|.  I  lb  I  and 
tills  a  sip  >vvs  that  I  In  ei  i  ■  <r  ini  mill  a  ed  ill  r  > )  is  less  than  TV  A 
similar  con ip.u  i  <  n  n|  I*  > |  -  '  I  Ml  and  (JOt  pivc  a  hound  lot  t  lie 
error  in  ■  "I  s'  <  'oic-cipn  ntlv  the  use  nl  t»,  if, a  1  is  jus 
t  died  |-  i.i  I  In  t  m.  a  e  him  i;  ,„i(,  n  ,  I  lie  upper  I  i  ill  1 1  nl  ill 
teKrntimi  m  I  u  i  n  In  a ild  cixe  a  iim  IiiI  ipproxiiiial mu.  The 
main  jo  1 1 ! n  ,i»i  ii  l ",  i  la  ippi ox imal o .-is  desi  rilied  in  bee.. 
II  III  n.  tin  m am '.e  i  am  xv 1 1 1 1  experiments  described  in  the 
sect  Ion  (li.it  ; ,  i  I.  i'.v 

IV.  OBSEIRVtD  SCATTFRING 

I  ills  si  i  l  inn  .  I  -  ,  1 1  le  pilot  i  i  aplue  obserx  al  lolls  ol  the 

er  1 1  teal  a  1 1 ;  •  I  <  v  at  ia  t  i.ais  m  the  -a  at  1 1  red  ill  tens  it  v  pi  edit  ted 
in  I  In-  pi  a  vi,  ue  ■  i  ,  i  inn.  I  more  it  -In  uv  -  a  diagram  ol  I  lie  a  p 
)iu  at  os  I  )i  l  died  wan  r  was  placed  in  an  alt  i  mi  nn  in  eon 
tamer  The  mm  .  wall  .  ol  llm  emit. liner  Ii.itl  been  blaeki  tied 
so  that  tin  x  ab  m  lie,  I  mm  h  ol  the  si  rav  lip  It  l  A  steel  needle 
passed  via  a  u  all',  i  In  ou"  h  a  rubber  septum  ill  the  lint  loin  ol 
I  he  eon  I  ii  ue  i  I  lu  o!  In  r  end  o|  I  lie  needle  xvas  eon  lb  ,  led  I  o 
a  sv  *  i 1 1 ; ; e  tli.il  ,  oiii .mu  d  air  I  lor i/ ontal  boles  I  d  cm  in  ill 
a  i  net  ir  pa-  i  I  Ii  i  em.b  tlm  wall .  ol  I  be  eont  alder  a  I  riplil  all 
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FIG  4  Photograph  dIUmi  critic. il  angle  tar  )..*i>!  tliOt.c.  tux'  inim  ,.  "*  m-  .  • 
hands  have  a  charartenshc  separation  ol  nppioonialnlv  . i< *  'h.  v  n.  m 

point  1/0(1  with  its  eioctoe  Hold  patallol  t'.i  Iho  si  aiiornnj  p!,i>u 


t m.; .  i ■  >, jo •  .  iii  ciease,  (iom  left  to  right  Tho  vertical 
■  ■a  , *i y  . * ■, i •  i-l  hy  iiiiliat  lion  Tlio  incident  beani  was 


gles.  These  holt's  worn  sea  It'd  with  O  mii:  •  in  |  ■. » ■  m  i  wat  ■ 

It'ii ks  and  wore  co\  cri •< I  by  ant  irt'llct  t  imi  t '  tali  t!  v  nn I  .v  A 
bubble  tiinld  lie  i-reiilcil  Hi  the  tint  tip  "I  'he  in  edit'  b\  ■  !> 
creasing  the  vtihtmt'  ol  the  syringe,  this  i it 1 1 ■  I •  It •  wa  stal'h 
provided  its  volume  was  kept  sulficielitlv  mall  The  biibble 
Wits  approximately  spherical  in  it>  eenii a1  o  ::t"t>  'IA.i  .■ 

discussion  ol' t lie  stability  id  bubbles  -a impended  1 . . .  1  d1'  ■ 

s«*t*  Kel.  II.)  The  base  nl  t  lie  I  nibble  wn  ■  in  n  no  a  a  a  it  !■  I  b 
needle  during  the  ohsers  at  ion 

The  polarized  beam  Iron)  a  h  i»W  lie  Ni  la...  i  w.e  ,  > 
pantled  to  produce  an  K  inni  diam  parallel  heats  I  be  e\ 
ponded  beam  passed  normally  and  hnn/nnlall\  1 1 > 1 1 >i ip b  .no 
of  the  windows  m  stub  a  wav  that  Ibe  central  rats  <>|  i  be  l»  am 
passed  close  to  the  tenter  of  the  bubble.  A  eamei ..  >.  a- wed  lie 
bubble  through  the  secolid  window  I  be  needle  and  ibe  inner 
wall  of  the  aluminum  container  were  blackened  in  ri  dnee  (In 
intensity  of  stray  lip. bt  that  would  have  nth  rw  .  .<■  nl>  mie.l 
•  he  scat  t  cri  np  nl  t  he  bubble.  I  .i  it  lit  seal  lei  etl  bv  t  lie  Iniblli 
near  I  he  crtl  ical  scat  lering  angle  had  an  angle  ul  me  a  ini'  <  ai 
the  wiiter  inwiliiw  interlace  ol  approMtnateK  'm1. 

This  light  is  retracted  at  the  window  air  mini  lace  Mo 
camera  telephoto  lens  (Nikkor  (.f,  manulactured  bv  Nikon 
I  ne, )  had  an  el  feet  ive  local  lent:  I  h  /  i  it  '.'nil  mm  I  le-  Imi  •  u  a 
placed  with  its  axis  horizontal  and  appmxiniaiclx  p  i'  allt  I  I. 
the  refracted  i  ritical  lavs  The  polai  i/alnni  I"  On  idee1 
mount  'brumal  ic  light  was  tisiiallv  chosen  wit  ii  t  lie  t-le.  1 1  a  lield 
parallel  In  tile  scattering  plane  that  n  a  n  spinel  i  ■  ’  'mill 
preceding  analysis  Tliis  choice  lor  the  polari/aiimi  go'e.s  tin 
greatest  magnitude  lor  the  curl  I  o  mit  ol  M  '  1 1  1  ■  t  'll' 

( 'onsetpii'llt  l\  it  gives  the  most  lapld  decrease  tin  |/,|  with 
i  tur  tuts  111),  i  it  to  I  it  should  lit  the  polal'i/.at  toil  Ini  vt  bu  h  I  In 
omission  nl  t  he  wave  In  ml  wit  ho  ’'Din  Ki|  r’'.!  I  will  pi  min'  e 
the  smallest  error  in  the  pirdi<  o-d  scat Icriup  I  I m  diainco  i 
ot  the  bubble  in  tin  ■cut  I  mm:  plui,,  w,i  n  "  'urd  with  a 
microscope  loeated  u,  the  air  direetlv  above  lie  buhnli 

Photographs  were  taken  with  the  camera  h  o  .  Inclin'd  o 


l.no-  'h  .l  i’i.  (iliii  plane  was  at  the  optical  transform 

I  i  ii. i  lie  hie.  l.ri  i)  •  s  .;i  where  <y  ami  i/>,  ’  denote 
ni.  ncs  .tile  i.o.'.  n  it  side  the  i  ell  i  hut  correspond  tor/i  and 

;i  Mu  '.'..is  i  i  i.'i.i  .|"eoi|\  increments  iii  the  distance  / 

•  ir  Ii.  i  1 1  i .  i  ;. i . • , •  h . .  a  ;'..iloe  ate  ft'lnletl  to  the  increments 

ii  \  , \  I'  im  \l  ■  ■  I  Numerical  calculation  ol  the 

.a  so  i  ela  ••  and  M  e  w  iter  air 'nlerfiires  show 

I I  a  •  ■  ■  m  in  a  m  c.  hi  related  to  the  corresponding  Arj 

le.  'n  H. A  ii  wiili  a  negligible  error.  1  he  success  ol  this 

i.„  ill...!  n  . . a  i lie  i  lion  e  ol  the  plane  of  the  viewed  window 

hihai  n  i .  ne. ii Iv  perpendicular  !(•  t lie  critical  ray.  These 
it  .nits  i'i  i\  tie  eoinl'iued  to  give  St]  ~  Al/nb  (Ot  course,  if 
iii.  i  !•  t  i r .  i <  i- 1" i •  i  \ t cl  mil  to  the  sphere  had  been  uniform  so 
■  I  is  i  I..  ■  w  "nl. .  a  v.  s  needl'd,  it  wnuh  I  not  lie  necessary  to  in 
\  -it'".  . '  i  *  i ".  s  i  >  1 1 . 1 1  n  nis  Then  ,\l  'n  h  could  be  used  as  the 

1 1  m  a  .  ..  o  n. Hie  ii  r  iii  Ki | .  r.’li  because  ol  the  Fourier 
'  i ..  1 1  I  •  1 1' : . 1 1 1 1  i  ii  i  >|  H'l  t  le  .  i  it  a  lells. 1  '| 

I.  ...  i,  |  i '  it  im  i  e|  'indiii  i  inn  nl  a  photograph  of  t  be 
niii  al  oicl'  •'  it  let  mu  wit  b  a  bubble  diameter  of  I  fi  -fc  <M).ri 

•  i : i  I  lie  •  ah  iiidu  all  a  At;  ielai  ive  to  an  urliil  rarv  refer- 
...  c  "  I  hi-  him  ii-cd  was  Tri  X  with  an  exposure  time  of 

I  i.o  fl.c  in . i  ul  \  ci  t ii  id  dark  bands  wit h  a  characteristic 
I’M  in:  I  i  pi  i  \  i  mat  cl  \  mrad  I  —  1 .1°  !  are  examples  of 

iln  mi . .  r  ia  1 1 1 1 1 1  i  predii  led  bv  hlip  ('!’.!)  The  concentric 

,  ii.  ul  .  ■  1 1  in;  i"  ib.it  ai  i'  i  entered  close  to  the  most  cent  nil  dark 
tic..'  .  i .  ii  1 1 .  ii  i  ni  l\  .m  a. '  if  act  asm  h  ialed  with  spurious  re- 
t  ii .,  i  • .  ■  1 1  i  ii  Ins  I  rum  i  In-  exit  window  or  from  I  lie  telephoto 
n  u  I  si  ilok  nr  "Hi  mi  the  let  l  ol  t  be  tip.  ure  is  indicative  ul 
tic  n  In.  inni  m  niieiisitv  (bat  occurs  when  ic  <  t). 

Mi  i  n|.  te-iiinio  ti-r  -.cans  were  obtained  from  two  negatives 
i  .  is.  1 1 1 1  ii.  a  ipi.mi italivi'  comparison  with  the  theory. 
I  t,  ,.  c. imp. in  s  i  in  relative  intensity  interred  from  a  scan 
:  i  m  ,-,sl  i  "  lor  a  Imlihli  diameter  nl  I  .lib  b  O.Ob  mill.  Two 
Ins  pai  mu  h  is  we  e  used  lor  each  negative  to  make  the 

■■  1 1  1 1  ■  .  m  w  il  h  t  he  I  hem  \  possible:  til  Tile  local  ion  cf  the 

. .  ■  a  ■  i  .  i .  si  1  iv  could  not  itc  determined  eX|N>rimenUilty 
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Ifomn  microrlonsitonwler  sc.m  ol  it photi  >•<  nmi.ilivc  It >< .  i u il it >li- 
diameter  war.  I  60  r»'*n  llio  solid  I  ww  i  .  mvii  t.y  I  a  t  *u I  which  r.  the 
result  ot  the  phy'.n  ,tl  '"ptu.’i  .t|>prn*tm.,li"n  I  we  i  .ii.irnrUirs  w*iie  adiustiHI 
to  facilitate  tins  comp  msnn 


so  tliilt  an  arltitrai y  constant  aii!:lr  was  added  to  the  abscissa 
of  the  measured  intensity  profile  in  optimize  the  agreement 
with  the  theorv.  In)  The  ahsohtle  intensity  could  not  lie  tie 
termined  experimentally  solhat  th  experimental  miensitv 
/,.  was  obtained  by  tnnll  i|>lvm;\  measurements  liv  a  constant 
factor  K  that  was  selected  to  optimize  ai'.ieeincnl  with  Ivp 
(HO).  'I'he  mat'll  it  ud  c  of  (» '  was  sill  ficienl  Iv  small  t  hat  il  was 
neglected  in  the  compulation  of  ihc  theoretical  prediction  -o 
that  f>f>*  could  lie  used  as  the  prediction.  The  relationship 
between  the  measured  transmittance  ol  the  negative  >  and  I, 
is'-': 

/„  =  Kir  -  r,,1'1)  (Tt) 

where  /o  was  measured  at  an  unexposed  region  ol  the  imitative 
and  the  eontrasl  index  >  was  estimated  Innn  manufacturer's 
literature  to  he  0.70  for  the  conditions  ol  lilm  development. 
The  t(ti)  measurements  were  smoothed  by  hand  to  reduee  the 
mnitnitnde  of  the  circular  frini;e  artifact  in  the  itilcind  /,., 
'I'he  scale  of  the  smoothmj'  only  slightly  exceeded  the  spa.  nit; 
of  those  fringes.  The  validity  ol  I  he  proeeduii'  was  su list  an 
liated  hv  eomparinj’,  parallel  hut  disluu  t  scans  from  the  same 
negative, 

Figure  f>  illustrates  .I  dost'  agreement  ol  i  lie  relative  pou 
t  ions  of  the  maxima  and  minima  ot  I  lie  miens  it  \  hands,  u  it  li 
those  predicted  hv  Kip  him.  Tie  n  values  weie  calculated 
using  Fqs.  I'd)  and  h.’O)  with  n  I  TT:;.;  and  'u  1  » >  mm; 
('(li)  and  .S’t'c)  were  olitamed  In. in  t .  1 1  ih  ■  The  i;.:,  *  values 
and  loeal  ions  lor  the  maxima  and  minuii.i  xh  ir  < >1 1 1 . lined  limn 
Fresnel’s  original  i iihnl.it  ion  lor  i  he  de  a  i  ipl  ion  ol  edge  dil 
fraction  (as  ipioted  hv  N ! a  . i r  l '  1 )  There  is  in  apparent 
discrepancy  ol  dim  I  hem  .  with  tie-  I.  .11.  ay  im  teal  ore.,  ol  the 
measurement;  III  there  appeal's  to  lie  a  iliileicni  e  ut  inten 
sities  ill  the  region  n  ■  I  tli.it  r.mnol  l ■  e\pl. lined  hv  the 
omission  ol  I  ’  I  nun  the  call'll  la  l  u  m;  In*  I  In  amplil  ude  ol  the 
intensity  oscillations  exceed  .  the  prctli.  I  Ions  I  ol  n  '  anil. 
( iiil  I  here  is  an  unexplained  drop  in  !  lie  i utensil  v  near  ic  a 
Allot  her  neg.it  t\  e  shows  1 1 1 1 . 1 1 1 1  at  i\  elv  aiailai  disagreements 
hut  also  a  similar  agree  me  u  I  lor  i  lie  lo.  a  I  ion  o|  the  maxim.i 
and  minima  with  the  th.->«rv.  The  'pi.miiiv  .'ml  quality  ..t 
measurements  t  hat  can  lie  ol  it  an  led  with  1 1 ..  ix  a  dahle  apt  i.i 
rains  are  not  adei  piale  lor  del  cl an  mu  i,:  tin-  •  I  ■  1  u  in  an  . .  a  i  his 
ealeulalion. 


A  comparison  of  photographs  taken  with  huhtiles  of  dif¬ 
fering  radii  shows  that  the  angular  spacing  of  the  intensity 
oscillation  decreases  with  increasing  n,  which  is  ox|M>ctcd  from 
the  lorm  ol  Kip  (  This  decrease  is  also  evident  to  the  eye 
if  the  radius  is  increased  as  the  bubble  is  viewed  thrmp.h  a 
teles*  ope  locus'  d  on  uilmity  Intensity  variations  similar  to 
thoge  displayed  in  Figs.  I  and  h  were  also  observed  with  the 
incident  polarization  of  j  -  1 

V.  FINE  STRUCTURE  IN  THE  SCATTERING 

Cnrelul  exainin.it  ion  of  the  photographic  negatives  revealed, 
in  some  cases,  an  *r. ■  illation  in  the  intensity  periodic  in  rj  which 
was  superimposed  on  the  diffraction  stiurlure  described  in 
See.  IV  1  1  For  example,  oscillations  in  intensity  with  a  period 
ol  ap|!i"'  n.  ol  1 1  ■  mi. id  are  visible  in  the  original  photo 
graphs  J •  i  in!  ■  Im  cd  hi  Fit!.  I.  t  It  is  unlikely  that  this 
st hi.  l ui i  x  Hi-  a  in  t  he  printed  version  appearing  m  t Ins 
I* •< t r t i.*l  ui. i  n  ..  *  i.l'. e.  I  in  Fig. 'j.l  Fxamipnt ion  ot  the 

negiii  .e  I.  ■  iii-.i  1 1..  1 1 •  i !}•«•  spacing  i)  decreases  wit h  in 
ere.r  in,  1 1  i  ■  up  nil.  i ose, ipe  was  used  to  measure  >1 

fr., m  ih  leg  ■  a.  i.'i  hobble  *  with  1 .7  mm  C  2«  $  T2  mtn. 
The  m  ul  u.  i.u  i  « i  inn  I  he  fitted  by  the  equation,  Fj  = 
\i.\  .ip  \  ,  '  I  t  o  t,  however,  it  was  difficult  to  oti 

tuii  *  >,ii  i  !>  .'in.  ui  elm  nl  s  from  each  tiegat  ive  because  ol 
ii.u-i  and  T  •  '  ,  oil, el  ol  the  oscillations.  K.nergv  lisso 
ciai  oil  with  • .  i  ‘  no  i  i\  •  incident  with  h  <  h,  will  eventual  I  v 
lenvi  il-'  '  T'1  .’im.  t .  ( |  with  i/i  —  From  the  linearity  o| 

t hi'  i  qua i  ei''  l'  uel  i  Im  magnitude  ol  A,  Il  is  likely  that  the 
oseiltat  .<••,  m  ,  .m  ,  d  l  x  t  he  interference  of  that  light  wit  li 
the  ci  it  ie, ill ,  o  lie  i.  .I  light;  however,  the  numlier  of  rel’lee 
twin-  .ual  /  x  Pm  •  ,i  i  he  path  of  this  interfering  light  has  n. >1 
lieeii  d.  ■■mi".  I  I!  a  ep  uat  ion  ol  t  lie  niaioi  vir!  Hill  soul  ee: 
foi  ■!■  ..  i  •  i  \ 

Sum  la  i .  i  !  i  ’ .  ,•  ii  ed  oscillations  can  he  observed  in  the 
rainbow  ,i i.  "ii. •  ii'iiu  liquid  drops  with  n,  >  > p,  and 
mm  uicl  ir >i..  .■  e  i 1  nil  i  c,  'i  :oii.  The  author  lias  found"'  I  ha  I 
then  -  put  mg  I  , e  , eiselv  proport ional  ton  and  for  drops 
of  xyleii.  1  1  "ii  in  water:  A  ~  U.ti.  These  fringes  were 

found  i"h,  ua!  i  dete*  i  mg  micron-amplitude  shape  ox. 
cilia  I  um  T  i  ■  The  ■  i  <  i  ■ !  i  a  i  net  i  -  drops  since  changes  in  t  he 
drop’-,  uipi  i  i.i-i  m  ingiil.tr  shift  in  ill,  posit  ion  of  these 
Iriip ;i  ■"  •  "0  '"Hi  "i  F  when  the  shape  oftlm  drop  oscil¬ 

lated  a.  i  i  piii.  •  a  i  ot  a  fexv  II/..  the  fringes  appeared 
to  I ili 1 1  i. ■  '  I  ,  a  ul .  - ..  illations  ot  the  shape  ut  the  Imlilile 
should  .  •  ll"  me  si  nut  tire  lo  blur,  lialiid, 

Inihti!  !  I,..- .  iias  u  e-..  ii*-l  hv  room  vibrations  tnav  have 
bln i o  T  i is  i  ■  ■  a  tin.  r  negatives  xvhets  none  were  i  h 

-leixi  T  a  xg.  a  "a  a.  iv  t vpically  I ,'TO si.  I’hesiniilaiitii  ■■ 

ami  diM. . ''-in  li  angle  scattering  with  rainbow 

m  it  ■  i  1 1  ic  ,  .  i  .  ■  '  o  ('■.  -  anparing  I  ig.  -I  with  photogr.iph:- 

nt  t  In  1  a  I  I  "  I  "  I  .  I  ',  I  ...  I  III  I  i 


VI.  DISCUSSION 

Tin'  tiitn.il  ii  'Jr  sr.il  taring  of  eight  l  ill  hides  was  phot" 

graphed,  the  ’ n i I ih|*.  . . . a  .  were  in  the  range  Mi  '.!.T  nun 

F.h  Ii  ph"  i, a  iph  i  mill. ir  in  appearance  to  Fig.  I  in  that  die 
ring,  to  ;  "i  luli'i'  at  x  nor  In  I  lie  rl'it  it  a  l  angle  is  1  ptile  visible  hi 
cithci  Pile  it  im  id'  ol  p.'lal  i/at  ion.  I  'niiscqueut  l\  llltle  m 
I  ,  n  1 1 1 1  .  "  I  ■  "  1 1'  H  I '  I  t  lie  pi  l'd  let  ion  ol  geni net  I'll'  opt  II  - 
lira'  I  I  '  I  l  ;  ’  i  1 .  r  ■  ■ !(  I  "  r  .III  I  M  lllll  It's  ol  til  Is  si  ,'r  II I  XX  .1 1 1  I 


Mill 
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Furthermore,  the  measurements  of  the  positions  of  the  in 
tensity  maxima  and  minima,  measured  relative  lo  an  node 
lermined  angle,  are  m  reasonalile  agreement  with  K(|.  (Hill 
An  experimental  lest  ot  several  aspects  ul  this  calculation  mav 
require  the  use  ol  photoelect  lie  det  ei  I  <  ns  hecause  of  limit  a 
tiolis  in  the  present  appaiat  u  These  include  t  he  following 
features:  (i)  the  amplitude  of 'the  ringing  ot  the  intensity,  (n) 
the  ahsolutc  positions  in  the  maxima  and  minima,  (in)  the  rale 
of  decay  ot  the  intensity  for  ir  $  1,  ( iv )  the  rang  •  ol  radii  for 
which  K(|.  CIO)  may  he  used,  and  (v)  ipianlila"  la-nre 
ments  of  the  interference  phenomena  descrihe  .  i  sec.  \ 
Unfortunately,  the  largeness  of  the  hulihle  radii  makes  it 
impractical  to  calculate  the  Mie  scatterini:, 

A  brief  discussion  of  this  use  of  the  physical  optics  ap¬ 
proximation  is  in  order.  Surface  wave  phenomena1  do  not 
appear  to  he  included  in  this  calculation;  Kir;.  I>  is  evidence  that 
they  are  not  needed  tor  an  approximate  description  ol  tliesi 
observations.  Surface  waves  can  he  included  m  the  Watson 
transformat  ion  method  which  has  been  used  to  ! annulate  the 
scattering  from  a  cavity. Is  t"  A  comparison  with  the  results 
in  Fi>;.  f>  has  not  been  performed  because  of  t  lie  complicated 
nature  of  predicted  scattering  amplitude  jc  g  ltd  1H,  Kps. 
(4.9'd)  and  (■!.  HU  )J;  indeed,  the  present  author  h  as  been  unable 
to  determine  it  the  ringing  is  consistent  with  those  amplitudes. 
A  second  short  com  i  nr;  of  the  approximations  presented  in 
Secs.  I  III  is  that  the  curvature  of  the  hulihle  orthogonal  to 
the  plane  of  scaltci inn  is  omitted  until  Kp,  (HI)  where  its 
geometric  effect  is  included  in  the  a •'  factor. 

Diffraction  may  have  a  major  effect  on  the  far  field  ap 
pearance  of  the  crit leal  nn);le  scattering  for  hi i hides  in  water 
illuminated  liy  collimated  white  light.  I.et  v>,  and  </•,  "  denote 
respectively  i/\  tor  violet  and  for  red  light.  The  dispersion  ot 
water  and  Kps.  (I)and  (2)  give:  </>, '  -  <;■,  '  ~  17  tnrwl;  enn- 
sequently  the  geometrical  scattering  (Fig.  I  Hof  Kef.  1),  as  seen 
through  a  telescope,  should  appear  liluish  in  the  region  </*,. " 
<  0  <  0,.'.  It  is  evident,  however,  from  Fig.  5  that  diffraction 
smearing  of  the  otherwise  abrupt  change  in  intensity  tor 
monochromatic  illumination  will  tend  I,.  reduce  significantly 
th'  coloration  predicted  by  geometric  opt  ics 

This  work  was  motivated  hv  an  investigation  of  optical 
methods  for  measuring  the  refractive  index  of  superheated 
liquid  drops  suspended  in  immiscible  host  liquids.  At  stiff i 
cicntly  high  teinperal ores,  such  drops  frequently  have  n,  < 
n„  and  they  should  exhibit  critical  angle  scattering  phenomena 
similar  to  that  described  for  hubbies.  Kven  though  the  dis¬ 
tinctive  critical  angle  scattering  predicted  by  geometric  optics 
is  altered  by  diffraction,  n  may  he  iulerreii  from  the  absolute 


angles  of  the  intensity  maxima  and  minima,  provided  the 
present  theory  is  found  to  give  the  ubstilulr  ungles  when  n  is 
known. 
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Scattering  by  bubble  in  water  near  the  critical  angle:  inter¬ 
ference  effects  (P.  L.  Marston  and  D.  L.  Kingsbury)  0.  Opt. 
Soc.  Am.  71,  192-196  (1981)  [An  Erratum,  reproduced  below, 
was  published  in  0.  Opt.  Soc..  Am.  71,  917  (1981).  This 
notes  errors  of  transcription  whicFTare  not  corrected  in 
the  reprint  given  here.] 


errata 


In  the  caption  to  Fig.  3  of  our  paper,1  q  ■  -J>  -  </>,.  should  be 
written  */  =  </>,  -  <l>,  which  is  the  definition  of  t?  given  in  the  text 
of  our  paper.  In  the  line  below  Kcj.  (If)),  ii\  «■  (8/nn)l/2cos  0r 
should  be  written  ns /i i  =  (8/nn)l/l!/cosW,.  Also,  the  top  line 
of  the  right-hand  column  on  page  194  should  be  written  as 
", . .  error  from  using  Eq.  (22) . . . 
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A  physical -optics  approximation  is  derived  fur  light  scaticring  by  dielectric  spheres  with  refract  ive  indices  less  than 
their  surroundings,  and  it  is  applied  In  air  hubbies  in  water  The  approximation  gives  the  coarse  structure  in  the 
scattering  when  the  scattering  angle  </i  is  near  the  critical  scattering  am  le  i /», .  where  </>,.  H.l°  for  bubbles  in  water. 
Diffraction  has  been  observed  to  be  important  in  the  critical  region  because  of  an  abrupt  change  in  the  amplitude 
of  the  relied  ed  wave  |I*.  I..  Marston.  -I.  Opt.  Sue.  Am.  1*9,  1205  12 1  I  1 1  OHO)].  Interference  that  is  due  to  a  refracted 
wavt  produces  oscillations  in  the  intensity  with  an  angular  quasi-period  of  magnitude  I \/ti  I1  ’  rad  near  i/i,  .  where 
A  is  the  wavelength  and  ti  is  the  radius.  Unlike  diffraction  related  ost  illations,  the  interference  oscillations  in¬ 
crease  in  magnitude  as  forward  scattering  is  approached.  < fptical  tunneling  thrcaigh  bubbles  is  also  discussed. 


INTRODUCTION 

An  approximate  description  of  the  scattering  of  light  by  large 
dielectric  spheres  may  be  obtained  from  ray  optics,  provided 
that  the  observation  angle  is  not  in  a  region  in  which  wave¬ 
front  diffraction  is  dominant.1  In  the  ordinary  case,  in  which 
the  refractive  index  of  the  sphere  exceeds  that  of  the  sur¬ 
roundings,  forward  scattering,  the  rainbow,  and  the  glory  are 
the  regions  in  which  diffraction  is  important.1  :*  In  the  less 
explored  case  of  a  sphere  whose  refractive  index  is  less  than 
the  surroundings,  a  new  region  appears  in  which  diffraction 
is  important.'*  'Phis  region  may  be  called  the  critical  scat¬ 
tering  region  because  it  is  associated  with  light  reflected  from 
the  sphere's  surface  when  its  angle  of  incidence  is  close  to  the 
critical  angle  normally  associated  with  total  reflection  from 
plane  surfaces.  Diffraction  effects  in  this  region  have  been 
observed  and  described  with  a  physical-optics  approximation. 1 
The  purpose  of  this  paper  is  to  discuss  the  effect  of  the  in¬ 
terference  of  the  reflected  wave  with  certain  rays  that  pene¬ 
trate  the  sphere  and  are  refracted  to  angles  in  the  vicinity  of 
the  critical  scattering  region.  This  interference  leads  to  a 
modification  of  the  previous  model*  that  becomes  significant 
in  the  critical  region  when  an  air  bubble  in  water  has  a  radius 
<0.-1  mm. 

I.et  the  scattering  angle  </>  denote  the  angle  of  the  distant 
observation  point  measured  with  respect  to  the  direction  of 
propagation  of  the  incident  plane  wave  and  the  sphere's 
center.  The  critical  scattering  angle  is  i/>,  =  it  -  20, ,  when 
0,  =  arcsin  (n  ')  is  the  critical  angle  of  incidence  and  n  = 
in, ,/n,)  >  1,m„  and  n,  being,  re> pectively,  the  refractive  indices 
of  the  outer  and  the  inner  media.  Air  hubbies  in  water  have 
n  ~  MW  and  </>,  ~  82.82°.  Figure  1  illustrates  several  rays 
that  have  a  deviation  angle  <A  of  50°.  Each  ray  is  characterized 
by  a  parameter  p,  where  /j  —  1  is  the  number  of  reflections 
from  the  internal  surface  and  p  =  0  has  only  an  external  re¬ 
flection.  Uavs  that  enter  the  bubble  below  the  center  line  are 


indicated  with  a  prime.  The  ungle  of  incidence  measured 
from  the  surface  normal  is  II,,,  and  the  corresponding  refracted 
angle  is  it,,.  Rays  with  p  ?  0  have  0/,  <  0 , . 

The  previous  physical-optics  approximation'*  described  the 
diffract  ion  of  the  virtual  wave  front  associated  with  t  be  p  = 
0  ray  only.  It  was  shown  that  diffraction  is  important  in  the 
region  |  (A/o  I1'-’ rads,  where  n  is  the  sphere  radius 

and  A  is  the  wavelength  in  the  outer  media;  consequently, 
diffraction  is  important  near-/),  even  if  A/ct  is  quite  small.  The 
present  paper  shows  that  the  interference  that  is  due  to  the 
p  =  I  ray  leads  to  a  modulation  of  the  scattered  intensity  as 


o  ?  i 


Kig.  I  Hay  paths  in  the  scattering  plane  with  a  scattering  angle  </> 
=  The  number  gives  the  ray  parameter  />  The  dashed  ray  is 

forward  scattering  that  is  due  to  tunneling  from  the  reflected  ray  and 
is  discussed  in  Sec  I. 
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if  function  ot  </>.  The  /)  >  2  r.iys  will  lie  omitted  since  they  lend 
to  virtual  wave  fronts  much  smaller  in  amplitude  than  the  /> 
=  0  and  p  —  1  wave  fronts.  The  p  =  2'  ray  leads  to  a  super¬ 
posed  low-amplitude  modulation  of  the  scattered  intensity 
with  an  angular  period  ^A(,Vn),  where  a/ A  is  the  lateral 
separation  ot  rays  0  and  2’  after  emergence  from  the  sphere. 
For  an  air  bubble  in  water,  ray  tracing  gives  A  ~  0.82  with  (/> 
*  «/>, ,  whereas  observations  close  to  0 (Sec.  V  of  Ref.  4)  gave 
A  =  0.70  ±  0.04.  The  angular  period  of  this  modulation  is 
when  a  »  A.  The  emphasis  of  the  present  paper 
is  on  the  coarse  structure  in  the  scattering;  consequently, 
printed  ravs  are  omitted.  Mie- theory  results'1  (to  lie  published 
separately)  substantiate  these  approximations  in  the  region 
(' />,  —  40°)  <  0  <  0,  for  n  —  4/H. 

1.  DESCRIPTION  OF  THE  VIRTUAL  WAVE 
FRONTS 

Diffraction  and  interference  effects  in  the  scattering  are 
modeled  by  first  using  ray  outics  to  describe  virtual  wave 
fronts  positioned  within  the  si  atterer.  These  wave  fronts  are 
then  allowed  to  diffract  to  th*  fur  field.  Fortunately,  it  is  not 
necessary  to  express  the  exact  amplitude  and  phase  of  the 
wave  fronts  in  terms  of  Cartesian  coordinates.  This  is  because 
approximate  dependences,  dciived  below,  justify  the  use  of 
the  stationary-phase  upproxin  a. ion  in  the  diffraction  inte¬ 
grals. 

Figure  2  illustrates  Cartesian  coordinates  used  in  the  de¬ 
scription  of  virtual  wave  fronts  in  the  scattering  plane.  The 
origin  is  chosen  to  be  the  point  where  the  reflected  ray  with 
f>»  -  0,  contacts  the  sphere.  The  positive  e  axis  is  the  path 
of  the  scattered  critical  ray  predicted  hv  ruy  optics.  The 
positions  of  the  virtual  wave  fronts  associated  with  rays  0  and 
I  are  denoted  by  the  functions  and  iqtu),  where  the  u 
axis  is  perpendicular  to  the  t1  axis  und  lies  in  the  scattering 
plane.  Phase  shifts  of  ray  0  that  are  due  to  reflection  beyond 
the  critical  angle  are  not  included  in  r„  but  are  accounted  for 
with  a  separate  phase  factor  in  the  description  of  the  ampli¬ 
tude.  The  reference  phase  is  selected  such  that  r„(0)  =  (). 
Ray  tracing  gives  in  =  r* ■  at  u  =  t),  and  the  leading  term  in  a 
Taylor  series  for  in  is1 

I'u  =  <ru  ",  (I) 

where  o  =  -In  cos  0,  )  The  Taylor  series  for  iq  may  he 
obtained  by  including  effects  of  ret  ract  ion  and  t  he  change  in 
velocity  within  the  sphere;  the  leading  term  is 

r  i  =  'hiui'.  (2) 

» 

The  resulting  wave  fronts  are  shown  in  Fig.  2.  The  calculation 
ol  the  interference  of  waves  0  and  I  uses  the  exac*  path-length 
difference  that  is  not  easily  expressed  as  a  function  of  u. 
Effects  of  wave-front  curvature  orthogonal  to  the  scattering 
plane  will  he  approximated  with  normalization  factors  to  he 
introduced  subsequently 

In  the  description  of  the  scattering,  it  is  only  necessary  to 
consider  two  cases  with  orthogonal  incident  polarizations, 
since  the  general  case  may  lie  obtained  by  a  linear  combina¬ 
tion.1  The  subscript, ;  is  assigned  the  value  I  when  the  inci¬ 
dent  electric  vector  is  entirely  perpendicular  to  the  scattering 
plane  and  the  value  2  when  entirely  parallel.  Waves  scattered 
to  the  tar  field  associated  with  p  =  0  and  p  =  1  rays  may  be 
simulated  by  a  wave  on  the  u  axis  with  the  following  ampli¬ 
tude: 


v 


Fig.  2.  Positions  of  the  virtual  wave  fronts  for  p  =  It  and  p  ~  1 
si  uttering.  The  dashed  rnv  is  the  incident  critical  ray  that  is  lellected 
to  become  the  e  axis.  A  portion  of  the  bubble's  surface  is  also 
shown. 


Up./Ut.t )  =  h,,j  exp|t(-  hvr  -  u)( )],  if) 

where  t  -  v  — 1  ,  w  is  the  frequency,  f  is  the  time,  und  h  = 
2tr/\.  The  exp(-  iuit )  factor  will  be  omitted  subsequently. 
Ray  tracing  and  the  physical -optics  approximation  give  the 
following  amplitude  factors1-1: 


ho j  =■ r, {»»),  u  <  0,  (4) 

/iiu  =  it  >  o.  tr>> 

/l|.,  =  0,  If  <  0,  Hi) 

h\,j  -  1 1  —  fy(tf|)”|  \  kl)  ,  u  >  0,  (7) 


where  the  r,  are  Fresnel’s  amplitude  reflection  coefficients 
for  a  plane  surface, 11 


r\(t),,)  =  : 


_  sintp,,  -  Op) 
sinlp,,  +  Op ) 

r,U)r)  =  , 

tantp,,  +  Op) 

and  Pp  is  the  refraction  angle  predicted  by  Snell’s  law. 


18) 

(!)) 


Pp  =  arcsintn  sin  Op).  (Ill) 

The  Stokes  symmetry  relations  have  been  used  in  the  deri¬ 
vation  of  Kq.  (7),  where  I)  is  Van  de  Holst’s  divergence  factor,1 
which  may  be  written’ 


I!  _  _ _sin  0 1  cos  0 1 

2 1 1  -  (n  cos  ff|/cos  pj)jsin  0 

0  is  the  scattering  angle  of  the  p  =  1  ray  predicted  hyruyop- 
t  ics, 


0  =  2(/)|  -  (f|l,  112) 

and  iv  is  a  positive  constant  normalization  factor  determined 
subsequently. 

I  he  phase  shi! tsM  ,1,  are  associated  with  the  (nearly)  total 
reflection  of  ray  0  as  If,,  exceeds  0,  (tunneling  through  large 
bubbles  can  lie  neglected;  see  Sec,  4): 

tnn(<V2)  =  (sin- n  - -)>  J/cos(l,„  lilt) 

tan(,V,/2)  =  n  -  tan(<i,/2).  ,  pp 

'i  he  !),  vanish  when  «„  =  II , ,  hut  they  become  positive  as  0„ 
exceeds  0, ,  which  represents  an  (ir/rannmrnl  of  the  p  =  0 
virtual  wave  front.  (The  reader  is  cautioned  that  some  an 
I  hors  write  A i  and  <Y.  with  an  incorrect  sign  for  their  choice  of 
the  sign  ol  w  in  the  time-dependence  factor/1) 

Approximate  dependence  of  hJK,  on  small  values  of  u  is 
obtained  by  first  relating  If,,  top  foru  ~  (I and  It,,  ~  II. .  This 
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leads  to  the  following  proportionalities:  1  —  ht))  a  (— u/a)tei 
for  u  <  0,  ,'i  i  j  u  (u/fj)  '  -  for  u  >  0,  and  1 

bj  pju"*,  (If.) 

where  f1\  =  (H/nti )'  •'  cos  ft,  and  fi,  =  n-fi |. 

2.  FAR-FIELD  DIFFRACTION  AND 
INTERFERENCE  OF  THE  SCATTERED  WAVES 

Let  R  denote  the  distance  to  the  observation  point  from  the 
center  of  the  bubble.  When  R  »  kn-,  the  scattering  ampli¬ 
tude  is  proportional  to  the  Fourier  transform  of  (/„,  +  U\ 
which  may  be  written  as  /,  =  /,,.,  4-  f\  ,,  where 

/,,,>(sin  7)  =  J'  Up.jiu)  expt-  iku  sin  7>dt/  (16) 

and  7  =  0,—  0.  An  approximation  for/,, ,  with  7  >  0,  which 
was  introduced  in  lief.  4,  is  to  let  =  0  for  u  <  0  because 
l  hyj  decreases  abruptly  as  u  becomes  negative.  It  is  not  clear 
if  ray  optics  can  be  used  to  model  the  virtual  wave  in  the  region 
in  which  \Hh{)j/ilu  |  is  large;  however,  its  use  leads  to  the  dif¬ 
fraction  of  reflected  energy  away  from  the  critical  region  and 
supports  the  approximation  U»,j(u  <  0)  -  '•  Most  of  the 
diffracted  p  =  0  wave  in  the  critical  region  comes  from  the 
virtual  wave,  where  u  is  small  and  positive,  and  an  approxi¬ 
mate  result  is  obtained  by  using  Eqs.  (1 )  and  (15)  in  the  entire 
positive  u  domain.  The  integral  may  be  approximated  bv 
using  the  principle  of  stationary  phase1"  with  the  result  that 
/,,  ,  =  (A a  cos  (lc),/2 go, j/2,  where4-1 1 

go  j  ^  expu'7„,;)  |F(u.)  -  /•’(-  oo )  |,  (17) 

w  =  [(a/A)cos  0,  |I/,J  sin  ij,  (18) 

7 nj  -  - cos 0,.)u*  +  — - — •  (19) 

4  a 

wheic  the  first  term  of  Eq.  (19)  is  the  phase  shift  of  the  geo¬ 
metric  ray  reflected  to  </>  and  l'(u)  is  Fresnel's  integral, 

F(w)  =  ^  exp[i(7r/2)2-|dr.  (20) 

When  the  observation  point  is  in  the  critical  region,  Eq  (2) 
may  lie  used  in  the  evaluation  of  [\ y. 


wave.  Consequently,  the  error  from  using  Eq.  (28)  in  the  total 
srat taring  should  be  acceptable  even  if  7  is  small. 

The  scattered  intensity  normalized  to  that  of  the  geometric 
result  for  a  perfectly  reflecting  sphere  of  the  same  size  is 

lj  =  'h  U’oj  +  gt,;|7  (24) 

where  the  normalization  was  obtained  by  noting  that  |g,,j  |  ► 

\  2  as  ic  ■■»<».  The  actual  intensity  is  t  he  incident  intensity 
multiplied  hv7  I,(a/R)-/4.  The  value  of  s  in  Eq.  (7)  is  deter¬ 
mined  to  be  2  by  noting  that  only  then  doesg,;  become  the 
value  predicted  by  ray  optics.7  Because  of  the  interference 
of  the  Hp  j,  the  intensity  is  sensitive  to  the  computation  of  <7oj 
-  71,).  Comparisons  with  Mie  theory  have  shown  that  the 
phase  difference  given  by  Eqs.  (19)  and  (28)  is  limited  in  its 
usefulness  because  of  errors  introduced  by  Eqs.  (1 ),  (2),  and 
(15).  The  agreement  with  Mie  theory  is  significantly  im¬ 
proved  when  ( 7,,.y  —  7 1  ,y )  is  taken  to  he  the  exact  value  pre¬ 
dicted  by  ray  optics17: 

7o  j  ~  7i  j-  2/m(cos()|  -  cost/,,  -  n  ~ 1  cos/q)  -  A, (()„), 

(25) 

where  (),,  is  predicted  bv  ray  optics  to  be 

(>„  =  (7T  —  0)/2  -  II,  +  (77/2).  (26) 

To  evaluate  Eqs.  (22)  and  (25)  for  a  given  value  of  0,  it  is 
necessary  to  solv’  the  transcendental  Eqs.  (10)  and  (12)  for 
II 1  and  (>\.  The  evaluation  of  Eq,  (24)  requires  the  evaluation 
of  Eqs.  (7)- (14),  (17),  (18),  (22),  (25),  and  (26). 


3.  DISCUSSION 

Figures  8-6  illustrate  the  result  of  three  approximations  for 
the  scattered  intensity  for  an  air  bubble  in  water  normalized 
to  the  ray-optics  value  for  a  perfectly  reflecting  sphere.  In 
each  figure,  n  =  4/8.  The  ray-optics  result  for  a  bubble  is 
equal  to  1  +  2.*i  1,/-’  when  7  >  0  and  rj  (((„)- when  7  <  0.  It  ex¬ 
hibits  a  cusp  at  the  critical  angle,  which  is  also  evident  in  the 
rav-optics  predictions  plotted  by  Davis7  and  by  Welford.1- 
The  predictions  for  the  diffraction-only  model  developed  in 
Kef.  4  are  given  by  |g„  ,  |  -72.  The  present  model  must  even¬ 
tually  break  down  in  the  region  where  7  <  0  because  of  the 


fij  ~  h]j  exp|-  i\ku  sin  7  +  (A<vu'-72)||du.  (21) 

The  phase  of  the  complex  exponential  is  stationary  when  u 
=  i7,  where  u  ~  a  cos  I),  sin  7.  Since  hij  is  a  slowly  varying 
function  of  u,  the  stationary-phase  approximation  may  be 
used  to  evaluate  Eq.  (21)  provided  that  U  is  well  within  the 
domain  of  integration.1"  The  result  is  / ,  ,  =  (An  cos  0,  ),r~ 
g  1  where 


1  (  7T  \j  , 

g  1  ,j  -  2h  1. ,((),)  exp 

17  +  71;  +  0(A/a)'-\ 

1  \4  ) 

(22) 

Msin  7)- 

Tt ./  = 

2a 

(28) 

In  Eq.  (22)  the  stationary-phase  point  that  is  used  in  the 
evaluation  of  h\  t  corresponds  to  the  II 1  given  by  Eqs.  (10)  and 
(12).  The  use  of  the  stationary -phase  approximation  breaks 
down  as  7  -  »  0,  since  t7  0;  however,  this  is  just  where  It, ,  * 


7  DEGREES 

Fig.  8.  Three  approximations  for  scattered  intensity  is  a  function 
ol  7  =  i/i  —  1 />,.  with  n  =  4/8,  the  electric  field  perpendicular  to  the 
scattering  plane,  and  An  =  10,000.  The  dashed  curve  is  the  result  of 
rav  optics,  which  adds  the  intensities  of  p  =  0  and  p  =  1  scattering. 
The  dotted  curve  is  the  diffracted  p  =  0  wave  only  and  is  from  Kq. 
1 17).  The  solid  curve  is  given  by  Kq.  (241,  which  includes  the  dil- 
lraction  oft  he  p  =  0  wave  and  the  interference  of  the  p  =  1  wave. 
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Kit;.  4.  Siune  ns  Kit;,  It  but  with  parallel  polarization. 


'  .'ll  to 

r,  OFORKf  5 

Kit;.  •>.  Same  as  Kit;.  r>  hot  .  it  h  parallel  polarization. 

assumption  that  f  hr  <  hi  =  0.  In  the  figures,  the  model 
hits  been  extended  into  t hat  region  by  evaluating  Kip  KM)  with 
g,  ,(t]  <  0)  =  I).  Fresnel’s  integral  F(al  was  evaluated  with 
a  numerieal  approximation.1 1 

The  bubble  sizes  in  these  figures  are  manifest  by  consul 
ering  illumination  by  a  He  Ne  laser  with  a  wavelength  in  air 
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of  (iH2.H  nnt  and  a  wavelength  in  water  X  ~  474.fi  tint.  Then 
ha  =  10  (MM)  gives  a  =  0.755  mm.  When  measuring  the  scat¬ 
tering,  refraction  corrections  may  be  required  to  relate  the 
scattering  angle  in  the  water  to  that  seen  by  an  external  in 
st  rumen!. 4 

The  diffraction-only  theory4  predicted  oscillations  in  the 
intensity  that  are  homologous  with  respect  to  bubble  size, 
polarization,  and  refractive  index.  They  decrease  tnonn- 
tonically  in  magnitude  as  r/  increases  and  have  an  angular 
quasi-period  £(X/u)1/a  rad.  Figure  H  shows  that,  for  j  =  1  and 
ha  =  10,000,  interference  alters  the  scattering  significantly 
for  Z  10°.  After  the  fourth  maximum,  the  oscillations  in- 
crease  in  magnitude.  It  is  noteworthy  that  the  quasi-period 
of  the  newly  predicted  oscillation  is  also  roughly  (X/n  )l/i!  rad 
|as  may  he  shown  by  Kqs.  (1)  and  (2)|.  Figures  4  -0  show  that 
the  effect  of  interference  is  even  larger  for  j  =  2  and  Ibr  ha  ~ 
1000.  'The  intensity  averaged  over  a  quasi-period  is  roughly 
the  ray -optics  result.  Consequently,  for  large  spherical 
bubbles  in  white  light,  the  ray-optics  results  are  useful  except 
a  lien  to  S  2,f>,  where  diffraction  corrections  become  essential 
(ic  =  2.5  corresponds  to  i]  =  4.4°  at  ha  =  10, (KM)  and  t)  =  14.0° 
at  ha  =  1000.) 

The  experimental  data  plotted  in  Ref.  4  were  for  ha  *  10 
950  and  j  =  2.  A  comparison  of  Fig.  4  with  those  data  indi¬ 
cates  that  it  is  not  surprising  that  the  diffraction-only  model 
correctly  described  the  relative  positions  of  the  first  four 
maxima  hut  that  it  underestimated  the  magnitude  of  the  in¬ 
tensity  oscillations  for  i)  >  11°.  A  new  apparatus  has  iteen  built 
that  permits  observations  in  the  region  r;  >  7°,  which  was 
obscured  in  the  original4  apparatus.  Visual  observations  for 
air  bubbles  with  a  ~  1  mm  in  monochromatic  light  suggest 
that  the  magnitude  of  the  intensity  oscillations  tends  to  in¬ 
crease  with  increasing  n  hut.  as  of  the  time  of  this  writing, 
photographs  of  a  quality  to  permit  quantitative  measurements 
have  not  been  obtained.  These  observations  preceded  and 
motivated  the  present  theoretical  effort. 

The  intensities  given  by  Kq.  (24)  have  been  compared  with 
Mie  theory1  14  with  n  =  4/!l  and  ha  =  25,  100,  KWH),  5000,  and 
10, (MM).  The  computations  were  performed  using  Wiscombe’s 
Mie-scnttering  algorithms1’’ with  minor  modifications.  The 
eomparision  shows’’  that,  for  0  <  rj  <  00°,  the  course  structure 
is  approximately  described  by  Kq.  (24).  Krrors  in  ihe  pre¬ 
dicted  locutions  of  intensity  oscillations  are  typically  less  than 
one  tenth  of  a  quasi-period.  Furthermore,  the  phase  differ¬ 
ence  of  the  scattering  amplitudes  for  the  two  polarizations 
given  by  Kq.  (28)  (below)  reproduces  the  coarse  structure  in 
the  Mie  results.  Kquntion  (24)  differs  most  from  the  coarse 
structure  in  the  Mie  results  when  j  =  2.  This  may  be  due  in 
part  to  its  omission  of  the  p  =  2  wave,  which  has  a  smaller 
amplitude  for  j  =  1  than  for  j  =  2.  For  example,  at  )/  =  110° 
the  intensities  of  p  =  2,  scattering  in  the  units  of  Figs.  II  (i,  are 
given  by  ray  optics’  to  be  0.012  and  0.015,  respectively. 

Kquat  ion  (24)  does  not  describe  fine  struct  are  in  the  Mie 
intensities  which  is  most  prominent  when  j  —  I.  This  si  rue 
lure  has  a  quasi  period  ~A(X/n)rad  (0.11°  for  ha  ~  1000), 
which  indicates  that  it  could  be  included  in  the  model  with  the 
addition  of/)  =  2'  scattering.  Ray  optics  verifies  that  the 
amplitude  of  2'  scattering  in  the  critical  region  is  largest  when 
/  =  1.  Kxcept  for  the  tine  struct ure,  tlnf  Mie  intensities  de 
crease  gradually  as  r;  decreases  I  hrough  zero.  The  error  in  Kq. 
(24),  with  )/  negative,  generally  decreases  with  ha  as  diffract  ion 
becomes  more  important . 1 
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Equation  (24)  may  Ik*  converted  to  the  units  commonly  used 
in  Mie  algorithms  through  the  following  transformation: 

=  (A-ul-7,/1,  (27) 

where  .S’,  are  the  scattering  amplitudes  in  the  usual  units.11, 1,1 
The  phase  difference  of  .S' i  and  .S' .  is  needed  if  the  scattering 
for  arbitrary  states  of  the  incident  polarization  is  to  lie  cal¬ 
culated. 1111  The  phase  difference  is  predicted  to  he1 ' 

nrg(.S\d  -  arg(.S'i)  =  argtA’n..'  +  am..')  -  urgtA’n.i  +  am  |), 

(2H) 

where  nrg  denotes  the  complex  argument  of  the  quantity  in¬ 
dicated.  The  present  calculation  uses  an  expt—  lull )  time 
factor,  whereas  the  ,S,  are  often,11'’  hut  not  always,11’ specified 
by  using  ex|>(+  i cof ),  so  the  right -hand  side  of  Kip  (28)  should 
he  multiplied  by  —  1  for  some  applical  ions. 

The  coarse  structure  in  the  scattering,  which  arises  because 
of  the  interference  of  the  reflected  and  refracted  light,  may 
he  a  useful  experimental  tool.  The  angular  position  of  the 
maxima  should  shift  if  there  are  changes  in  the  size  of  a  bubble 
or,  for  the  case  of  a  drop  with  n,  <  ri„,  if  the  refractive-index 
ratio  changes.  Shifts  in  the  structure  of  rainbow  scattering 
have  been  useful  for  detecting  shape  oscillations  in  drops,11* 
and  a  similar  technique  could  be  used  for  detecting  small 
changes  in  the  scattering  from  hubbies. 

Diffraction  associated  with  the  critical  angle  and  structured 
scattering  that  is  due  to  interference  should  also  he  present 
in  the  scattering  by  cylindrical  bubbles  or  dielectric  cylinders 
with  n,  <  n„.  for  a  circular  cylinder  with  a  symmetry  axis 
perpendicular  to  the  incident  beam,  the  critical  scattering 
angle  and  ray  tracing  in  the  sealterinp  plane  are  the  same  as 
those  described  here  for  spheres.  The  divergence  factors  used 
in  the  description  of  the  virtual  wave  fronts  differ  from  those 
of  a  sphere.  The  geometric  optics  of  a  large,  perfectly  re¬ 
flecting  cylinder111  show  that  (neglecting  tunneling)  |hn.,(n  > 
0)  |  is  no  longer  constant  but  is  "  |siti(«/>/‘2)| 1  The  normalized 
magnitude  and  positions  of  the  intensity  oscillations  should 
differ  from  those  of  a  sphere.  Critical  angle  diffraction  and 
interference  structure  should  also  be  present  in  acoustic 
scattering  from  fluid  spheres  and  cylinders. 

4.  ELECTROMAGNETIC  TUNNELING 
THROUGH  SPHERICAL  BUBBLES 

The  present  model  hits  assumed  that  the  reflection  is  total 
when  II  ii  >  II, .  The  existence  of  evanescent  waves  within  t  he 
bubble  will,  however,  frustrate  the  total  reflection,  as  in  the 
case  of  plane  dielectric  slabs  separated  by  an  air  gap.  ’’  The 
effect  of  tunneling  on  nearly  total  reflection  at  curved  di- 
elect  ric  interfaces  has  been  studied  to  determine  losses  from 
curved  optical  waveguides”;  however,  the  sense  of  the  cur 
vat  lire  differs  from  that  in  the  present  case.  The  effect  of 
tunneling  through  bubbles  can  be  estimated  by  using  the 
plane  slab  results.  1  in  which  the  gap  <1  is  equal  to  the  bubble 
diameter  and  II,,  becomes  the  angle  ol  incidence  at  t  he  plane 
surface.  The  result  of  this  approximation  is  that  the  wave 
that  has  tunneled  through  the  bubble  will  always  be  trans¬ 
mitted  in  the  forward  direction.  The  dashed  ray  in  fig.  I 
shows  the  rav  transmitted  because  of  the  trust  rated  reflection 
ol  I  he  / i  =  ll  rav.  For  each  p  =  tt  ray  with  II,,  >  H, .  t here  is  a 
forward  transmitted  ray  with  r  =  II,,.  The  intensity  Irans 
mil  lance  I  hroiigh  a  plane  gap  is  found  to  he  of  order  (/.•</) 
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when  II,,  =  II,,  and  it  decreases  rapidly  as  II,,  exceeds  II, 
Tunneling  should  have  a  negligible  effect  on  the  reflection 
from  a  bubble  unless  ka  is  small.  This  conclusion  is  also 
supported  by  the  success  ol  the  present  model  in  reproducing 
Mie  results'’  when  ka  =  2,r>. 
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Mie-scnllering  algorithms  wart1  list'll  In  Compute  scattered  inti'iisitii's  and  phase  diiicrences  lor  air  lailiiili's  in 
water,  ltesnlts  arc  pli  >1 1  <*ci  as  a  function  ol  liic  scattering  angle  i/»  in  I  lie  general  range  "I  ill*  ttU°  lor  si/e  paraineters 
1,0  of '.id.  Hill.  1.IHI0.  and  III. (Kill  Icnrrcspomling  lo  radii  a  '  l.d  am  In  II. si  mini.  As  decreases  lielow  the  critical 
scat  ter  inn  angle  at  Mii.H".  I  lie  intensity  increases  and  undergoes  liroadly  spaced  osci  Hat  ions  l  lull  are  descrilieil  tiy 
a  physical-optics  approximation  developed  in  tile  aulliors  in  a  separate  puli1ic.il mu  Mie  scattering  also  exliilnls 
finely  spaced  oscillal ions 


INTRODUCTION 


( fbservnl  ions  nf  light  scat  I  eretl  by  nil'  Imlihles  in  water  in  I  lie 
vicinity  i it'  the  ci'it  ical  scattering  angle1  indicate  I  hat  dil'frai 
t  inn  is  significant,  even  it'  I  he  Im  hide  radius  a  is  as  large  as  I 
mm.  The  cuarse  st  met  tire  in  the  critical  scattering  region  has 
been  described  with  a  physical-optics  approximation  that 
accounts  tor  I  lie  dill  ruction  of  light  relied  ed  from  t  he  bubble's 
surface  and  the  interference  ol  certain  refracted  light.''  In  this 
Letter  we  use  Mie  theory  1  •’  to  calculate  tar  field  scattered 
intensities  and  phase  differences  in  the  critical  region  of 
spherical  bubbles  in  water  and  compare  these  exact  results 
with  those  of  t  he  physical  optics  approximation. 1  Previous 
applicat  ions  of  Mie  t  henry  to  hubbies  emphasized  t  he  caleu 
hit  ion  ot  the  scattering  efficiency  factor'1  *  and  t  lie  radiation 
pressure, H  In  contrast  with  the  scattering  by  spherical 
drops,1'"11'  the  exact  angular  structure  ir  the  scattering  by 
bubbles  is  not  wei!  explored. 

Consideration  of  the  approximate  scattering  '  introduces 
the  features  to  he  expected  in  I  he  exact  Mie  results  and  fa 
cilitates  the  choice  of  the  range  and  frequency  of  the  angles 
needed  for  a  descript  ion  of  t  he  angular  st  ri tel  til  e.  Figure  I 
i  I  lust  rates  several  rays  with  a  scat  taring  angle  i/i  of  50°.  In  t  his 
I  igure  and  the  Mie  coinput  at  ions  I  hat  follow,  the  relative  re 
I  ract  ivc  index  n,  In,,  is  taken  to  lie  ii/'-l,  where  n,,  is  t lie  refi  ne 
live  index  of  the  water  and  n,  is  the  ri  tractive  index  of  the 
bubble  contents.  Hays  are  characterized  by  a  parameter  /j. 
where  p  —  I  is  t  he  number  of  reflections  within  t  he  huhhle  and 
p  =  I)  has  only  an  external  reflection:  It,,  is  the  angle  ol  inci 
deuce  at  t  he  /Hit  ray  when  it  lirst  touches  t  he  surface  of  I  he 
bubble.  When  sin  II  ”  //,  in,,,  t  he  p  =  O  ray  is  totally  reflected 
I  till  effects  of  tunneling  and  surface  curvature  are  neglected  I 
The  corresponding  condition  on  the  scattering  angle  ol  the 
reflected  ray  is  "  H'_>. X'J".  where  ,  t  he  cril  ical  scat 

ter  lug  angle,  is  related  lo  n,  in,,  through  Kps.  I  I )  and  I  21  of  Kel. 
I .  I  For  air  hubbies  in  seawater.  •;»,  may  he  as  large  as  S  t".l 
I  lie  physical  opt  ies  approx i mat  ions 1  lead  to  l  he  predict  ions 
ol  st  ructure  in  t  he  angular  scat  lering  wit  h  an  angola i  spread 


of  the  order  of  1  A/a  1 1  rad  |  ~  I  Tit  An  I  1  <leg|.  where  A  is  t  he 
wavelength  m  thr  uutcr  mriliitm  and  A  =  'Jir/A.  Structure 
with  a  <|  i  las  i- pc  r  iod  of  that  magnitude  will  he  referred  tons 
coarse  slnuiurr,  in  coni  cast  to  /me  si  nut  it  re  in  the  scattering 
that  has  a  quasi  period  in  t  he  crit  ical  region  of  /ess  llnm  A/a 
rad  (or  H(i()°/Au  I.  The  f  ine  st rttel ure  original es  in  part  from 
I  lie  interlerence  of  I  be  p  =  2  ray  with  I  lie  rest  ol  the  scattering. 
To  describe  the  general  effects  of  the  critical  angle  oil  the 
scattering,  Mie  coin  put  at  ions  are  needed  wit  h  a  spread  i  it  | 

-  i/j|  somewhat  larger  than  (A/a  I1  with  angular  steps  some 
what  smaller  than  A  '2a  to  prevent  sampling  errors. 


COMPUTER-PROGRAM  CONSIDERATIONS 


The  Mie  solutions  were  computed  by  using  the  i-'uutuan 
Ml KVO  subroutine  developed  by  W  iscoinhe. '  1 1  chosen  in 


l‘  ij:  I  Kay  pat  hs  in  tin1  >c;il  Icrin^  plane  with  ;i  scat  tcrinj:  anjde  if* 

•  '<»'  I  lie  number adiarent  tn ea<  It  ray  i:ives  tin*  rav  parameter  /> 

I  In-  apprnximal  n>n  pM'ii  in  Ki  t.  includes  the  ini eileri'ine  and 
Hit  t  rat  !  a  hi  i>l  v\  ,i  \  i  s  ;t  'St  trial  r<  I  with  ravs  n  and  I 
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preference  to  *\is  lust  it  M1KV  I  mill  i  he  because  nl  I  he  inemorv 
limitatiims  of  the  HI’  MMX)  computer  sx.Uein  used  in  the  work. 
We  found  it  necessary  lo  romrrl  nil  Wiscnmhe’s  internal 
vnrinliles  lo  III1  doiilile  preeision  ( I  I  sigmlicanl  digits).  and 
we  were  I  hen  nhle  In  duplicate  his  'allies11  lor  scattering  hy 
a  drop  with  n, In,.  =  1.5.  iWisconihe's  suhrottlines  assume 
I  I  (lit; it  single  preeision.  I  Most  complex  aril h met ie  was  re 
writ  ten  as  real  aril  hniet  ie  hei  utse  of  I  he  ahsenee  of  eon  ;<le\ 
double  preeision.  The  logarithmic  derivative  of  the  lit 
ealti  Bessel  fulict ion,  commonly  denoted  l>v  wtiseaho 
lated  hy  downward  recursion,  which  is  known  lo  he  stable. 
The  program  was  tested  against  Wiseomhc's  tables  for  //,//'.. 
=  l.f>  of  the  complex  sealiering  amplitudes  and  effiil.ncy 
factors  and  against  the  scattering  elTieiem  y  factors  for  snta.'t 
size  parameters  when  n,/n„  ~  0.7.r>  as  ;>  veil  m  Refs.  7  and  S. 

Bel  ,V  denote  the  mini  her  of  terms  included  in  I  he  ennipu 
tilt  ion  of  the  Mie  series  and  /V,,,,,,  denote  the  tut  in  her  of  angles 
for  which  the  Mie  scattering  is  computed.  Since  MIFVt) 
uses''1'1 1  N  =  Att  +  f  ,05t  An )'  -'  +  2  and  /V,m|.  must  also  lie  in 
creased  in  |>roportion  to  A  a  to  resolve  the  fine  strut  e,  exe¬ 
cution  time  is  roughly  proportional  to  (An ) '  for  a  fixed  spread 
of  angles,  Whereas  the  results  for  An  =  25  rei|itired  a  fraction 
ol  a  minute,  those  for  At/  =  ID.IHIU,  whete  the  .ingle  step  size 
was  (l.(>(>5°.  retplired  2.5  h  for  HtMHl  angles.  I  These  times  were 
realized  alter  vectori/.ation  of  two  loops,  which  speeded 
ti)i  execttl  ion  hy  a  factor  of  approximately  I.) 

Higher  values  of  .V  than  that  given  in  Wiscottihe’s  equal nai 
above  were  tested  for  the  bubble  ease  with  An  =  11)1)  and  1000 
and  found  to  yield  identical  results  to  six  decimal  places. 

In  contrast,  the  physical -optics  approximation  program  was 
executed  in  seconds,  regardless  of  size  parameter.  Kquatioiis 
used  in  that  program  are  listed  at  the  end  of  Section  II  of  Hof. 
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( ‘nictitated  intensities  /,  as  a  fund  ion  of  t  lie  scattering  angle 
i ;>  are  shown  in  Figs.  2  I.  The  normalization  was  chosen  so 
that  /,(.;>)  =  I  represents  perfect  reflection  according  to  geo 
metric  optics.' 1 '  At  a  distance  /i‘  »  An  '  from  the  center  of 
the  sphere,  the  actual  /  polarized  intensilv  is  the  incident 
/  polarized  intensity  nniltiplied  hy  /, t d //>’ I  7 1.  where  /  =  I 
denotes  polarization  of  the  electric  vector  perpendicular  to 
t  he  scat  ter ing  plane  and  /  =  2  t  he  parallel  case.  The  van  de 
1 1 u  1st  normalization1  used  in  Wiscninbe's  routines  required 
that  the  squared  modulus  of  the  scattering  amplitude  com 
pitted  with  MIKVt)  he  multiplied  hy  (2/Aul  '  to  he  expressed 
in  I  he  units  of  Figs.  2  I.  ’ 

As  is  shown  in  the  figures,  t he  exact  solution  suggests  that 
a  fine  ripple  sl  met  lire  is  superposed  upon  the  coarse  structure 
that  is  described  hy  the  physical-optics  approximation.  The 
ripples  have  a  quasi  period  of  about  0.H2  l  \/<t  I  rad  near  t /■, . 
decreasing  slightly  with  ■;>.  To  avoid  severe  distortion  re 
stilling  from  sampling  errors,  il  was  necessary  to  use  till  angle 
step  size  less  than  <V7«  rad.  This  fine  st  met  lire  is  die  to  t  he 
interference  of  ravs  not  included  in  the  approximation,  pri 
tnarilv  the  //  =  2'  ray.  Consequently,  the  ripples  should  he 
damped  out  for  hubbies  containing  a  gas  that  absorbs  light. 
Fine  structure  from  air  bubbles  hits  been  observed.1 

The  plots  in  Fig.  5  show  !  lie  phase  difference  of  the  scat 
tered  waves 


a.'itt 


(b) 


Fig  2.  Calculated  normalized  scattering  intensities  for  Aa  =  in.lMHl, 
The  electric  vector  is  parallel  lo  the  scattering  plane.  The  solid  curve 
is  trout  Mie  tlioorv  The  dash  curve  is  the  phvsical  optics  ap 
proximal  nm  go en  in  Met  2. 

a  =  arglN  d  -  argt.S’ | ),  1 1 1 

win  re  the  S,  are  Mie  scattering  amplitudes  for  a  time  tie 
peiidcnee  of  expt  -tu/ ).  where  a'  is  t he  frequency  and  /  t he 
t  ime.  This  phase  difference  can  he  used  along  with  / 1  and  / 
to  compute  the  scattered  intensity  and  polarization  for  arhi 
tiiiry  states  of  incident  polarization,1  The  physical -optics 
approximation  Ibr  a  is  given  by  F.q.  I2H)  of  Kef.  2.  As  is  de¬ 
scribed  in  Kef.  2,  it  was  necessary  to  reverse  the  sign  of  a  front 
I  lie  output  of  t  he  M IF. VI I  program  because  Wiscombe  assumes 
a  time  dependence  of  expt  t  iu't  I. 


DISCUSSION 

For  ■;>  >  <;>, .  the  agreement  with  the  approx i mat  ion  is  best  tit 
low  size  parameters  where  dill  fact  ion  accounts  for  the  scat 
taring  in  this  region  in  the  approximation. 

The  coarse  structure  oscillations  at  angles  <  «/•,  are  similar 
to  the  predict  ions  of  diffraction  theory,  which  gives  a  relied  ed 
intensity  proportional  to  t he  squared  inodtl lus  of  Fresnel  iti 
tcgrals  |K,q.  CIO)  of  Ref.  l|.  The  first  maximum  of  the  re 
fleeted  intensity  occurs  at 1 

-  </.  -  1 .21  \/.t  cos  //,  1 1  1  Irttd).  (2) 

T  he  t'ornii  spiral  behavior  of  the  Fresnel  integrals  predicts 
a  decrease  in  the  amplitude  of  I  lit  se  intensilv  oscillations  with 
decreasing  o  and  also  a  decrease  in  their  quasi  period.  The 
amplitude  decrease  is  most  evident  at  lai ge  size  parameters 
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Kit».  :t.  bike  Fig,  2  Ittit  with  (ii)  Iw  =  ltMH)  mul  the  electric  vector 
perpendicular  to  scattering  plane  and  Oil  kti  =  1(M)  and  the  electric 
vector  parallel  to  the  scattering  plane. 

(Fig.  2)  where  the  diffraction  pattern  is  lean  spread  out;  at 
small  size  parameters  the  diffraction  decrease  is  lost  because 
of  the  increasing  amplitude  of  the  intensity  oscillations  as¬ 
sociated  with  the  interference  of  the  p  =  1  ray.  These  oscil¬ 
lations  are  also  of  quasi-period  5(A/n)l/'-!  rad,'2  diminishing 
with  </>•  In  addition  to  the  j  polarizations  of  Figs.  2  4,  the 
physical-optics  approximation  has  been  verified  lor  the  j  not 
shown  and  for  both  j  at  kn  =  5000.  Both  the  approximation 
and  the  Mie  theory  predict  that  the  positions  of  the  coarse 
maxima  depend  only  weakly  on  j;  the  coarse  maxima  lie 
slightly  closer  to  </>,.  when  ;  =  I.  The  approximation  predicts 
that  the  intensity  variations  of  the  course-structure  oscilla¬ 
tions  are  greatest  when  j  -  2,  which  agrees  with  Mie 
theory. 

Mie  theory  shows  that  the  finc-ulructurr  intensity  varia¬ 
tions  are  larger  for  j  =  1  than  for;  =  2;  this  is  probably  because 
p  *  2'  scattering  is  greatest  when  j  =  I  .a  For  lit)0  <  •/>  <  00°, 
thi.  positions  of  the  fine-structure  intensity  oscillations  for 
j  =  2  are  typically  shifted  from  those  for  j  =  I  by  one  half  of 
a  fine -struct ure  quasi-period,  except  when  </>  is  slightly  less 
than  i />,.. 

It  is  evident  from  Fig,  5  that  the  coarse  structure  of  the 
phase  difference  for  Mie  scattering  roughly  follows  that  of  the 
physical-optics  approximation.  This  has  also  been  verified 
at  ka  =  25.  The  gently  sloping  curve  of  the  approximate  re¬ 
sult  near  <1 v  may  he  understood  in  terms  of  the  phase  advances 
b)  of  the  (totally  reflected)  p  =  0  ray;  f>>  is  approximately 
(n„/n,Y2/> i  near  </>,  . I>a  Comparison  of  Figs.  5  and  •!  shows  that 
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the  maxima  in  the  coarse  structure  o!  |A|  are  associated  with 
minima  in  the  coarse  structure  of  the  /,  An  anomaly  is 
present  in  Fig.  5(h)  for  •/>  ~  .("'J.  Mie  computations  of  /,  in 
dicate  that  /|  is  quite  small  at  this  anomaly;  for  ha  =  10(1, /, 
has  a  minimum  value  of  H  X  It)  1  in  the  units  of  Fig.  4(h)  at  </i 
=  Itti. 705°.  Tire  plot  of  5  lor  An  =  25  reveals  a  similar  anomaly 
at  0  ~  2:t°. 

The  range  of  huhhle  sizes  associated  wit  h  Figs.  2  5  will  he 
illustrated  hv  considering  green  illumination  with  a  wave 
length  in  air  of  550  nm  and  a  wavelength  in  water  A  of  4  Id  mu; 
then  An  -  10,000  gi\  -s  ii  =  (552  pm  and  kn  =  25  gives  a  =  1.0 
pm. 

It  is  remarkable  that  the  physical  optics  approximation 
describes  the  general  features  of  the  broad  intensity  maxima 
at  </>  ~  45°  for  kn  =  25.  A  stationary  phase  approximation 
was  used  in  Kefs.  1  and  2  to  eliminate  a  phase  term  |i /,(.'.  a  ) 
in  Ftp  *24)  of  Kef.  !|  from  the  diffraction  integral.  The 
arguments  given  in  Kef.  I  that  justify  that  approximation 
would  appear  to  break  down  for  kn  -  25.  The  rough  success 
of  the  model  indicates  that  the  maximum  can  lie  attributed 
in  part  In  the  first  diffraction  maximum  and  in  pari  to  the 
interference  of  the  refracted  wav  associated  with  the />  =  1 
ray;  however,  25  may  be  close  ,o  the  lowest  value  of  kn  for 
which  this  model  is  useful. 

The  principal  purpose  for  developing  and  lasting  the 
physical-optics  approximation  was  to  see  if  the  coarse 
structure  in  the  scattering  could  be  related  to  the  diffract  ion 
and  interference  of  reflected  and  refracted  waves,  The  model 
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Fig.  I  l.fkc  Fig.  2  hut  with  It  a  =  2'>  and  the  electric  vector  (it)  per 
pendit  ular  mul  (lit  parallel 
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Kin  A  I  'I  lasc  iliHei'enre  Inr  l  hi'  sial  It-riiii:  amplitudes  in  t  lit*  l  %v. . 

split's  Ini-  lit)  l/n  -  1 1  Mill  mill  I  It)  An  =  Inn.  Tliesulid  curve 
l-  Irmti  Mic  theory  mil1  Ki|.  I  1 1.  The  <  l:islu<  I  rune  is  I  lie  plivsii  ;il 
optics  n|>|>r< >x i mui it >n  hv  Kip  r.lHi  ui  |<h. 


mny  Iiiivc  practical  applications,  however.  even  tlnni|i!i  ii  i  anils 
I  hi'  line  structure  in  I  hi-  srnlli'iiiip.  For  «>xmii|>lc.  in  the 
si'iilti'i'inn  by  a  spherical  polydispcrsion  of  Ini lililfs.  the  lim- 
structure  wuulil  l  >»•  lust  anil  I  lit'  course  structure  retained  it  1  lit- 
spread  ui  linlililc  sizes  nr  wavclenplhs  was  mil  inn  larpe. 
l‘'nrl  hcrnitirc.  I  lie  inmlcl  inav  lie  tiscliil  in  I  lie  il"si|;n  nf  in 
st  mincn  Is  |  lint  depend  on  I  lie  seat  Iclin):  of  liphl  by  bubbles. 
Kxamplcs  inrlllilr  lal  devices  Ini  in  ■  iln  nplical  measurement 
nl  iniimlnililili'  pnpulal inns  al  si*a;  ilil  lasta-  Doppler  nnc 
ninnii'try  of  bubbles  in  liipiids;  anil  It)  I  lie  deled  inn  nl  Imlilili' 
inipi'i'li'i'l inns  in  pluss.  .Inlinsmi  anil  ('nnlci,l:l  liavi*  nsi'il  la) 
In  nil!  nil)  evidence  nl'  a  link  nl  I  in  III  lies  al  sea  with  radii  nl  less 
than  approximately  ill)  pm.  in  apparent  disutpeement  with 
similar  acoustic  measiiremenls.  Their  inslrumeiil  made  use 
nl  seallei  ilii'  wil  Ii  «/.  ~  !)()“,  wliieli  is  slmwn  by  Fi|;s.  H  and  I  In 
lie  miieli  smaller  Inr  mierulnilililes  I  hall  Inr  a  penmelrir  re 
Heel  inn  from  a  perleel  rcflcdnr.  Hnlli  the  pliysieal  nplies 
approximation  and  the  Mie  (henry  sujtpcsl  that  hiihhle  ile 
I  eel  iull  insl  in  nielli  s  rel  viiii;  nn  I  list  ill  ie  seal  I  eri  lip  silt  mid  ai 
ii'pt  I  ip  lit  wil  h  less  I  lian  I  hal  nl'  I  he  first  dill  rail  inn  maxi 
nnilli  i's  predieled  by  Kip  KJ).  ( ienniel  l  ie  npl  ies1 "  pives  /,  > 


1  Inr  all  i/i  <  ill,;  ennseipient ly.  it  significant ly  nverestiniales 
I  he  seallcred  intensity  elnse  In  <;>,  when  An  <  101)0. 

Keller1 1,1  has  altempted  In  use  lipllt  seallerinp  with  f/.  ~ 

1 10“'  In  measure  I  lie  si/e  sped  rum  nl  mierulnilililes  in  water  fur 
the  purpnse  nl  lestinp  models  nl  eaviliilinn.  References  I  I 
ami  In  pive  anpiilar  seallerinp  pal  terns  (reported  to  lie  oil 
I  aiiicd  I  min  Mic  theory)  with  A  ~  O.liadH  pm/l  .Til!  and  a  = 

l .  a,  7.a,  and  10  pm  (which  correspond  In  Ini  ~~  III.  (ill,  100,  and 
K!dl.  Comparison  nl  the  results  nl  our  Mie  prupram  and 

m. idel  show  that  Keller’s  patterns  err  significant  ly;  they  do 
mil  include  line  structure  and  they  incorrectly  describe  the 
coarse  si  ruclurc. 

We  are  prateful  In  W.  .1.  Wiscotnbc  Inr  pruvidilip  the  iuitita! 
cninpuler  prop ra ill  from  which  I  lie  proprnm  used  here  was 
derived  anti  In  II.  Medwin  Inr  brinpinp  I  be  work  nl  •Inbnsnii 
and  ( 'nuke  lunur  allenl  inn.  This  research  was  suppurled  in 
part  by  I  lie  Washington  Slate  l  hiiversily  Kcsearcb  and  Arts 
( 'mil mil  tee  and  by  I  lie  ( )|Tice  of  Naval  Kesearch.  I '  I  .  Mar 
ston  is  an  Allred  I'.  Sloan  Research  Fellow. 

‘  A u I  bur  In  wlioni  correspondence  should  lie  addressed. 

RKI'KKI'.NC'.F.S 

I .  IV  I ..  Mar-luii,  "I  't  il  leal  ancle  si  al  lerinu  by  a  linlililc:  pliv-iciil 
tipi  ies  appiex  anal  ivn  and  observations,"  -I.  Dpi.  Sue.  Am.  lit), 

I'.’ii;.  i  .o  1 1  pi', 'lip  7(1,  a.'idiK.i  ipism. 

I’.  I.  Alar-lmi  and  I).  I.  Kinpsliuiy,  'Kcalterinp  In  a  Iniliblc  in 
waler  near  I  lie  via  I  leal  anple:  iiilerlereiice  el  I  eels."  .1. 1  >pt .  Sm 
Am.  71,  I!)'.!  I'.Mi  l  I'.IMl  i. 

I  i  Mie.  "Heilrape  /nr  Oplik  Irnlier  Meilien.  spe/iell  kniloidaler 
Melallusunpen."  Aim.  I  'In-.  1 1, cip/ ip l  25. 1177  I  la  i  I  III  IK  i 
I  1 1  ( van  di-  1 1  nl-l .  I,u;li  I  Sitilh'fini!  In  Sinn  1 1 I'nrtn  h  s  I  Wilev, 
New  York.  11157) 

a.  W.  .1.  Wi-niinhe.  "Improved  Mie  -eallermp  aliioril Inns,"  Appl 
Dpi  III,  I  at  la  I  al  II I  HIIKIII 

ii  II.  II.  Hull.  H  < )  .  (  lilliipreclil .  and  ( '  M,  Sliepc  eviili,  "Tlieorel  leal 
liplil  si  at  term);  eoelln  ieiils  lur  relative  refractive  indexes  less  than 
mill x  and  lur  totally  rellcrlinp  spheres."  .1  dpi  Sue.  Am,  II, 
|M  .0  I  Ida  1 1. 

7.  It.  II.  Hull  ii  nl  .  Tahir*  nl  l.iuhl  Si  nl  Iiiiiii;  /■'uinlmii.-  U' 
Muhipan  I’ress,  Anil  Arbor.  Mil'll.,  lll.aKi 
s.  1,1,  /.elniaiiuvieh  and  K.  S.  Slid  rill,  Tahir*  o/  l.iuhl  Sail  I  mnu. 
Wirl  III  l  iii'llii  ti  nl >  a/  h.iliin  linn  Sml  Irnin:,  mill  l.wht 
l‘rr**urr  1 1 1  vilmmeleorolopinil  I’ldili-limp  House.  I.cninprad, 
I  inis  i. 

P  M  Kerker.  77m  St  nttri'iin:  nj  l.nilil  'iinl  iltlirr  /'.7e,  t rnniiiniirt t> 
Ihntuii inn  I  Academic.  New  York,  IIMilll. 

In  I  VI  lave.  ' 1  Seal  lemu;  ul  visible  liplil  by  larpe  waler  spheres." 
Appl  dpi  h,  i mi  i imiKi. 

II  W.  .1  Wlscoeihe.  "Mie  Srallei  ill);  I 'alvillul  iolls;  Advanees  in 
I'ei  hniipn  and  Fast.  Vector  Speed  Cuinpuler  < 'odes,"  NI'AK 
I 'eel  i  nival  Nute  N<  ’AH  TN  1  III  I  STI(  I  Nal  ional  ( 'e  liter  lur  A I 
niusptieric  Heseareli.  Huuldei  .  ( ’ulu.,  107111, 

I'.1  (i  K  I  lavis,  "Seal  I  ei  iln;  ul  I  in  III  hv  all  air  bubble  m  waler."  .1. 1  >pl 
Sue.  Am.  15.  .v;.;  hi  i  i;i;,:,i 

Kl  H  I ).  duluisun  and  l\.  t  ’  t  'nuke.  "Hubble  pupil  lal  inns  and  speel  ra 
ill  ruaslal  waler-:  a  plinluprnp'.ic  approaell,"  .1,  I  ieuphvs.  lies, 
St,  :17C,I  :(7lili  1 11)7111. 

II.  A  Keller.  “Kill  Slreuliehl  /.aldverlalueh,  auprwandl  rnr 
Hesl i ii i Him i);  des  Kavilatiuiiskeiinspeklriinis."  dpiik  .'12,  Ilia  170 
1 10701. 

I  a.  A .  Keller,  "The  inHuei|ve  ul  I  lie  ea  vital  inn  in  irleus  speel  rum  oil 
cavil  a!  inn  iueept  i,  in.  invest  ipated  with  a  seal  (eled  liplil  i  unlit  ill); 
Inel hud,"  .1,  Kasir  Klip  III,  1*1.  ‘I'.’a  I  I!I7'..'I. 


SUPPLEMENT  TO  PAPER  NO.  4 


This  supplement  compares  unpublished  Mie  computations  (solid  curve) 
with  the  physical -optics  model  from  Paper  3  (dashed  curve).  The  refrac¬ 
tive  index  ratio  is  n./n  -  0.75.  The  notation  is  as  in  Papers  3  and  4 

i  o 

except  In  Figures  SI  and  S2  we  give  the  total  normalized  scattering 
I  =  (1^  +  1^)  / 2  for  the  case  of  unpolarized  incident  light.  For  Figures 
S3  -  S5,  scattering  for  the  same  ka  but  for  the  orthogonal  polarization 
is  given  in  Paper  No.  4;  the  associated  figure  from  that  paper  is  noted 
in  parenthesis. 

FIGURE  CAPTIONS  FOR  SUPPLEMENT 
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=  25. 

Fig. 

S2 

Log 
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10)  of 

Iu  for  ka 

=  100. 

Fig. 

S3 

F 

for  ka 

=  100  ( 

see  F ig. 

3b). 

Fig. 

S4 

>2 

for  ka 

=  1000 

(see  F ig. 

3a ) . 

Fig. 

S5 

5 

for  ka 

=  10000  (see  Fig 

.  Za ) .  See  a  1  so 

Fig. 

4  of  Paper 

1  for  <j) 

from  75“ 

’  to  85°. 

Fig. 

S6 

Pha 

se  difference 

for  scattering  amplitudes 

(as 

in  Fig.  5a, 5b) 

for  ka  =  25.  Note  that  the  model  does  not  give  the  fine 
structure  present  in  the  Mie  result. 
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One  of  the  causes  of  light  scattering  in  glass  is  the  pres¬ 
ence  of  small  gas  hubbies.  In  critical  applications  such  as 
optical  fibers  this  scattering  is  not  negligible;  a  microbubble 
trapped  at  the  core-cladding  interface  may  scatter  a  signifi¬ 
cant  fraction  of  the  light  in  the  fiber.1  Here  we  use  Mie  theory 
and  physical  optics  to  describe  the  main  features  of  bubble 
scattering  in  glass.  We  anticipate  that  these  results  will  be 
useful  for  the  detection  of  bubbles  in  quality  control  appli¬ 
cations. 

We  are  unaware  of  prior  publication  of  the  Mie  intensity 
for  0  180°  scattering  from  bubbles  in  any  dielectric.  Struc¬ 
ture  in  the  Mie  scattering  from  bubbles  can  be  related  to  dif¬ 
fraction  and  interference  near  the  critical  scattering  angle  and 
to  the  Brewster  null  in  the  reflectivity  of  parallel  polarized 
light.  For  scattering  from  a  sphere  with  a  refractive  index 
which  exceeds  that  of  the  surroundings,  the  critical  scattering 
angle  is  not  present;  also,  the  effect  of  the  Brewster  angle  is 
less  pronounced. 

We  approximate  the  bubble  shape  as  a  perfect  sphere  and 


assume  the  glass  has  a  refractive  index  of  1.46,  as  does  the  glass 
used  in  the  cladding  and  outermost  core  layer  of  a  typical 
optical  fiber  preform.  (In  the  resultant  fiber  a  bubble  trapped 
at  the  core  -cladding  interface  is  broken  up  into  a  sequence  of 
much  smaller  bubbles  elongated  in  the  direction  the  fiber  is 
drawn.  ‘ )  We  take  the  bubble  gas  to  be  air,  so  that  the  relative 
refractive  index  n,/n„  ^  0.685.  A  bubble  is  described  by  its 
size  parameter  ka,  where  k  =*  2ir/A,  A  is  the  wavelength  of  the 
light  in  the  glass,  and  a  is  the  bubble  radius.  Here  we  con¬ 
sider  bubbles  with  size  parameters  5,  25,  and  100.  For  A  = 
0.43  pm  (which  corresponds  to  632.8  nm  in  air  for  visual  in¬ 
spection)  ka  =  100  gives  a  =  6,9  pm;  ka  =  5  gives  a  *  0.35 
pm.  For  A  =  0.56  pm  (0.82  pm  in  air,  a  typical  operating 
condition  of  fibers)  ka  =  100  gives  a  »  8.9  pm. 

If  the  bubble  is  illuminated  by  a  collimated  light  beam,  the 
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Kig.  1.  l.ognrithm  (base  10)  of  the  normalized  scattered  intensity 
predicted  hy  Mie  theory  (solid  curve)  for  ka  =  100  and  the  electric 
field  perpendiridar  to  the  scattering  plane.  The  dashed  curve  is  the 
physical  optics  approximation  given  hy  Kq.  (24)  of  Kef.  It  which  is 
useful  when  scattering  angle <l>  is  close  to  the  critical  scattering  angle 
it,,.  ~  94°.  The  approximation  fails  to  describe  forward  region  ( </■  < 
20°)  and  the  backward  regt  >n  ( >T>  ~  180'’), 


Fig.  2.  hike  Fig.  1  but  with  electric  field  parallel  to  the  scattering 
pintle  (y  =  2  scattering). 
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Fig.  It.  Scattering  for  ka  =  25.  Mie  theory  for  j  =  I  (dolled  curve 
and  7=2  (solid  curvet.  Dashed  curve  is  /  j  from  physicul  optic 
approximation. 
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angular  distribution  of  the  scattered  intensity  is  described  by 
Mie  theory. -  At  a  distance  R  »  ha-  from  the  bubble  center 
and  a  scattering  angle  </>,  the  ./'-polarized  scattered  intensity 
is  given  by  the  incident  /-polarized  intensity  multiplied  by 
/;(</>)(«//f  )'-'/■•.  where  j  =  1  denotes  polarization  of  the  elect  ric 
vector  perpendicular  to  the  scattering  plane,/  =  2  the  parallel 
case.  The  normalization  of  the  functions  //(0)  is  chosen  so 
that  /,([/>)  =  1  represents  perfect  reflection  according  to  geo¬ 
metric  optics/1,4  Using  a  program  described  in  Ref.  5  incor¬ 
porating  a  modified  version  of  the  MIKVo  subroutine  devel¬ 
oped  by  Wiseombe,"  we  obtain  the  //(i/>)  functions  shown  in 
Figs.  1  4. 

The  coarse  structure  in  the  vicinity  of  the  critical  scattering 
angle  </>,  may  be  understood  with  reference  to  a  physical  optics 
approximation  developed  by  Marston. 11,7  Most  of  the  scat¬ 
tered  intensity  in  this  region  is  due  to  three  types  of  rays  which 
we  denote  by  the  parameters  p  =  0,  1 , 2,  where  p  is  the  number 
of  chords  within  the  bubble  and  the  p  =  0  ray  has  only  an  ex¬ 
ternal  reflection.4  When  the  local  angle  of  incidence  0  satisfies 
sinfl  >  sinW,  =  n,/n„  the  p  =  0  ray  is  nearly1'  totally  reflected; 
the  abrupt  boundary  of  total  reflection  at  0,.  results  in  a  far- 
field  diffraction  pattern  of  this  ray.  Rays  with  II  <  II,.  are 
mostly  transmitted  inside  the  bubble.  Rays  with  p  —  1  are 
refracted  twice  and  interfere  with  p  =  0  rays  tn  the  region  </» 
<  i/i, where  </>, .  =  ?r  —  211, ■  is  the  scattering  angle  of  rays  inci¬ 
dent  at  llv.  Mnrston’s  physical  optics  approximation  con¬ 
siders  the  diffraction  and  interference  of  p  =  0  and  1  rays. 
Small  quasi-periodic  fine-structure  variations  ride  on  the  large 
quusi -periodic  coarse-structure  variations  and  are  due  mainly 
to  the  interference  of  certain  neglected  p  =  2  rays  (p  =  2'  rays 
in  the  notation  of  Refs,  II  and  5).  If  the  average  intensity  over 
an  angular  range  containing  several  fine-structure  variations 
is  desired,  the  approximation  may  be  used  to  predict  / ,  in  the 
range  20°  5  </>  <  »/>,.  =  9!i.f>4°.  This  model  breaks  down  for 
small  bubbles,  as  can  be  seen  in  Fig.  4,  due  to  its  use  of  the 
Fresnel  reflection  and  transmission  coefficients  of  planar  in¬ 
terfaces  and  the  stationary  phase  approximation  of  diffraction 
integrals. 

The  physical  optics  approximation  predicts  a  quasi-period 
5  (X/«)I/'J  rad  for  the  coarse  structure,  and  a  fine-structure 
quasi-period  of  0,95(X/a)  at  </>,- ,  decreasing11  with  < />.  These 
features  are  clearly  evident  for  ha  =  25  and  100,  as  is  the 
predicted  broad  decline  in  intensity  for  </>  >  </>,.  - 
arcsin|1.2(X/u  cosW,  )l/2|  (identified  in  Ref.  5  ns  the  approxi¬ 
mate  location  of  the  last  coarse-structure  maximum). 

The  broad  minimum  in  l  >  near  1 10°  corresponds  to  light 
scattered  from  the  Brewster  angle:  parallel  polarized  light 
incident  on  a  plane  glass  air  interface  is  totally  transmitted 
when  II  =  114.41°.  Thus,  when  ha  is  large  enough  to  use  the 
Fresnel  coefficients,  we  expect  a  negligible  contribution  to  /•.> 
at  0  =  180°  -  2(84.4°)  =  1 1 1.2°  from  the  externally  reflected 
(/)  =  0)  ray.  This  </»  can  also  be  written  as  2  nrctun(n„/n, ). 

For  ha  =  5  (Fig.  4),  the  line-structure  quasi-period  is  ap¬ 
proximately  equal  to  the  coarse-structure  quasi-period.  The 
minimum  in  !■<  at  102°  and  the  slight  dip  around  45°  are  the 
last  remnants  of  this  structure.  As  ha  is  reduced  below  5,  we 
find  that  the  scattering  pattern  approached  that  of  a  dipole 
radiator,  as  predicted  by  Rayleigh  scattering  theory:  /t(</>) 
becomes  a  flat  line,  and  /;>(</>)  is  symmetrical  about  a  minimum 
at  90°. 

We  have  compared  these  plots  with  ot  hers  made  for  n,Ai„ 
=  0.75  (air  bubbles  in  water)  for  the  same  values  of  ha  and 
found  them,  as  expected,  quite  similar.  The  most  noticeable 
difference  is  the  slightly  larger  quasi-periods  for  bubbles  in 
glass.  In  both  cases,  the  best  angles  to  look  for  scattering  will 
be  those  angles  where”  */’  "5  <b,  ~  (X/u)l/2.  Corrections  are 
sometimes  needed7  when  computing  the  effect  ive  scattering 


Fig,  4.  bike  Fig.  8  but  with  ha  =  5. 


angle  observed  by  a  detector  due  to  refraction  of  the  scattered 
light  as  it  leaves  the  glass  and  enters  the  detector."  The  coarse 
and  fine  structures  present  when  ha  Z  25  may  be  used  to 
obtain  independent  estimates  of  a.  That  diffraction  is  im¬ 
portant  at  i/»,  can  be  seen  bv  comparing  the  Mie  results  to  the 
geometric  prediction,4  including  only  the  p  =  0  ray,  that  //(</>,  ) 
=  1  and  with  the  physical  optics  prediction  that  //(</>,.)  =  0.25; 
the  Mie  results  lie  much  closer  to  the  latter  than  the 
former. 

The  result  that  scattering  is  best  observed  when  0  5  0,.  - 
(X/a)'7’  rad  can  also  be  used  in  certain  cases  for  scattering  by 
nonspberical  bubbles  if  their  approximate  size  is  known.  For 
example,  if  the  bubble  is  spheroidal  and  either  the  major  or 
the  minor  axis  lies  in  the  scattering  plane,  the  critical  scat¬ 
tering  angle  predicted  by  ray  optics  is  the  same  as  that  pre¬ 
dicted  for  a  sphere.  The  diffraction  reluted  shift  in  the  first 
maximum  is  roughly  (X/ti)l/2  rad,  where  a  becomes  an  average 
(which  depends  on  the  orientation  of  the  spheroid)  of  the 
semimajor  and  semiminor  axes.  The  details  of  the  coarse  and 
fine  structures  will  depend  on  the  bubble’s  shape. 

The  approach  to  the  scattering  described  here  is  to  allow 
the  light  to  enter  and  leave  the  glass  via  plane  surfaces  and  to 
neglect  surface  reflections."  Scattering  near  90°  has  been 
used  to  detect  bubbles  in  glass.111  For  glasses  with  «„//»,  =  1.5, 
0,  =  90.4°  and  the  Brewster  scattering  angle  is  112.6°.  Pre¬ 
vious  applications  of  Mie  theory  to  bubbles  are  reviewed  in 
Refs.  5  and  I  I  Approximations  for  scattering  by  inhomo¬ 
geneities  within  fiber  waveguides  have  recently  been  pub¬ 
lished, 12,1 11  but  these  do  not  describe  the  angular  structure  for 
scatterers  with  ha  >  I . 

'Phis  work  was  supported  by  t  be  Office  ol  Naval  Research. 
1’.  1,.  Marston  is  an  Alfred  1*.  Sloan  Research  Fellow. 
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H.  There  may  he  an  exception  to  this  guideline  for  large  spherical 
hohhles  heeause  backseat  tering  can  he  enhanced  hv  the  axial 
focusing  of  p  =  It  glory  rays.  Observations  of  this  enhancement 
for  hohhles  in  liquids  are  described  in  I*.  1..  Marston  and  1).  S, 
l.angley,  .1.  Opt.  Soc.  Am.  70,  l (>(17  ( I9H0I. 

9.  Kefruetion  corrections,  such  as  those  described 'in  Kef.  7,  are 
simplified  when  the  light  enters  and  leaves  the  scattering  volume 
via  plane  interfaces,  and  a  lens  is  used  to  place  the  detector  in  the 
far  field.  This  may  he  achieved  by  immersing  the  glass  in  an 
index-matching  liquid  and  illuminating  and  observing  it  via 
windows. 
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Description  of  computer  codes  for  Mie  and  model  computations 
(D.  L.  Kingsbury  and  P.  L.  Marston)  [This  is  a  slighlty 
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in  Water  and  Glass,"  M.  S.  Thesis,  Washington  State  Uni¬ 
versity  (1981 ).] 
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Append i x 
COMPUTER  CODES 


Three  programs  were  used  to  generate  our  plots.  MIESC, 
incorporating  a  modified  version  of  Wiscotnbe's  MIEVO,  produced 
the  Mie  and  CHIT  the  approximation  results.  The  versions  given 
here  store  results  in  disk  data  files  where  they  can  be  accessed 
by  the  plotting  routin  MTEPL  for  display  on  the  HP  2847A  monitor 
or  for  paper  plots  on  the  HP  9872A  plotter.  All  programs  are 
written  in  Hewlett-Packard  RT-IVB  Fortran.  CHIT  evaluates  the 
approximation  developed  by  Mnrston  and  Kingsbury  (1981). 

MIESC 

MIESC  asks  for  the  following  inputs: 

1.  REF  INDX  is  the  complex  relative  refractive  index  n. /n0 . 
The  imaginary  part  must  be  non-positive. 

2.  BETA  is  the  size  parameter  ka  ~  2rra/A,  where  A.  is  the 
wavelength  in  the  outer  medium. 

3.  AST  is  the  smallest  scattering  angle  for  which  scattering 
angle  for  which  scattering  functions  are  computed,  DPHI 
the  angle  step  size,  and  NUMAMG  -  Nan  is  the  total  num¬ 
ber  of  angles. 

If  NUMANG  is  even,  results  are  calculated  for  the  angles 

j  ~  AST  +  i(DPHI),  i  =  0  to  NUMANG /2 ,  and  their  complementary 
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angles  (180°-  <pt  )  .  (If  NUMANG  =  1,  90°  must  be  the  only  angle. 

If  NUMANG  =  0,  only  the  extinction  efficiency  factor  Q Jt>  is 
computed . ) 

For  example,  if  n,  =  1.33  -  i(.08),  n  =  1.0,  and  scat¬ 
tering  functions  are  desired  for  4*  ~  ‘10*  to  90°  with  increments 
of  .5°,  one  should  type  in 

/. 7508, -.0282, 40. , .5,101 

MIESC  then  asks  for  the  name  and  cartridge  number  of  an  existing 
datafile  in  which  to  store  the  output.  What  data  is  stored 
depends  on  what  the  output  buffer  RBUE'  lias  been  equivalenced  to. 
In  the  version  included  here,  each  record  v/ill  consist  of  three 
real  numbers  giving  4',  I,  (<£) ,  and  I^^),  and  the  records  are 
entered  in  order  of  increasing  <t> ,  except  that  the  record  for 
(180*-  $ )  immediately  fe  .lows  that  for'  4*  ■ 

If  the  phase  difference  DIFA  is  not  desired,  the  code 
that  computes  FH1,  PH2 ,  and  DIFA  is  skipped. 

To  bring  subroutine  M1EV0  down  to  more  manageable  size, 
the  original  comments  were  largely  omitted.  The  commented  ver¬ 
sion  is  given  in  Wiscombe  (1979).  Most  complex  variables  have 
been  converted  to  2 -dimensional  arrays  or  to  separate  variables 
prefixed  with  "R"  or  "I."  For  example,  NUM  became  NUM(l)  ana 
NUM(2),  denoting  the  real  and  imaginary  parts,  respectively,  and 
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ZET  became  RZET  and  IZET.  (This  was  done  because  of  the  absence 
of  double  precision  Fortran  in  RT-1VB  Fortran.)  Most  complex 
arrays  were  similary  converted  to  separate  arrays  for  the  real 
and  imaginary  components. 

Subroutine  names  prefixed  with  "DV"  refer  to  functions 
in  the  HP  Vector  Instruction  Set,  which  considerably  speeded  up 
the  execution  of  two  former  DO- loops. 

To  keep  the  loaded  program  within  a  32K  size  limit,  the 
original  MIEVO's  complex  BIGA  array  was  discarded.  (BIGA,  which 
stored  the  A„ ,  contained  as  many  elements  as  the  number  of 
terms  used  in  the  Mie  series.  Thus  for  ka  =  10  000,  a  total  of 
2(10  088)  =  20  176  double  precision  must  be  stored,  requiring 
161  408  bytes.)  Required  An  values  are  instead  stored  on  and 
retrieved  from  a  system  disk  using  EXEC  calls. 


The  use  of  the  Hewlett-Packard  minicomputer  instead  of  a 
regular  "pay"  computer  system  was  absolutely  necessary.  Extra¬ 
polation  from  Table  8  of  Wiscombe  (1980)  indicates  a  MIEV0  exe¬ 
cution  time  of  over  100  seconds  on  a  CRAY-1  for  ka  =  10  000  and 
NUMANG  =  8000.  The  several  size  parameters,  repeated  testing, 
and  dozens  of  plots  required  made  our  choice  of  computer  systems 
inevitable. 
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CRIT 

Required  CRIT  inputs  are  the  relative  refractive  index 
RI  (now  assumed  real);  the  size  parameter  BETA;  PTS  and  STRT, 
which  determine  the  angles  for  scattering  computations;  and  LAST 
ANGLE,  the  largest  PH  (c p)  .  N  =  n0/ni  ~  1/RI  is  used  to  calcu¬ 
late  the  critical  incident  angle  : 

THC  =  arcsi n ( 1/N ) . 

(Arctangents  are  used  in  the  program  due  to  the  unavailability 
of  the  arcsin  function.)  The  critical  scattering  angle  <f>  is 

then 


PHC  =  Tr  -  2  ( THC  )  . 

For  a  given  incident  angle  TH-THC,  the  angle  PH  at  which 
the  refracted  ray  is  scattered  is  given  by  an  equation  derived 
using  Snell's  Law  and  elementary  geometry 

IH  -  2(arcsin[ sin(TH)/Nj  -  TH). 

For  a  desired  smallest  PH,  one  chooses  STRT  such  that  TH 
=  STRT(THC)  substituted  in  the  above  equation  gives  this  PH. 

As  CRIT  executes,  TH  is  stepped  up  in  increments  of 
E  =  (THC  -  THiltit  ) /PTS .  As  TH  approaches  THC  the  corresponding 


step  sizes  for  PH  grow  steadily.  Consequently,  at  selected 
angles  TH2 ,  TH3 ,  and  TH4 ,  the  step  size  in  TH  is  reduced. 


The  function  G1  is  the  geometrical  divergence  factor 


sm(TH)  cos(TH) 

G  I  rr  — i - - - - - - 

2 1  1  'L^co»(TH)/cos(R)]|si,,fPH) 

in  the  expression  for  the  intensity  of  refracted  radiation 
(Davis,  19S5;  van  de  Hulst,  1957),  where  R  --  arcsin(Nsin(TH)  )  . 

IP  and  IS  are  the  intensities  of  the  parallel  and  per¬ 
pendicular  refracted  radiation,  respectively,  computed  using 
Fresnel's  intensity  transmission  coefficients  (1  -  FP )  and 
(1  -  FS)  for  plane  surfaces.  The  coefficients  are  squared 
because  this  radiation  is  refracted  twice  --  at  entrance  and  at 
exit . 


The  variable  W  is  the  upper  limit  of  integration  on 
the  Fresnel  integral  used  in  the  reflected  radiation  calculation 
(Marston,  1979).  Subroutine  FRES  evaluates  this  integral  using 
a  numerical  approximation  (Abramow.itz  and  Stegun,  1955). 

The  phase  differences  B1  and  B2  between  the  reflected 
and  refracted  rays  for  perpendicular  (Bl)  and  paral  lel  ( B2 ) 
polarization  have  three  components,:  the  phase  shift  the 
reflected  ray  undergoes  at  the  bubble  surface  (DELI  and  DEL2 ) , 
the  phase  path  difference  E’J’AC,  and  a  Tr/4  phase  difference  aris¬ 
ing  from  the  curvature  of  the  scattered  wavefront.  (For  the 
reflected  ray  this  phase  shift  is  already  included  in  the  Fres¬ 
nel  integral . ) 
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CRIT  is  now  ready  to  compute  the  scattering  functions  at 
angle  PH: 


s,  -  VTs  +  (l/,{2){  FC  +  i F3 ) exp ( - iBl ) 
S»  =  ffp  +  ( 1/Y"2  )  (  FC  +  iFS )  exp  ( - iB2 ) 

where  FC  and  FS  are  the  output  of  subroutine  FRES. 


M1EPL 

Output  was  plotted  by  program  V.TEPL  using  the  HP  1000 
Graphics  package.  MIEPL  is  not  included  in  the  program  list¬ 
ing  which  follows. 
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IF<  1ST  .  GO  .  I  H y  .»  ChLL  C  LOSE',  I  DCS  > 

WRITE'-  LU  IN.  2  007  »  I  ERR 
FORMAT*  "  I  ERR*" , 15  > 

USE  THIS  FORMaT  FOR  MORE  DETAILED  OUTPUT 
FORMAT':  1H5 ,5;-;,  HIE  SITE  PARAMETER  =  FS  .  2  ,  1  5:: REFRACTIVE  , 
k  '  INDEX  ='  .D14.6,E12.3//  ANGLE  Mix,  S-SUb-1  ',21a,  'S-SU6-2', 
♦  1  5:-' .  II ITEMS  I  TV  '  .  2X  ,  OEGPOL  .  S,  .'Ill  2/  12'  .» 

FGRMhT  <  F7 . 2 . 5E  i  4 , 6 ,  FS  .4.204.6,  E6 . 3  ,< 

FORMATS  F7 , 2, F 1 2 , 3  > 

END 


SUBROUTINE  MIEVOc  IU  > 

COMPLEX  SBACK ,  St,  S2 
DOUBLE  PRECISION  BET  A , NMU , P I 
REAL  H2CUT 

COMMON/ 1 NOUT/P I C 2  •  . SBACK . SIC  512  v,  S2c  51 2 >, BETA, XMUC 512 > , N 2 C U T , 

*  HUM A NO , NMP 1 , :  0 ! . DENT , lUOUT , KNUD. I LEM, I SECT, MR 0, IOP1 0, IOF20 
LOGICAL  NO  I  MAG ,  NOuNGS 

COMPLEX  ANN i . BNM 1 

D I  MENS  I  ON  RSPC  256  > ,  I  SP*  256  > ,  PSt'K  256  > ,  I  SM>  256  .<  .  PSPSC  256  .> ,  I SPSC  25' 
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*  2 1  N  VC  2  >  .  FFC  2  > ,  AC'.  2  >  .  DECK  2  .) ,  HUM';  2  >  ,  EPS  1  .  EPS2  ,  TEMC  2  )  ,  TNAC  2  >  , 
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*  S  1  '■  257  > ,  P  I  NM  1  1  >  > ,  ■'  S2C  25  7  >  ,  TMP'  1  >  ’ 

DATA  EPS  1  /  1  ,  0-2/ ,  EPS2/1  .  0-9/  .  MAMIT.-  10000/ 


IFCBETA.LT,  0.  .OR.  PIC  I  >  ,  le  .  0  ,  .OR.  RI<  2  > ,  GT  .  0  ,  :>  STOP  10  00 
HOANGS  =  HUM Am G . EO . " 

NO  I  Mac,  =  PI'  2  > ,  CiE  ,  H2C.L1T 


CALCULATE  NUMBER  OF  TERMS  IN  MIE  SERIES  •  A  LEAST  UPPER  BOUND 
USING  EMPIRICAL  FORMULAE  FITTED  FTP  SIZE  PwFAMETEPS  UP  TO 
2  0, 0  0 0 


IF-:.  BETA.  LT.  4200.  >  NT  =  BETAM  05+BFTa  *  1 
IFC  BETA  .  GE  ,  4200  .  >  NT  =  BE  Th  M  » BET A  !  ♦  c  1  ,  , 
NTP 1  =  NTH 


COMPUTE  BIGh 
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r'  •- 

-OF  oh  ;  .'stem  rise  FOR  aRRA  6  c 

C  m  L  L 

C  ■  4  . 

i  TRw!  I':IH  ,  I L,  ISL  ,  I  :>E'  T  :• 

ICNWD 

- 

IDK 

r 

ILFN 

2 

TfiL  1 

l  '  • 

--  \ 

.  1  »•  1 

-  t 

I  or  2 

-  i 

’■  r  c 1 

ChLL 

ix  o 

i  .i , 

1  L‘  o  0  '  ,  ki  ,  T£M  1 

67 


2IHV-:  1  .)  =  TEM-;  I  VBETm 
2INVC2.')  =  TEM*.  2  VBETm 

PREPARE  FOR  DOWN -REC URREMC  G - 

COMPUTE  INITIAL  H I GH-ORDER  BICA-M  USING  LENTZ  METHOD 

FF<  I  >  =  NTPt  -ZINV< )  > 

F F <  2  >  =  HTP1  ■*ZINV1'  2  > 

MM  -  -1 
KK  =  2* NT +3 

AK*.  1  '  =  <  MM*KK  .tttZIHV-;  1  .> 

AK<  2  >  =  <  NM-rKK  :>*Z1 NV*:  2  .< 

DEN1'  1  >  =  aK<  1  > 

DEN*  2  .1  =  AK*  2  ) 

CALL  DC D V *'  0 ,  10  0,0,  ,  FT,  TEN  * 

NUM<  1  >  =  DEN<  1  ')  +  TEN,  I  > 

NUM<  2  >  =  DEN1'.  2  >  +  T  t’i'k  2  ’> 

KOUNT  =  I 


KuUNT  =  KOUNT  +  1 

IF< KOUNT , GT , MAX  1 T >  GO  TO  40 

CALL  DCDV<  -  1  ,  HUM<  1  > ,  NUM<  2  * ,  hII  ,  TEM  > 

CALL  DCDV< - !  , DE1K 1  ) . DEN-  2  > , *<K , TNw > 

IF  >;  TEM1'  1  >  ,  GT  ,  EPS  1  ,  hMD  .  THA*  1  i  .  GT  .  EPS  1  >  GO  TO  30 

ILL-CONDITIONED  CASE- -STRIDE  TWO  TERMS  INSTEAD  OF  ONE 
MM  =  -MM 
KK  =  KK+2 

AK<  1  >  =  MM *KK  )*2IMV>.  I  > 

AK*  2  >  =  <  MM*KK  > '* 2  I  NV*  2  > 

NTHu  1  >  =  h1'1  1  >  *; N U M ' ’  1  t  -  mK<  2  n"NUM*:  2*  +  IDO 

NTNC  2  -  -  Ak*  I  >*NUM<.  2  >  +  nt  *  2  >.*NUM'.  1  > 

DTDC  1  >  =  Ml  1  >*DEN<  1  >  -  hK<  2  WDEIK  2  >  ♦  IDO 

DTO<  2  I  =  A>  .  1  )  *0EiK  2  >  +  Ak*  2  >*REf«  !  > 

CALL  DCDV'I  U  .  NTN'  I  >  ,  NTM'  2  > ,  DTD  ,  TEM  '> 

TMA<  1  >  -  TEM'  1  <  >  I-  F •;  1  .  -  TEN'.  2  >  >  FF<  2  > 

FF<  2  >  =  f  EM';  1  » f.  •  >  +  TEM' 2  >*t  F<  1  > 

FF<  1  >  =  TMA*  1  i 
MM  =  -MM 
K  K  =  k  K  >  2 

Ak'<  1  )  -  \  MM  ft  V  '■  '  Z  I  N  V '  t  .> 
hK  *  2  !  -  '  Mil-fKk  1  KZIH7*  2'  ■ 

CAL.L  DODV':  0  MI  Ml-  I  ‘  I II.  IN'  2  .» ,  NTN  ,  TEN  .> 

HUM*'  1  ■'  =  Ml  '  1  i  «•  TEM'  1  ) 

NUN''  2  •  =  i-O  >  2  >  *■  TEM'  2  • 

CALL  0 C D V ■;  0  ,  DEN’,  )  >,  DEN*  2  ) .  DTD,  TEM  * 

C'l  N-  t  »  =  hi.'.  I  ■  +  TEM  t  * 

Dl  M*  2  ’  =  Ht  2  >  ♦  r £M'  2  * 

KOUNT  =  I  OUMTf I 
GO  T0  20 

CALL  DC  07'.'  0  MUM'  t  '  .  ill  IM'  2  ■  .  DEN  TT  i 

TEM*'.  1.-  =  TT’  1  i * F F \  t  ..  -  TT*  2  >«*FK-.2  > 

FF’  2  •'  =  TT  *  !  **FF-  2  *  +  TT'  2;*FD  1  > 

FF<  i  ■  -  7i;.:  t  : 

'  ".r  C I '  T  0-'  CONVt'F  '  •  «C  t 

I  F'.  DhCV'  T  ;  '-1 D  i.i  'Ll.  EF  22  ,  hND  .  DhD  D  T  ■  '  2  >  •  l.T  EPS  2  ')  GO  TO  SO 
MM  =  -MM 
y  K  =  I  t  •*  2 

•:  I  »  -  1  l'?H  2-1  I  *  2  IN’-'  I  > 

A V 1  2  *  -  '  MM  i  j  »;  >  *  -  i  N ..  •  2  ' 
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C  hLL  C> C C> V 0  1  [>  0  ,  0  .  ,  HUM  ,  TEM  ) 

HUM'  1  >  =  nt  •  1  >  ►  TEM'  1 

HUM1'  2  i  =  hi  <  2  +  TEN1!  I  1 

cull  c>c c> o ,  t  c> o ,  o , ,  den  ,  tem  :> 

DEN>'.  1  :>  =  nt  1  >  +  TEI'K  1  .> 

DEN<  2  >  =  HT<  2  >  +  TEMT  2  > 

GO  TO  20 

40  UPITE'  LLIOUT  .  :?.Ori|  :>  HT,  BET  A ,  P I  ,  AK  ,  NUN  ,  DEN ,  TT  ,  FF 

jOI  FORMhK/,V  com  1HUED  FRACTION  FOR.  w-GUB-HT  FhILEO  TO  CONVERGE 

*  NTs  ,  hi/  '  '  =  '  , E2  0 .  8/  REF  ft  INDEX*  '  ,  2E20 , 3/  wK  =  ' ,  2E20  .  8/ 

*  HUM-  ,  T.E20  '  DEH-  ,  2E2  0  .  8/  '  TT='.2E20  8/  FF  =  '  ,  2E2  0 . 8  ) 

STOP  1 002 

5  0  RBI  Gw  =  FFv  1  .> 

IB  I  GO  =  FFl  2 
BUFF-:  31  :>  =  FBI  Gw 
BUFF'  32.)  =  IBIGh 
MB  31 

DOWNWARD  RECURRENCE  FOR  BIGW-H 

DO  70  H  *  NT ,2,-1 
TEM<  I  :>  =  M  *  2 1  HV‘-  1  ’  +  RBI  Gw 
i  Ef'h  2  >  =  M  ♦  2  I  HV-  2  >  *  IBIGH 
CALL  BCD  Vi  0.  1D0, 0.  ,  TETLTMA  > 

RBI  Gw  -  H+ZIHWt  1  >  -  T  rl  h  •  1  > 

IBIGH  =  toZIriO'.'2>  -  TMW'  2> 

MB  =  TIB -2 
BUFF':  MB  >  =  PPIGw 
BUFFLMB+1  >  =  IBIGh 
IF  '"  MB  .  HE  .  1  ■  GO  TO  7  0 

ChLL  ELEC'  2  I C HOD , BUFF . J LEN , 1 OP  1 , 1 0P2  > 

IF  ■  rl  E!0 . 2  .'  GO  TO  8  0 
MB  - 


1  Til  '  =  l  OF'  1  -  1 
I  OP 2  I  SECT -2 

7  0  CONTI  liUC. 

ChLL.  ELF!  Ci  2,  IrH,;C  .BUFF.  ILL'  I  OF!  ,  I  OF  2  > 
80  MB"  *  MB 

I  OP  10  •=  IOF't 

I  np'  i  =  I  OP.: 


:ui.'  cum  thui: 

Ml:  -  MB  ij 
I  DPI  -  TOPIC 

I  OP 2  s  l  OP 2  ci 

C  hLL  EXEC  '1,1 12  MI.1C* .  LUFF  .  I  LEN  ,  I  OP  1  0  ,  IOP2  0  < 

RIOPIV  -  ID-i. 'PI'  l  ' 

inn  iwl.  T  C  i "  >hHT  I  Hue  USED  FOR  EFFICIi'.UT  OhLlULwT ION  OF 
HUM!- P  I C  mL  CO  i  2FIC  I  ENTS  IH  Hit  OGP  f  EG 
ru  =  U'ii 
PH  -  ID" 
i  1 M  ~  1 

INI  TIuL  t2ii  RKmITI-BCSBEL  FUNCTION  CETh  -lip  UPHWRD  RECUFPEMC 
PC  IH  =  P  ;  l-LTm  . 

CHIN  =  D'-'V  •  -  FT.  > 


r-r.INF-1  PS1U.-BF.Th  -  CHIU 
CHI  HP  !  =  CHID. "6).-  In  f  PS  ID 
P  Z  E  T  r  i  =  P  S  1 1  ) 

IZCTD  =  CHID 
PZTNP!  ---  PS  I  up  1 
IZTHPt  =  >.  H  I  Mi-  1 

IMTTIal  ITS  PPE'v' IOuS  CuE'FIC IEDTS  .  m 

F OR  USE  It)  mS  .'MMETF-'.  FhCI  OP  SERIES 

hNI-I  !  =  >0.0  t.t .  !'• .» 

L<DH  1  =  >'  0  .  0  ,  U  .  0  1 
OEXT  =  0 . 

IIIITIOI.T'E  Hiir,Ul.>;P  F- '  H  PJH  hUD  SUMS 


C  ALL 

I.-7M0V' 

o . 

r»  0 . 

n  . 

r"'  P  ,  1 

.  MM  ) 

CALL. 

DVAOV'. 

0  , 

[-•ii . 

\\  , 

I  ;.P  ,  ! 

,HN  » 

CALL. 

07m  O1"1' 

0  , 

,  D  0  . 

h  , 

ASM,  1 

,  Mfl  > 

CALL 

c>Vf!.-;.v> 

.  [•■  0 , 

0 , 

!  ■;  a  .  i 

,  nn  > 

Call 

i'i'/MOV* 

n 

['  0 

M  , 

.  I;"  ;.l'  S  , 

)  mu 

Cai.L 

PVM07. 

i*t 

.  D'l 

ij  , 

,  i  -  *  '  . 

)  MM 

CALL 

CiVTIC'.V.: 

!  M 

DO 

'!  . 

P'-IV 

1  MM 

CALL 

DVKOV-: 

.  DO  . 

U  J 

,  i 

j  MM 

CALL 

DVtTOV- 

,  D  O  . 

,  o  , 

,  r- 1  MM! 

,  !  fi 

CALL 

07  MOV' 

1 

.  DO  . 

J.l  , 

,  P  I  1  !  ,  1 

.  MM  • 

m 


WRITE1'  LUOUT  ..2  0'  S  -  DIMOUDT.  1 1.  ,  tmt-lntlS 

forma  r>  i .  e  i  .2 . 3  > 


DO  5  00  )i  -  :  t)T 
RE'.  I  Oh  -  DCPP  MS 
IBICH  =  ru.I'F-  Mb'  I  > 

COMPUTE  T'-'E  VARIOUS  tl  WEPIChL  COFI'F  to  IFU 

run  r  ;  -, ♦  -.  j 

i  IJOI  (•••  l  -  Uti  ■  tip  t 
;  tlf.;  -  ',[■!.  1 

I.  “ r  =  i':  !  I  t  ;  r  ; 

HP  1  DM  =  *  P -  * U 1 1 


CwLCUI.hTE  THE  Mil-  SERIES  C  Oi£FF  I C  I  CUTS  L  I  T  T 


IF'.  MO  IM 

•  GO 

TO 

7  i 

1 1 0 

OF  HE 

P  M  L  L  H 

:.E 

call  do 

■  i  i  . 

r.  i  o 

H 

I  L:-  r 

: j  >-t  .  ,; 

I  Tf". 

TEA'  1  ' 

TO".' 

1  > 

P 

TM.'* 

MET  m 

TriAL  i  :> 

=: 

TEA 

1 

I.. 

... 

71W‘ 

1  - 

TEC'  2 

>  1~ 

n-iM"  C  > 

~ 

TE"*' 

1  >  i 

I 

TMP 

1  f 

TEA'  2 

‘•PC 

h  . !  L  !  '  = 

>. 

'  TEA 

’  1  • 

* 

r . 

:  j  m 

r 1  •-.* 

.IP 

* 

Ah 

' 

T  Ah>  i  ■ 

»■  t 

► 

i  !■* 

M*  ! 

i  ij 

Ml)1.  2  1  - 

• 

.  TEA 

V  ‘ 

* 

■  IM 

~  \  >  ¥ 

Mil  1 

) 

- 

TH 

1 W 1  i 

t  t' 

t 

v-i  •  J;  • 

MJ 

TEA'  - 

= 

Pin 

!  \  ?’ 

r 

- 

Gw 

f-  ! 
i  1 

■  ,J  :  t-  ? 

Tt 

i  L.M 

TEA'  2  i 

= 

Pin 

r¥  ; 

\ 

!  j  H 

f  P  i 

1  ^  ‘  ' 

i  •• 

!  C.H 

TiTm1  i  ' 

- 

t"  M  f'.  • 

!  ■  : 

!-. 

... 

T : :  »* 

t  - 

I  L ’  . 

-  I  E 

Trtu'  J  . 

- 

T r  • 

1  •  » 

1 

‘ 

1  f 

Tl"' 

.  F  ” 

r  i 


i  r  i  1  —  i  '■  i  ■  -  •  i  U  u.  i  '  ' 


j  i  P  [ 


EE  Th 


1,  B-S'JE'-N- )  .> 

.  HT  hL.L  hHGI.  E 


NEEDED 


lE-h  hHD  LITTlC-B 

i  ; 
t: 

PS  I  DP  t  .-  HMw  2  •  ' 
!  i r  .  -P'S  I  t)  '■♦'TI'Ih'  2  1  > 


( I:  '  . 


h 

i  F”‘  1  Hr  1  '-tTMp'  2 


3  0  0  CONTINUE 

NO-AB 0 R P T  I  0 H  C h s C 

TEM.  1  .)  =  ■  PI  OR  I  /’♦■PPIuhM'II.-'P.ETh  m-PSTIIP  1 -P.TTN 
TEN''  2  >  =  '•  P I OR1V*  P':]  UnfFM/BL  !  A  >*I2TNr-  1  - IZETN 
TMA<  I  i  =  %  R  l  OP  IV  *  P.B  I  GA-*  FN/6ETA  **PS  I  HP  1  -RSI  H 
ChLL  BCD'  \  0  .  rt  1»  :■  t  .  i.i  ,  ,  rEH  hH  • 

TEH'  1  >  a  •.  f  I  1  '  i  RL  I  Gn +  F  N, 'BE A  I  +  R2TNP  1 -F-ZETN 
TEH1'  2  >  *  PI',  1  T  *R.BI  Gh+FH,'BETA  )+ I  ZTNP  1  -  I  ZE  TH 
TMTK  1  >  ~  ■  P  T  >  1  •  rPB  I  On  *-FM/BETA  )  ♦PSI  HP  1  -  PS  I N 
CALL  DCD’A  i.i  .  Trln*  1  > , 0 , , TEH , BN < 

350  QErfT  a  OE.,r  t-TUi.irlp  l  •-Hi'.L  n»K  1  >  '♦SNGL'  Bit  t  >  >  » 

IF<  HOANGS  >  GOTO  4*5  0 

PUT  HIE  SEP  IE'-;  COEFFICIENTS  IN  FORM  NEEDED  FOR  COMPUTING  S+, 

RmNP  =  COEFF  ¥<■  wN'  1  't-PN*  1  >» 

IANP  =  coerr  Ht!'  2  »♦•£»•!•.  2  >’> 

RC  HP  a  COE  FT  +*  MIC  1  >~PIK  t  >  < 

I3NP  a  C  OEFF  «•>  hM<  2  •  •DIK  2  >  > 

RANPM  -  MM^RhNP 
IaNPM  =  MM* I hNP 
RBNPM  a  NH*fl'llP 
I  BN  PIT  a  MMflPHP 

HDD  UP  SUMS  WHILE  UPWARD  PEC  UP':'- 1 NG  mN'JLmP  FUNCTIONS.  LITTLE  PI 
AND  LITTlE  ThU 

CALL  DVHPV'  ;;NU,  1  .PIN  1  .'/TEH,  1  .  NN  > 

CALL  DVSUB'  V  TEN  ,  I  .  P  I  Mil  1  ,  t  ,  THP  ,  1  .  NN  > 

CALL  DV:  MV>  FM,  THP,  1  ,  I.'.UN,  1  NN  • 

CALL  DVSUB'.  Till  IN  1  ,  P  1 1 . 1 1 1  .  t  TaUH,  1  ,NM  > 

call  dvho /•  l  i m ,  i  . pi mm i  ,  i , hm  « 

ChI.L  DVP  I  V-  M  I  T'N  ,  Ti-|p  1  VTEH,  1  PIN  I  .  NN  ' 

ChLL  1 1 ., .  iDC"  T  INN  l  ,  1  ,  ThUII,  !  .  /ILU,  1  ,  MM 

CALI.  DVP  IV  iMP  .  V  TLU  t  P  ;■  \  ,  p ,  I  Nil  ■ 

CALL  DVP I"1  ! hMP , VTEM , !  1 : P , 1 , I SP , I .MU' 

C  ALL  DVP  I  V'.  PI  M'-H ,  V  T EH  ,  I  ,  F  SMC. .  1  ,  BSMS  ,  I  .  NN  • 

ChLL  DVP  IV'  !  I.'.UT  I! ,  VTLH  ,  1  ,  I 'SI  IS  ,  1  ,  I  SMS  ,  I  ,  NN  ' 

CALL  D  VS  NR'  P  1 1  ill  t  I  .  T  m  IJ  M  ,  I  ,  VTEM,  1  ,  NN  » 

call  dvt- i V'  P'Pmp ,  v n.H ,  i  ,r-.h,  t  ,pc.m,  i  ,  nii  > 

Cull  UVPtV'  !  PI  IP  VIP  1-1,1  t:U.I  i  CM  ,  I  NN  * 

ChLL  D7PIV-,P».-N!--M  V  1 1  T-l  .  t  .  PSP--. ,  '  ASP-.  .  '  ,  (IN  • 

call  U-..-P r,1'  i,.NPM.vren  t  i .  ,  i  isps,  i  . mm  - 

UPI-hTL  RCLEVhIM  OUmNT  l  TIES  Pi.R  HUNT  PaSS  1HFOUGH  LOOP 
45 U  TIM  -=  -  MM 

FN  ---  FMPI 
RN  P.'l it-  I 

A lll-l  I  -  '  MCI  ' ■ . M ' •  l  mI!'  I  ■  '  :  HGL.'  hII>  J  >  '  > 

BUM  I  a  i  IT'!.  Id,!  ■  pH  !  -i  PNGI  '  PM',  ,?  >  >  ' 

i’  hIT  1  M  h  1 1  Til  ' Ei.  Fi.itli  T  ION '•  BY  l.if  i  '■  r.  r  RECURRENCE 

PEL  T  :  '.  TUijlir  !.  I'-E’w  ri.f-i -RELrn 
IZrT  a  •  T  i.iOl  :r  l i  E  r  h  iMZTHP! -  IZETN 
pzr  ni  «  p  z  T1  '■  t 

IZEIN  -  IZT'lPt 
R.  rur-  i  i.  zt-I 
I ;;  r  r  ,  i  •  : .  r 
PS  I  II  =  I  '  I  Hi  ! 
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i 


PS  IMP  1  -  PZ7NP1 
MB  =  MB  *2 

IF  MB  ,  LT  .  33  >  GO  TO  500 
MB  =  1 

IOP2  =  I 0P2+2 

IF  < I0P2.LT. ISECT)  GO  TO  490 
I OP  I  =  IOP1+1 
I OPS  *  0 

490  CmLL  EXEC-:  1  ,  I  CNMD ,  BUFF ,  ILEN,  IOP1  ,  I  OPS  > 

500  CONTINUE 
r 
C 

GENT  *  2 . +0EXT/BETH*+2 
IF-!  HOANGS  >  RETURN 
C 

C  RECOVER  51  hUD  S2  FROM  C-.+ ,  S- 

00  300  ,1  =  1  ,  NM 

Si<  J>  =»  0.5*CMPLN-  SNGL-  RSP-,  J  >*PSM<  0  >  .» ,  SNGL<  ISP*;  J)>> 

Ss<  O'  *  0 , 3*CNPLX-.  SNGL-  RSP-.  J  >-RSH<  0  > .» ,  SNGL-  1  CP-  J)-I  SM<  J  >  >  > 

S  1  <  MNP  1  -  J  '*0 . 5  t-CMK.X'.'  SNGL<  RGPS-.  J  .'  +  RSNS-;  ..I  >  .'> ,  SNGL<  ISPS<  <J  .>+ISMS<  J  )  .)  > 
800  S2<  NMP  1  -  J  >*  0  . 5  *-C MF'LXX  SNGL-:  PSPS-  J>-RSMs-,  J  >  > ,  SMGL<  ISPS-:  J  >- 1 3HS<  *#>.»> 
C 

RETURN 

ENO 

SUBROUTINE  DCDV-'N,P1  , R2,3,G  > 

C  PERFORMS  COMPLE:;  MULTIPLICATION  8,  DIVISION  IN  DOUBLE  PREC 

DOUBLE  PRECISION  R 1  ,  R2  ,  S-:  2  > ,  0-'  2  ' ,  0 1 V 
IF  <  R  2  ,  E  0 , 0  .  .OHO.  S-  2  >  .  EG  0.  >  GO  TO  10 
DIV  *>  S-:  1  >-t  *2+S*' 2  >  *  +  2 
IF  (N.i  1,2 

1  C-;  1  .i  ~  OSG'-T*  -  Pi  *  *2+P2**2  >,DIV  T 
RETURN 

£  G-  1  j  »  <  R 1  *  >•  1  >  *-r  ♦  S',  2  »  ),  [i  I V 

I-1-'  2  *  *  -:  R2  *  s-  t  -  S’  t  *'S-  2  - c- :  v 

RETURN 

10  Q<  1  >  =  R 1  /S-,  1  .1 
CK  2  :»  =  0  . 

RETURN 

ENO 


"***r*mmmmm 


T a 00 0  04  1C  ON  CP 0  0 023  INC  00013  BLKS  P»<?000 


c 


£000 

2  0  0 1 
2  002 


'TH4  ♦OW 1  C.HT  ,  23  OULV 

PROGRAM  CP  IT 

CIMEM-ICH  PBUF '•'  3  > .  !C-CB<  H4  '.HuME*.  - 
rLl  H,M3,IP.I3  M3  1HTEH 
DOUBLE  PPEC  I  I  ON  N  ,  T  .  THC 

ssR  1  ‘e« ,t“. %i- \ - ;».«» > . «*.*< 

„  ,  |-,gi  Ipv  4  » ,  '.«;*>  *  . '  peuf  ‘ .5  ’  •  -•*•' 3  ’ 

OrtTfl  I-i-C/O.'  iOPTrl.  M.'  .  1L.  1  0.' 

LU  IN  =  LOuLU1  I  DU  Mil  •’ 

WRITE'  . . .  _  ^  T  h(1,tlE"  > 

FORMu  l  •  "ENTER  K  I  .  L-.fc  •  n  ,  F  •  . --.RT.LP 

REhDvLUIH,  f  '  PI  ,  BE  Th  ,  PTS  , ':'TRT  ,  r-TuHu 

WRITE-.  LU  IN.  2  0m  i  ■  ,  ,-H.:  »•, 

FORMAT <  3h.':  .12' 

PI  a  .?  ,  1  4  l'5  2  6 5  3  5  3  C>  0 
p  4  a  p  1  (•'  4  . 

$02  a  '-.QRTt.  2.  ' 

M  a  1 , 00/Rl  .  r  . 

THC  3  PATrtM-;  t  .  C'O.'C'SORTi  M++2-1  Do,'  ' 

TH2  ~  THE  - U  •  u£  2 
THI  *  THC--0 . 3E~2 
TH4  a  THC-l.?e-3 
PHC  3  P  I  -  2  .  *-THC 
H'i  a  HMI 

A  =  BE  r  R.O  2  DO' PI  ' 
hK  a  C-CiPT-  h  -  C  03-  SNijL-'  THC  -  1  » 

IH  a  $  TRY  f THC 
c  •-  ■  THC  -  TH  i.  PT 

,•  -i  -  c  ■  •  J 

C.  -•  '  ' 

b°OPEUeC-MTH  TILE  DESIHGATED  e. 

CHLL  OP EM-,  IOCB  •  IEPP  .  Nwl-lE  ,  l'>  TH  ■ 

IF-  IGP’R  ■  ET  .  <T  >  00  TO  333 
rBUF',  I  >  *  BETh 
PB'.iPC  2  •  *  1  ■  -  M 

p.euF’  3  >  r-  , .  ,  .Tpi)  •  »|  pemt  ecuivhi  f-uce  sthtememt 

'•.TC'P  FUME  T  1 1 DC  -  Uf'i-  .  «n 
CHLL  MR  l  TP'  Ii'CB  lEPP  ,  PBUP  '  1  L  1 
I  0  '"'.T  "a  DO  IIP  Cn-.LE*  TP  •  ' 

-  •  DmThU-  HfjT.'i-  PT-  ,.do-'N^t  >•*•♦4... 

R-TH 

PH  =  .T’.K' 

\  T  a  CUB’  TH  ' 

PD  a  PH  IH  DECREE'- 
Pf>  -  PH  *37 

°  P,-,MTUTELVhh"dE  HULST  DIVERGENCE  FhCTOP 
.  ..  T  •  '.T  '•  I  .  ■'*'"■  1 H 1  PM  '  ' 

J  I  “ 

=  p  *  TH 

FP  r  TmM'  D  ./ThM-  ' 

;=\3,‘. Ij!-'  ^TfoVnuPoF  the  tfhUumtted  ,.p,’r(„-TrC  ’  lay 


R 
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IP  =  Cl  1  *  4  .  ♦ '  |,-FP»*i 
I  S  =  G M<4  .  1  .  -FS 

ETH  =  PHC-PH 
W  =  S  I  N<  ETh  yfHh, 

TH8  -  THC+ETH/2 

DELI  »  HTuN',  SORT'  S I NC  THE  ':•*  *2  - 1  .  ,'NS  >/C  uS<  THB  >  ) 

DELS  =  2 , *HThN< NS > THNC DEL  1  »> 

DEL  1  =  2  . +DEL1 

ETHC  =  2  ,  >EETh-*v  COS'  THB  >-C  T  ♦  COS-*;  R  >/N  > 

B  I  -EThC  *P4  ♦•C  EL  1 
B2»EThC+P4+DFL2 
S  1 1  "SORT*  I  • 
s  1 2*SQRT <  ir  :• 

CB  I  =  COS'  B 1  - 
CBS* COS.  t;2  • 

'SB  1  *S I N <  B  1  .> 

SB 2= SINf  B2  > 

GOTO  40 

THIS  CODE  USED  WHEN  PH-VHC 
30  ETHD  =  ETHD-DETh 
ETH  =  ETHD, "57 , 2353 
PH  =  PHC-ETh 
PD  =  PH*57 . 2S58 
W  *  SI  NO  ETH.)*HK 
S  1 1  *  0  , 

51  2=0, 

CB  1  =  1. 

CB2= 1  . 

SB  1=0, 

SB2=  ii , 

4  0  C  HLL  PRES'-  M ,  F S  ,  F C  '> 

FS  =  FS/'S02 
FC-F  C/SQ2 

S  |  F'  =  FC  »  I  +  F S  ♦  S C  »  -t  I  I  1 
S  2 R - FC  +  C  B2-*  F  S «  S. C , .  S  I  2 
S  1  I  =FS*C  h  1  -FC  *":•!:  » 

5 2  I =  F S  +  C B2-FC  >  SC  2 
COMPUTE  PHwSE  DIFFERENCE 

PHI  =  57 , 2353  *h  I  HM-.  S  !  I  /S  1  ft  '/ 

RH2=57 . 2  SOS'  ♦  hThIK  S2 1  .-'S2P  '» 

IFC  S  I  P  ,  LT  .  0  .  .)  FH1  =  1  SO  .  +  PH  I 
IF'  F H  I  .  G  T  ISO,  i  F  H  1  =PH  1  -'16  0  . 

I  F(  S2P  .  L.  T  ,  0  ,  >  P H 2  =  t  3  0  ,  +PH2 
IF0FM2.GT.  i  0  0.  ,i  PH2  =PH2-7k  O  , 

D I Fh  =  PH2 -PM  I 

IFCDIFH.GT,  ISO.  .*  D  t  Fw-P  I  Fh-38  0  . 

IF1:  DIFk.LT.  ■!  3  0  .  ’  [■  1  F  h~D  I  Fh  +  3t  0  , 

COMPUTE  TOTHL  INTENSITIES 
T  I  S  = '  S  I  ft  *  *  2  +  S  1  I  ♦  *■  2  > 

T  I  P  =  '  S2P  *  *2 ■*•12  I  '  *  2  ( 

I NTEH=HLOG T ■  5*\ TIG+TIP >  > 

I  F<  Hl’IODC  PD  ,  !  E  -4  < .  EO  ,  n  ,  .  i.'p  j  TE<  l  U I N  ,  1  0  0  n  >  F  D  r  I :  ,  T I P ,  D I  Fh 
0  FORMmT' c3 .4  2E 1 4 . S  F I ft . 4  > 

CHLL  UP I TF0  I  DC & >  I  ERR , PBUF ,  I L  > 

IFOH.LT  THC  >  GOTO  SO 
IF-:  PD  .  LT  .  3  I  pMG  >  GOTO  3m 
GOTO  .>33 

S0i  IF<  7H  ,  GE  .  T 1 1 '  •  E  - C 2 
I  F'  r  H  ,  GE  ,  TH3  •  L  -  L"  3 
IF0  OH  ,  GE  .  TH-I  '  E  F4 


th-th+e 

1F<  TH  ,  LT , THC  >  GOTO  10 
E7AD*0, 

BETA3  ,  05 
GOTO  30 

999  CALL  CLOSE*  1DCB> 

l.lRITEi'  LU1N,  2003  >  1ERR 
£003  FORMATS  "IEPR«,M4  > 

STOP 
El  10 

SUBROUT  I  HE  FRES<  WE ,  P'S ,  PC  ) 

THIS  SUBR  CALCULATES  THE  FRESNEL  INTEGRAL  FROM  0  TO  WE 
W  »  MBS'!  WE  > 

F  a  i  i  ,  +  ,  92E*W  £  .  +  1  .  7‘S£^W+3 . 1  0  4 -*U  > 

G  ■  )  .  A‘2,+4,  t 42*W+3 . 492*W**'2+6 . S70*W**3  > 

A  a  1 , 57 079S*W**2 
C  «  COST  A  > 

S  -  S  I  N<  A  > 

FC  =  , 5+F^S-G^C 
FS  =  ,5-F*C-G*S 
IF<  WE  .  LT  .  0  ,  .i  GOTO  3 
FC  *  FC+.5 
PS  3  F S  +  .  9 
RETURN 

3  FC  3  ,5-FC 
FS  3  ,5-FS 
RETURN 
END 
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Glory  in  the  optical  backscatteriru)  from  air  bubbles  (U.  S.  Lamjley 
and  P.  L.  Marston)  Accepted  ior  Publication  in  Physical  Review 
Letters 


Abstract 

Observations;  of  light  backseat  to  rod  from  air  bubbles  I  n  a  viscous 
liquid  demonstrate  an  enhancement  due  to  axial  focus  in,;.  A  physical-apt  ies 
approximation  fo’*  the  cross-polarized  scattering  correctly  describes  the 
spacing  of  regular  features  observed.  The.  non-cross-polarized  scattering  is 
not  adequately  described  by  a  single  class  of  rays. 
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The  Mle  solution  for  electromagnetic  scattering  by  a  sphere  frequently 
does  not  lead  to  direct  interpretation  of  the  angular  scattering  pattern.  Con¬ 
sequently,  models  have  been  developed  to  facilitate  an  understanding  of  the 
structure  in  the  scattered  intensity  present  where  intensity  is  plotted  as  a 
function  of  the  scattering  angle  <j>  or  the  size  parameter  x  =  ko.  (k  =  wave- 

number;  a  =  sphere  radius).  These  models  have  emphasized  the  angular  regions 

2  3 

where  diffraction  is  important  for  a  drop  of  water  in  air:  the  rainbow,  ’ 

<fi  =  180°,^  and  <p  ~  0°.^’^  in  the  scattering  of  light  by  a  spherical  air 

bubble  in  a  liquid  or  in  glass,  the  real  part  of  the  refractive  index  of  the 

sphere  is  less  than  that  of  the  surroundings  and  the  models  must  be  signifi- 

7  8 

cantly  modified.  New  phenomena  appear,  such  as  diffraction  ’  in  the 

region  of  the  critical  scattering  angle  ■]>  .  Here  we  report  the  first  detailed 

observations  of  backsc.attering  by  air  bubbles  in  liquids  and  give  a  model 

which  describes  some  of  the  observed  features.  We  refer  to  this  as  glory 

3-5 

because,  as  in  the  case  of  drops,  the  ij>  180°  scattering  is  enhanced  when 
x  is  large. 

3  4 

Van  de  Hulst  ’  gave  a  partial  explanation  of  the  enhancement  for  crops 
by  noting  the  axial  focussing  of  those  backsce ttered  rays  which  have  a  non¬ 
zero  impact  parameter.  Winn  modeling  this  focussing  in  the  far  field,  dif¬ 
fraction  provides  an  essential  correction  to  ray  optics  because  the  facior  in 
the  scattered  intensity  which  accounts  for  geometrical  divergence  of  the  rays 

goes  to  00  as  <J)  -►  180°.  Examination  of  this  factor  in  ray-optics  models  of 

9 

scattering  by  bubbles  shows  that  this  «>  is  not  restricted  to  drops.  We  have 
modeled  the  backscattering  with  a  physical-optics  approximation.  The  proce¬ 
dure  is  to  (a)  compute  amplitudes  in  an  exit  plane  in  contact  with  the  bubble 
via  ray  optics,  and  (b)  allow  this  wave  to  diffract  to  the  far  field  where  the 
distance  from  the  bubble’s  center  R  >>  kci‘\ 
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Fig.  1  illustrates  several  rays  which  lead  to  backscattering .  The  paths 

are  determined  by  the  number  of  chords  p  and  m  *  m./m  where  the  refractive 

i  o 

indices  of  the  inner  and  outer  media,  in  and  m  ,  are  taken  to  be  real.  Fig.  1 

i  ( i 

is  drawn  with  m  ^  =  1.403  which  corresponds  to  an  air  bubble  in  the  dimethyl- 
siloxanc-polyiner  liquid  used  in  the  experiment.  All  rays  satisfy  sin0  =  msinp. 
For  <f>  “  180°,  tiie  off-axis  (or  glory)  rays  have  u  ~  t)  and  p  =  p where: ^  0  = 
pp  +  (2g+2-p)90°,  g  is  a  non-negative  integer  (g  =  0  for  rays  in  Fig.  1)  and 
m  <  .1  requires  that  p  >  3.  The  exit  plane  (dashed  line  in  Fig.  1)  touches  C 
with  its  normal  parallel  to  the  propagation  direction  of  the  incident  wave. 

Our  description  oi  the  Held  in  the  ex  i  t  plane  is  iacililuted  by  consid¬ 
ering  the  propagation  ol  a  wavelet  do  which  lies  close  to  the  hacksc.attered 
path.  Fig,  1  shows  de  for  p  e  3;  it  emerges  as  curve  dV.  This  curve  appears 
to  come  from  a  ring, -like  source  at  F  known  as  the  local  circle  in  the  analo¬ 
gous  p  =  2  scattering  from  drops*  with  /?  m  •  2.  Thu  source  is  ring-like 
because  the  figure  may  be  rotated  around  the  '  axis.  The  radius  of  the  ring 
is  b  =  rtslnO.  After  the  incident  ray  crosses  t  tie  dashed  vertical  plane  (the 
entrance  plane),  the  p  r  op  ijj . » rjlo  r  i  phase  delay  for  reaching  the  exit  plane  is 
H  ■  ka[l  -  cost)  +  (1  -  c-osfi) sec  (0-P)  +  2mpeosp|  .  The  ray  crosses  the  exit 


plane  at  a  radius  a  from  C'  with  s/a  =■  sinl'  -  (1  -  cosP)  t  an  (3-6)  .  The 

radius  n  of  arc  d'e*  follows  from  the  curvature  at  s  -  b:  r*.  =»  k(d  i|/ds*")  ~ 

a[l  +  *.! (px—  1 )  ^cosU]  where  i  -  tanp/tanO.  Tin*  spreading  of  the  wavelet  Is 

characterized  by  q  ••=  11m  d,er/de  as  de  ►  0  where  tiie  bar  denotes  the  arc  length 

An  equivalent  expression  for  q  is  ]  1. i m [ b  -  s(0))/(b  -  asin(J)|  as  0  >  0; 


Its  value  from  [.'Hospital's  rule  is  a/  (a  --  a).  Vectors  (9.  *  1,2)  denote 

orthogonal  hasis  vectors  in  both  the  entrance  and  exit  planes;  e  is  chosen 
parallel  to  the  polarization  of  the  incident  wave's  electric  field  E^exp  (-iurt.) 

O 

Tn  tiie  exit  plane,  the  (  l  eld  E^e^  of  t  ue  outgoing  ptlr  glory  wave  is 
computed  by  applying  Van  de  Hu  1st’ a  method  of  first  decomposing  the  fields 


4 


perpendicular  and  parallel  to  the  scattering  plane.  ’  Exit-plane  polar 
coordinates  centered  on  C'  are  (s,i^)  where  is  the  angle  relative  to 
and  s  and  4)  denote  local  basis  vectors.  We  assume  x  >>  1  and  use 
Fresnel's  coefficients  r^  for  the  internal  reflections  where  j  ■  1,2  for 
fields  parallel  to  i|'  and  s,  respectively.  If  )  s— b  |  <<  CL,  the  multiple 
internal  reflections  give: 

F.p  -  Ejq-'3  F^explin  +  ik(s  -  b)2/2u)  (1) 

where  n  ■  p  +  n(p  =  p),  F  (t|0  =  c^sin>  +  c^cos  'll,  F  (i|>)  =  ~  c^)sin2i^, 

and  Cj  ■  (-1)^^  ^rj  “  r^)  •  The  new  phase  term  p  accounts  for  the 
crossing  of  causLics  or  "focal  lines";  its  value  is^’^  -ir(p+g)/2.  The  r^ 
are  evaluated  at  8:  r^  *  sin(0  -  p)/sin(8  +  p)  *  “  tan(0  -  p)/tan(0  +  p) . 

The  sign  factor  in  Cj  accounts  for  a  geometrical  inversion  (present  when 

j  =  2  and  p  is  odd)  which  is  not  evident  in  descriptions  of  p  =  2  glory 

.  a  3,4 

in  drops. 

£ 

The  field  E^  at  a  distant  point  Q  is  computed  as  follows.  The  left  exten¬ 
sion  of  the  CC'  axis  makes  an  angle  y  with  C'Q.  When  y  is  small  and  C'Q  = 

2  12 
R'»ka  ,  scalar  diffraction  theory  and  the  Fraunhofer  approximation  give: 


us,  *,m'  * r,) 

2tf  iR '  q  J 


.,£  ik(s-b)  /  2ct , 
>W  e  ds 


f£  e-iks  sinycos(^-f;)d(jj 


where  £  is  the  angle  between  e^  and  the  projection  of  C'Q  on  the  exit  plane. 
In  Eq.  (2),  the  approximation  given  by  Eq .  (1)  has  been  extended  beyond  its 
useful  domain  in  anticipation  of  the  stationary  phase  approximation  (SPA)  of 
the  integral.  Direct  evaluation  of  Eq.  (3)  gives  W  (Y,0  =  W  (y,f,,s  =  b) 
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“  TT[(c1  +  c^)J^(u)  +  (CjL  -  c  ^ )  J ( »i )  cos2fj  and  W''  -  11  (e^  -  c^)  (u)  sln2^  where 
u  ■»  kbsiny.  The  SPA  of  Kq .  (2)  Rives  the  pth  glory  contribution  to  the  scat- 
tered  field  when  kb  /it,  and  thus  x,  are  large.  In  the  experiments  to  be 
described  x  .>  A000  and  the  SPA  is  applicable. 

I 

The  total  field  may  be  approximated  bv  summing  the  K  from  Kq .  (2) 

with  the  fields  due  to  axial  reflections  and  surface  waves.  Surface  wave 

contributions  shouLd  be  small  for  the  observed  bubbles  due  Lo  the  largeness 

of  x.  To  determine  which  glory  and  axial  terms  are  important  to  the  total 

field,  and  for  other  heuristic  reasons,  consider  the  Ji.- polarized  intensity 
t 

1  of  the  pth  field  taken  alone.  The  SPA  ol  Kq .  (2)  gives: 

P 

1.^  -  (2/ n)xl  f  1 WV*  ( ,  ,.<  )  ]2  (A) 

2  2 

where  I  -•  l  a  / AR  is  tlie  total  intensity  at.  a  distance  R  =  CQ  from  a 
K  1 

9 

perfectly  ref loctjnjj  sphere  of  radius  a  predieted  by  ray  optics,'  1  is 

2  t  2  j 

the  incident  intensity,  and  I  b "  u/rt  q  b  (a  -  ct>/u  .  in  Kq.  (A),  K 

P  *  K 

has  replaced  R1  from  (2)  and  y  becomes  180°  -  J;  because  R  it. 

3,9 

Geometrical  optics  gives.  Lite  intensities  I  of  separate  axial  reflections 

P 

> 

(e.g.  p  =  0  and  2  in  rig.  1)  wb ;  h  ate  propottion.il  to  a".  The  strongest. 

reflection  Iras  p  ~  0  and  1’.  77  1;  for  y  -  0,  I*  :  I  (in- 1  V  /(mil)*'  wit  lie 

0  K 

«2  li,  ( 

1_  “  0.  Since  f  does  not  depend  on  a,  I  "  lot  and  glorv  terms  dominate 

0  P ,  g  p 

the  backseat tering  when  a  is  large. 

Cons  idet  a  bubble  with  x  AOOD  and  m  r-  l.AOi  '  .  The  strongest  glory 

terms  have  g  -  0  and  p  77  3,  A,  and  3;  the  I  V I  for  y  -  0  are  respectively 

1.03,  0.A3,  and  0.16.  The  1  ^  decrease  with  increasing  p  due  to  the  partial 

reflections  in  the  bubble.  The  strongest  axial  ray  gives  I^/l  =  0.028.  The 

(J  R 

interference  ot  the  fields  depends  on  a  and  our  Hie  computations  verity  that 


the  backseat t ered  intensity  is  not  simply  proper t Iona l  to  it  even  for  this 
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large,  value  of  x.  The  2  “  2  (cross-polarized)  scattering  is,  however, 

nearly  dominated  by  the  p  =  3  glory  term.  Due  to  symmetry,  It  -  2  scattering 

2 

vanishes  as  y  -*■  C.  The  I  (y  ^  0,£.)  have  maxima  at  f,  =  *’45°  and  4135°  and 

P 

they  vanish  at  f,  -  0°,  +90°,  and  180°.  Let  y  =  j  locate  the  first 

P 

2  2 
maxima  of  I^(y,C  “  45°).  The  largest  2  ~  2  terms  have  I  (y  ,  (,  ■=  45°)/I^  = 

0.53  and  0.10  for  p  =  3  and  4.  To  the  extent  that  p  4  '3  scattering  may  be 

neglected,  the  2=2  intensity  will  be  quasi-per indie  In  v. 

We  have  numerically  verified  the  validity  of  Eq .  (4)  by  using  Debye's 
3  4 

localization  principle  ’  to  modify  Mio  theory  so  that  only  partial,  waves 
associated  with  p  =  3  rays  were  included  in  the  Mie  series.  Furthermore,  when 
Eq.  (4)  is  applied  to  spheres  with  certain  m  >  1,  the  resulting  l^(y  =  0) 

p 

agree  with  the  glory  "analog"  tabulated  in  Ref.  11.  This  analog  was  derived 

by  applying  the  Watson  transformation  to  the  y  =  0  Mie  series. 

Fig.  2  diagrams  the  experiment.  A  syringe  injected  bubbles  into  the 

liquid.  The  liquid  had  a  high  kinematic  viscosity  (~600  000  cS)  and  a  single 

bubble  could  be  observed  for  hours.  The  laser's  power  output  was  5  mW  and  the 

beam  diam  was  5  mm.  The  wavelength  in  the  liquid  2n/k  was  6  32.8  nm/.l.403; 

e^  lay  in  the  splitter's  plane  of  incidence.  The  camera  was  iocused  on  m  so 

7  1  2 

the  photographs  recorded  the  far-field  intensity  pattern.  ’  Photographs 
were  made  with  a  -  0.3-0. 8  mm  corresponding  to  x  -  4000-11000.  Exposure  times 
were  typically  5s  for  TrlX  film  and  a  200  nun  focal  length  camera  lens. 

Fig.  3  demonstrates  that  the  scattering  has  roughly  the  dependence  on 
£  predicted  by  Eq.  (4);  £,  =  0°  corresponds  to  scattering  toward  the  top  of 
the  photographs  and  y  =  0°  corresponds  to  the  center  of  the  symmetry.  Fig.  3(b) 
shows  that  the  9.  -  1  scattering  for  y  >  0.2°  is  significantly  stronger  for 
t,  =  +90°  than  it  is  for  f,  =  0  .  This  agrees  with  the  following  modeL  resulLs: 
(i)  (Cj/c2)2  >>  1  (for  p  =  3  we.  predict  c^/c^  -5.2);  and  (ill  for  this  x. 
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vi  1 

the  1^  depend  only  weakly  on  £  and  are  dominated  by  the  I  .  One  prediction 

of  Eq.  (A)  could  be  quantitatively  checked:  when  both  siny  =  y  and  u  >>  1, 

2 

the  minima  in  I  should  be  spaced  by  Ay  rad  such  that  kbAy  -  tr  where  for 

p  a  3,  b/a  =  0.AA7.  Fig.  A  compares  this  with  the  mean  spacing  of  ~A0  dark 

rings  lying  outside  the  9th  ring  from  the  center.  The  error  bars  combine 

uncertainties  in  measured  a  and  Ay  with  those  of  corrections  due  to  refraction 

at  the  cell-air  interface^  and  the  tilt  of  the  cell.  Fig.  A  shows  that  p  =  3 

rays  dominate  the  1-2  scattering.  The  modulations  of  the  intensity  along 

£  =  ±A5°  in  Fig.  3(b)  show  that  otlier  rays  contribute  to  f  =  1  scattering 

1  2 

since  the  predicted  I  «  [J  (u))  . 

In  conclusion,  backseat tering  from  bubbles  can  be  enhanced  by  axial 

focusing.  The  number  of  significant  glory  terms  depends  on  m.  The  main 

3 

contributions  differ  from  those  for  water  drops  where  surface  waves  and 

other  diffraction  related  terms'*  play  an  essential  role.  If  focusing  were 

7-9 

not  present,  scattering  by  large  bubbles  would  be  <<  1  in  the  region 

K 

(<f>c  +  10°)  <<{)<.  180°  where  <}>  =  2  cos  ^m  =  89°  for  m  ^  =  1.A03.  We  also 

find  evidence  of  p  =  3  glory  in  Mie  computations  for  bubbles  in  water . 

This  work  was  supported  by  the  Office  of  Naval  Research.  P.  L.  Marston 
is  an  Alfred  P.  Sloan  Research  Fellow. 
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Fig.  Captions 

Fig.  1.  Rays  which  contribute  to  hackscattering.  The  local  angle  of  inci¬ 
dence  is  0  and  C  is  the  bubble's  center. 

Fig.  2.  Apparatus  for  observing  backscattering  from  bubbles. 

Fig.  3.  Photographs  for:  (a)  crossed  polarizer  (£,  =  2  scattering); 

(b)  uncrossed  polarize  (Z  ■  1);  and  (c)  no  polarizer.  The  incident  polar¬ 
ization  was  vertical,  a  °  0.49  mm  and  x  =  6830. 

Fig.  4.  Measurement  and  model  for  the  angul.,r  separation  of  the  dark  rings 
in  the  Z  -  2  scattering. 
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